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I.  SUMMARY 


A  variety  of  hydrocarbons  has  been  pyrolyzed  in  a  single -pulse  shock  tube 
over  the  temperature  range  of  1000  to  2400  K  and  for  dwell  times  of  500  to  700 
microseconds.  Gas  samples  of  reactant,  intermediate,  and  final  products  were 
collected  and  analyzed  using  gas  chromatography.  Experimental  data  were  used 
in  conjunction  with  a  computer  model  to  develop  (or  confirm  existing)  detailed 
chemical  kinetic  models  for  several  of  the  hydrocarbons  that  were  pyrolyzed. 
Model  results  agreed  well  with  experimental  data  not  only  for  the 
decomposition  of  the  parent  compound  and  formation  of  products.  A  major  focus 
of  the  modeling  effort  was  to  describe  accurately  the  formation  of  aromatic 
species.  Specific  mechanisms  and  rate  coefficients  have  been  proposed  in 
order  to  describe  the  experimental  data.  These  latter  processes  are  believed 
to  lead  to  the  formation  of  polycyclic  aromatic  hydrocarbons,  soot  precursors, 
and  eventually  soot.  The  work  also  has  focused  on  the  pyrolysis  of  acetylene, 
since  this  species  appears  to  be  the  dominant  species  in  soot-forming  regions 
of  flames  and  has  been  shown  to  play  a  major  role  in  the  soot-forming  process. 


II.  INTRODUCTION 


The  production  of  soot  in  diffusion  flames,  such  as  that  existing  in  gas 
turbine  engines,  is  directly  related  to  the  chemical  nature  of  the  fuel.  In 
addition,  research  has  shown  recently  that  chemical  kinetic  processes  control 
the  growth  of  polynuclear  aromatics  and  probably  the  formation  of  soot 
precursors.  The  production  rate  of  soot  precursors,  in  turn,  limits  the 
production  rate  of  the  number  of  soot  particles  and  hence  effects  the  total 
mass  of  soot  produced.  This  work  focuses  on  obtaining  experimental  data  on 
the  pyrolysis  and  oxidative  pyrolysis  of  hydrocarbons  under  soot- forming 
conditions.  A  single-pulse  shock  tube  has  been  selected  for  the  experimental 
tool  since  it  is  one  of  the  few  devices  capable  of  obtaining  detailed 
information  on  product  distribution  over  the  range  of  1200  to  2400  K  in  an 
environment  free  from  complications  due  to  diffusional  effects.  Using  the 
experimental  data  and  single -pulse  information  obtained  from  other  work, 
detailed  chemical  kinetic  mechanisms  describing  the  pyrolysis  of  a  variety  of 
hydrocarbons  are  being  developed.  Chemical  mechanisms  describing  the 


formation  of  aromatic  hydrocarbons  (from  aliphatic  hydrocarbons)  are  included 
in  the  model. 


III.  RESULTS 
(A)  Experimental  Effort 

Seventeen  separate  series  of  single-pulse  shock  tube  experiments  were 
completed  during  the  first  year  of  this  contract  using  both  contract  (AFOSR) 
and  Corporate  (UTRC)  funds  (see  Table  I).  During  the  second  year,  ten  series 
of  experiments  were  completed  under  contract  sponsorship  and  eight  under 
Corporate  sponsorship.  (See  Table  II).  During  the  final  year  of  the  program, 
seven  new  series  of  experiments  were  performed  under  the  AFOSR  contract  and 
six  under  corporate  sponsorship.  Initial  reactant  concentrations  for  these 
thirteen  experiments  are  listed  in  Table  III.  Each  series  represents  a 
collection  of  individual  shock  tube  experiments  along  with  chemical  analysis 
of  reactant  and  products  (via  gas  chromatography)  for  a  given  initial  reactant 
condition  (species,  initial  concentration).  Shock  strengths  were  varied  to 
produce  a  range  in  initial  post-shock  temperature  from  approximately  1100  to 
2400  K. 

Examples  of  the  type  of  experimental  data  obtained  from  the  shock  tube  are 
shown  in  Figs,  la- 3a  for  the  pyrolysis  of  ethene  and  in  Figs.  4-8  for  a 
fuel-rich  oxidation  (-2%  oxygen)  with  about  3  1/2  percent  ethene  in  argon. 
Preliminary  modeling  results  for  the  pyrolysis  are  included  in  Figs.  lb-3b, 
using  a  chemical  kinetics  sequence  slightly  modified  from  that  presented  in 
Appendix  D.  The  formation  of  aromatic  species  varies  substantially  with 
oxygen  concentration.  In  Figs.  9-11,  the  formation  of  benzene, 
phenylacetylene  and  styrene  are  shown  for  three  different  oxygen 
concentrations  (2,  0.5  and  0%).  It  is  clear  that  as  oxygen  concentration 
increases  from  zero,  the  production  of  aromatics  occurs  at  lower  and  lower 
temperatures.  Other  Interesting  features  during  oxidation  include  high 
formation  rates  of  methane  and  a  high  CO/Cj^  ratio.  The  presence  of 
methane  indicates  the  importance  of  methyl  radicals  during  oxidation  whereas 
the  high  00/0^2  rati°  substantiates  suggestions  that  ethene  and/or  vinyl 
radicals  decompose  via  addition  of  oxidative  radicals  (0  or  OH)  rather  than 
extraction  of  H-atom  by  such  radicals. 

Other  experiments  include  pyrolysis  of  styrene  (an  alkylated  aromatic)  and 

1.3.5- hexatriene.  Styrene  was  selected  due  to  its  poential  role  in  the  growth 
from  single  to  double  ring  aromatics.  The  linear  six-membered  aliphatic, 

1.3.5- hexatriene,  was  selected  since  it  is  the  parent  of  two  radical 
intermediates,  i.e.,  1 , 3 , 5-hexatrien-l-yl  and  1 , 3 , 5-hexatrien- 3-yl  which  have 
been  postulated  as  intermediates  in  the  production  of  benzene  from  alipatics. 
As  much  as  55%  of  the  initial  hexatriene  is  converted  to  benzene  at  low 
temperatures  (1300-1400K) ,  thus  confirming  that  a  rapid  conversion  from  the 
linear  to  aromatic  compound  will  occur.  There  is  some  difficulty  in 
interpretation  of  the  experiments  on  styrene  due  to  a  15  percent  uncertainty 
in  the  ability  to  accurately  measure  styrene.  Below  1400  K,  there  is 
relatively  little  decomposition  of  the  aromatic  ring  although  styrene  does 
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decompose  principally  to  benzene,  phenylacetylene ,  hydrogen,  ethene  and 
acetylene.  At  these  temperatures,  there  is  a  slight  preference  to  form 
phenylacetylene  but  the  trend  is  opposite  at  higher  temperatures.  At  1450  K 
and  above,  there  is  rapid  loss  of  the  initial  styrene  although  there  is  not  a 
compensorate  production  of  low  molecular  weight  products.  Presumably  the 
styrene  is  converted  to  unobservable  high  molecular  weight  species. 

Mechanisms  describing  this  phenomena  are  under  investigation  and  they  may  be 
related  to  the  formation  of  soot  precursors. 

Related  to  the  experimental  work  is  a  computer  code  which  has  been  written 
to  help  organize,  document,  and  plot  experimental  data.  During  the  course  of 
this  research,  the  analytical  capabilities  were  upgraded  (under  corporate 
sponsorship)  in  order  to  obtain  analysis  of  not  only  the  to 
Cio -hydrocarbons ,  hydrogen  and  oxygen,  but  also  the  analysis  of  the  carbon 
oxides  (CO  and  CC^)  at  parts-per-million  concentration  levels  using 
catalytic  hydrogenation.  In  addition,  the  ability  to  quantitatively  analyze 
methane  and  the  C_ -hydrocarbons  has  been  enhanced  by  a  duplication  in 
measurements;  that  is,  these  species  are  separated  and  measured  on  capillary 
and  packed  columns . 


(B)  Model  Effort 


(1)  Code  Development 

Enhancement  to  the  computer  codes  for  development  of  models  used  in  this 
contract  have  been  made  under  a  Corporate -sponsored  program.  Yet  the 
advancements  enhanced  the  progress  in  the  present  contract  work  and  therefore 
will  be  described  briefly. 

Chemical  kinetic  modeling  is  performed  using  CHEMKIN,  LSODE,  CHEMSEN  and  a 
shock  tube  code  for  CHEMKIN  originally  written  by  personnel  at  Sandia  National 
Laboratories.  As  described  in  the  proposal  for  this  contract,  the  original 
shock  tube  code  had  already  been  modified  in-house  in  order  to  model  the 
quenching  process  in  a  single-pulse  shock  tube.  Plotting  packages  have 
been  developed  under  corporate -sponsorship  to  facilitate  interpreting  data 
from  CHEMKIN  and  CHEMSEN.  Plots  of 

1.  Species  concentrations 

2.  Reaction  contribution  to  the  formation/destruction  of  each  species 

3.  Principal  initiation  and  termination  reactions  and  their 
respective  rates 

4.  Reactions  which  are  (nearly)  equilibrated. 

5.  Sensitivity  coefficients 
can  now  readily  be  obtained. 


These  new  capabilities  may  not  appear  as  highlights  in  a  publication 
or  a  presentation,  but  have  been  of  great  value  in  meeting  the  overall 
research  objectives  of  this  contract. 

(2)  Chemical  Kinetic  Modeling 

During  the  course  of  this  program,  detailed  chemical  kinetic  modeling  has 
been  performed  for  the  pyrolysis  of  acetylene,  ethylene,  vinylacetylene ,  and 
benzene.  In  addition,  the  pyrolysis  of  biacetyl,  acetaldehyde,  and  pyridine 
has  been  performed  under  corporate  sponsorship.  To  perform  this  work,  it  was 
necessary  to  include  in  the  modeling  a  great  number  of  other  species, 
including  radicals  (see  Tables  VI  and  II  of  Appendices  A  and  C,  respectively). 
Throughout  this  effort,  there  was  a  strong  attempt  to  ensure  that  models 
(mechanisms  and  kinetics)  describing  the  growth  to  higher  molecular  weight 
species  were  consistent  via  thermochemistry  with  models  describing  the 
decomposition  of  these  larger  molecules.  The  constraints  of  describing  the 
decomposition  of  the  different  reactants  while  simultaneously  predicting  the 
formation  of  many  products  was  very  severe.  Consequently,  we  have  a  lot  of 
confidence  regarding  the  general  validity  and  rate  constant  expressions 
determined  in  this  study.  A  depiction  of  the  model  showing  formation  and 
destruction  of  a  single  ring  aromatic  is  shown  in  Table  IV  and  specific 
reactions  and  rate  constants  are  in  Table  V. 

A  focus  of  this  three -year  program  was  the  kinetics  and  mechanisms  of 
acetylene  pyrolysis.  The  importance  of  this  work  was  two-fold.  First, 
acetylene  has  been  identified  as  the  dominant  hydrocarbon  species  present  in 
soot- forming  regions  of  flames  as  well  as  suggested  to  be  the  main  growth 
species  of  soot  precursors  and  soot  itself.  The  chemistry  of  acetylene 
pyrolysis  clearly  plays  a  major  role  in  the  growth  of  soot  precursors. 
Secondly,  recent  published  articles  (see  Refs.  6  and  7  of  Appendix  A)  have 
proposed  that  acetylene  pyrolyzes  not  by  a  chain  radical  involving  H-atoms  as 
a  chain  carrier  but  rather  via  vinylidene,  a  diradical,  as  an  intermediate. 
These  proposal  could  lead  to  a  relatively  simple  model  describing  the 
formation  of  soot  precursors  (and  eventually,  soot) .  Unfortunately  if  the 
proposal  of  vinylidene  is  inaccurate,  the  use  of  such  mechanisms  in  a  soot 
model  could  lead  to  serious  errors  in  predictions  with  slight  changes  in  the 
environment . 

These  arguments  led  to  a  major  study  on  acetylene  pyrolysis  incorporating 
data  obtained  from  low  temperature  studies  in  flow  reactors.  The  modeling  of 
acetylene  pyrolysis  has  been  completed  (described  in  Appendix  A)  and  the  major 
conclusion  from  ths  analysis  has  been  that  acetylene  pyrolyzes  indeed  via  a 
radical,  chain  mechanism  which  is  initiated  by  the  presence  of  impurities  of 
acetone.  The  acetone  is  a  stabilizer  added  to  acetylene  and  is  difficult  to 
remove  below  levels  of  1  part  in  1000.  An  early  version  of  the  manuscript  did 
not  include  the  contribution  of  acetone  and  was  weakened  by  a  requirement  for 
an  unusually  high  A- factor  for  the  initiation  rate. 

In  a  second  paper  (See  Appendix  C)  on  acetylene  pyrolysis,  mechanisms  for 
growth  to  3-ringed  aromatics  were  analyzed.  By  using  rate  constants  and  a 
model  which  predicts  concentrations  of  low  molecular  weight  species 
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accurately,  basic  concepts  describing  ring  growth  were  found  to  be  similar  to 
that  described  previously,  although  some  specific  mechanisms  were  found  to  be 
substantially  different  from  those  proposed  previously.  As  described  in 
Appendix  C,  a  dominant  reaction  sequence  for  the  growth  from  one  to  two  rings 
is  shown  in  Table  VI.  At  the  present  time,  there  is  little  experimental  data 
which  supplies  evidence  to  confirm  the  existence  of  such  a  mechanism  and  the 
conclusion  is  drawn  solely  from  thermodynamic  and  chemical  kinetic 
considerations.  Despite  the  lack  of  data  with  which  model  calculations  can  be 
compared,  there  is  an  optimism  for  the  ability  to  predict  the  formation  of 
soot  precursors  and  the  formation  rate  of  incipient  particles.  Reasons  for 
this  optimism  are  described  in  Appendix  B.  A  plausible  kinetic  mechanism  for 
describing  these  processes  is  presented  in  Apendix  C. 

(C)  Model  of  Soot  Particle  Formation 

Most  practical  applications  of  combustion  processes  utilize  the  diffusion 
flame  concept  because  of  the  simplicity  and  broad  range  of  stability  provided 
by  this  combustion  mechanism.  Since  the  diffusion  flame  does  not  promote 
rapid  mixing  of  fuel  and  oxidizer,  however,  the  fuel  vapor  undergoes  high 
temperature  pyrolysis  before  it  is  oxidized  completely,  resulting  in  a 
tendency  for  the  intermediate  formation  of  soot  particulates.  Despite 
numerous  research  efforts  to  better  understand  soot  particle  formation  in 
combustion  processes  for  the  purpose  of  developing  soot  emission  control 
strategies,  soot  formation  in  practical  combustors  is  still  a  problem. 
Moreover,  theoretical  modeling  of  soot  particle  formation  within  combustion 
processes  has  not  been  very  successful  because  of  the  complexity  of  the  many 
chemical  and  physical  processes  involved. 

Since  we  believe  that  the  results  of  this  program  provide  a  basis  for 
understanding  of  the  pre -particle  chemistry  of  large  aromatic  species  which 
form  during  the  pyrolysis  of  fuel,  the  development  of  a  theoretical  model 
which  describes  soot  particle  formation  in  combustion  processes  has  been 
initiated  at  UTRC.  Model  development  will  be  completed  under  separate 
funding.  This  theoretical  model  will  simulate  the  evolution  of  a  spatially 
homogeneous  soot  particle  (aerosol)  as  a  function  of  the  time -dependent 
particle  size  disribution.  The  physical  phenomena  to  be  included  are:  (1) 
particle  nucleation  or  inception  from  a  time  varying  source  of  gas  phase 
chemical  species,  (2)  particle  growth  due  to  heterogeneous  condensation  of 
gas  phase  chemical  species  (often  called  gas-to-particle  conversion),  and  (3) 
coagulation  of  particles  via  particle  collisions.  By  judicious  use  of  the 
theory  developed  under  this  program  on  the  chemical  formation  of  pyrolysis 
products  from  hydrocarbon  fuels,  models  will  be  developed  for  use  in 
predicting  the  gas  phase  composition  and  the  rates'  of  particle  inception  and 
gas-to-paarticle  conversion  as  functions  of  time.  Numerical  solution  of  the 
equations  which  govern  soot  particle  formation  will  be  accomplished  by 
utilizing  an  existing  computer  code  named  MAEROS ,  a  general  multicomponent 
aerosol  computer  model  which  was  developed  by  F.  Gelbard  at  Sandia  National 
Laboratories,  Albuquerque,  New  Mexico.  A  copy  of  this  code  has  been 
obtained  by  UTRC  personnel  and  modifications  to  the  code  for  describing 
soot-formation  have  been  initiated. 


IV.  LIST  OF  PUBLICATIONS 


A  paper  entitled  "The  Pyrolysis  of  Acetylene  Initiated  by  Acetone"  by  M. 

B.  Colket,  H.  B.  Palmer,  and  D.  J.  Seery  has  been  submitted  for  publication  in 
Combustion  and  Flame.  (See  Appendix  A) 

A  paper  entitled  "The  Pyrolysis  of  Acelylene  and  Vinylacetylene  in  a 
Single-Pulse  Shock  Tube"  by  M.  B.  Colket  has  been  accepted  for  publication  in 
the  Twenty-First  Symposium  (International)  on  Combustion  (see  Appendix  D) . 

An  article  entitled  "Single-Pulse  Shock  Tube  Examination  of  Hydrocarbon 
Pyrolysis  and  Soot  Formation"  by  M.  B.  Colket  has  been  published  in  Shock 
Waves  and  Shock  Tubes,  Proceedings  of  the  Fifteenth  International  Symposium  on 
Shock  Wave  and  Shock  Tubes,  Edited  by  D.  Bershader  and  R.  Hanson,  p.  311, 
Stanford  University  Press,  Stanford,  California  1986  (see  Appendix  F) . 

Under  corporate -sponsorship,  a  manuscript  entitled  "Shock  Tube  Pyrolysis 
of  Pyridine  and  the  Thermochemistry  of  the  Pyridyls"  by  J .  C.  Mackie  and  M.  B. 
Colket  is  nearly  complete  and  will  be  submitted  to  the  Twenty- Second  Symposium 
(International)  on  Combustion. 


V.  MEETING  INTERACTIONS 

The  Principal  Investigator  of  this  program  attended  the  following 
meetings  during  this  program. 

1.  Eastern  Section  of  the  Combustion  Institute,  Gaithersburg,  MD, 
November  2-6,  1987.  M.  B.  Colket  presented  a  paper  entitled  "Some 
Thoughts  on  Pre-Particle  Chemistry".  (See  Appendix  B) 

2.  Division  of  Fuel  Chemistry,  New  Orleans  National  Meeting,  American 
Chemical  Society,  August  31-September  4,  1987.  M.  B.  Colket 
presented  a  paper  entitled  "Kinetic  Mechanism  for  Pyrolysis  of 
Acetylene  Near  1000K"  at  the  Symposium  on  Advances  in  Soot 
Formation  Chemistry.  (Appendix  C)  Also  see  ACS,  Division  of  Fuel 
Chemistry  preprints,  32  (3),  417  (1987). 

*3.  Department  of  Energy,  Combustion  Research  Contractors'  Meeting, 
Mills  College,  Oakland,  California.  June  2-5,  1987.  M.  Colket 
was  invited  by  Dr.  A.  Laufer  to  be  an  observer  and  participant  in 
discussions  at  the  D.O.E.  meeting.  (See  Appendix  D) 

4.  Twenty-first  Symposium  (International)  on  Combustion  in  Munich,  W. 
Germany,  August  1986.  A  paper  (See  Appendix  D)  entitled  "The 
Pyrolysis  of  Acetylene  and  Vinylacetylene  in  a  Single -Pulse  Shock 
Tube"  was  presented  by  M.  B.  Colket.  In  addition  M.  Colket  served 
as  Poster  Chairman  and  member  of  the  Program  Advisory  Committee, 
the  Program  Sub-committee  and  the  Publications  Committee. 
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*5.  The  1986  Spring  Technical  Meeting  of  the  Central  States 

Section  of  the  Combustion  Institute  in  Cleveland,  Ohio,  May  1986. 

A  paper  entitled  "Formation  of  C^ -hydrocarbons  and  Benzene  from 
Pyrolysis  of  Biacetyl"  was  presented  by  M.  B.  Colket. 

*6.  Department  of  Energy,  Basic  Sciences  Contractors  Meeting, 

Arlie  House,  Va,  May  1986.  M.  Colket  was  invited  by  Dr.  A.  Laufer 
to  be  an  observer  and  participant  in  discussions  at  the  D.O.E. 
meeting. 

7.  Division  of  Fuel  Chemistry,  New  York  City  National  Meeting, 
American  Chemical  Society,  April  13-16,  1986.  A  paper  entitled 
"Pyrolysis  of  C  H  ,"  by  M.  B.  Colket,  III,  was  presented. 

(Appendix  H) .  Also  see  ACS,  Division  of  Fuel  Chemistry  preprints, 
31  (2),  p.  98,  1986. 

8.  Eastern  Section  of  the  Combustion  Institute,  Philadelphia,  PA, 

Nov.  4-6,  1985.  A  paper  entitled  "Pyrolysis  of  Vinylacetylene , " 
by  M.  B.  Colket,  III  was  presented.  (Appendix  G) . 

*9 .  International  Conference  on  Chemical  Kinetics,  National  Bureau 
of  Standards,  Gaithersburg,  MD,  June  1985.  A  poster  paper 
entitled  "Pyrolysis  of  Acetylene  in  a  Single-Pulse  Shock  Tube,"  by 
M.  B.  Colket,  III.  was  presented. 

10.  Eastern  Section  of  the  Combustion  Institute,  Clearwater  Beach, 
Florida,  December  1984.  A  paper  entitled  "Addition  of  Alkyl 
Radicals  to  Aromatic  Rings,"  by  M.  B.  Colket,  III  was  presented. 
(Appendix  F) . 
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TABLE  I 


SERIES  OF  EXPERIMENTS  COMPLETED 
DURING  FIRST  YEAR  OF  PROGRAM 


Initial 

Reactant  Concentration  (% 


Benzene ,  C ,H, 

0  0 

0.12 

Benzene,  C.H, 

6  o 

0.012 

Vinylacetylene ,  C^H^ 

1.0 

Vinylacetylene ,  C^H^ 

0.11 

Vinylacetylene  C^H^ 

0.01 

C2VC6H6 

0.82/0.06 

C2H2/C6H6 

0.115/0.0085 

C2VC6H6 

0.016/0.0012 

Acetylene 

4.9 

Acetylene 

0.6 

Acetylene 

0.069 

Toluene ,  C-,HQ 

*  0.05 

Toluene,  C^H0 

*  0.035 

Toluene,  C-,H„ 

*  0.016 

Toluene ,  C^Hg 

*  0.0175 

C2VC7H8 

*  3.3/0.51 

Acetaldehyde,  CH^CHO 

*  0.105 

All  fuels  were  diluted  in  argon  and  (except  for  the  experiments  on  pure 
toluene)  were  conducted  at  total  pressures  of  approximately  six  to  nine 
atmospheres . 


*  Research  conducted  under  Corporate  -  sponsored  program. 


TABLE  II 


SERIES  OF  EXPERIMENTS  COMPLETED 
DURING  SECOND  YEAR  OF  PROGRAM 


Initial 


Reactant 

Concentration  ! 

Ethene ,  C_H. 

2  4 

3.5 

Ethene , 

0.3 

Ethene , 

0.05 

Styrene ,  C.H 

O  O 

0.4 

Styrene,  CgHg 

0.2 

C2V°2 

3.3/2.05 

C2V°2 

0.22/0.13 

1 , 3 , 5-Hexatriene ,  CgHg 

1.0 

1 , 3 , 5-Hexatriene ,  C^Hg 

0.09 

Acetaldehyde,  CH^CHO 

* 

0.01 

Biacetyl,  (CH3C0)2 

* 

0.22 

Biacetyl,  (CH^CO)2 

* 

0.022 

c2h2/(ch3co)2 

* 

3.5/0.18 

c2h2/(ch3co)2 

* 

0.3/0.15 

Toluene ,  C^Hg 

* 

1.0 

Toluene ,  C-,H0 

* 

0.1 

Toluene,  C^H0 

/  O 

* 

0.01 

Research  conducted  under  Corporate -sponsored  program 
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TABLE  III 

Series  of  Experiments  Completed 
During  Third  Year  of  Program 

Reactant  Initial  Concentration 


trans-3-penten-l-yne,  CcH, 

J  O 

trans-3-penten-l-yne,  C  H 

J  o 

cis-3-penten-l-yne,  C^Hg 

cis-3-penten-l-yne,  C^Hg 

ethylene /oxygen, 

ethylene/oxygen,  02^/02 

acetylene/benzene,  C.H./C  H, 

Li  DO 

Toluene ,  C^Hg 
Toluene,  C^Hg 
Pyridine,  C^H^N 
Pyridine , 

Pyrrole,  C^N 
Pyrrole,  C^H,_N 


1.0% 

0.05% 

1.0% 

0.052% 

3. 5V0. 52% 

2  5/ 0.03^ 
3%/0.3% 

*+1.0? 

#+p  .  OP.V'". 

*0.7% 

*0.lU8% 

*0.5% 

*0.075% 


All  fuels  were  diluted  in  argon  and  were  conducted  at  total  pressures 
of  approximately  six  to  nine  atmospheres. 


*  Research  conducted  under  corporate  -  sponsored  program 

+  Sampling  system  and  gas  chromatograph  modified  to  collect  and  analysis  for 
polyaromatic  species 


TABLE  V 


REACTIONS  AND  RATE  COEFFICENTS  FOR  STEPS 
LEADING  TO  FORMATION  OF  AROMATIC  RING 


Low  Temperature  Sequence 

Log_A 


H  +  C2H2  « 

C2H3 

12.74 

C-H-  +  C-H.  s* 

n-C, Hc 

12.05 

2  3  2  2 

4  5 

n-C4H5  +  C2H2 

f‘C6H7 

12.81 

-  c-C6H7 

11.36 

C-C6H7 

a  C.H,  +  H 

0  O 

13.12 

High  Temperature 

Sequence 

H  +  C2H2  ^ 

c2h  +  h2 

13.6 

C2H  +  C2H2  “ 

n-C4H3 

13.56 

n-C4H3  +  C2H2  * 

K6H5 

11.71 

a 

c-C6H5 

10.22 

c-C6H5  +  H2  a 

C6H6  +  H 

12.39 

c-C6H5  +  C2H2  a 

C8H6  +  H 

12.00 

k  -  A  exp  (-E/RT) 

Units:  A,  cc,  mole,  sec 


E/1000 

2.5 

4.0 

9.0 

0.4 

24.6 


20.5 

3.0 

-0.1 

1.4 

9.5 
4.0 


E,  cal/mole/K 


TABLE  VI 


Reactions  Leading  to 
Formation  of  Polycyclic  Hydrocarbons 


+  H 


+  c2h2  - ►  £ 


+  C2  ^2 


#« 


:C  =  c\ 

^  H 


3.5%  ETHENE  PYROLYSIS,  AUPHATICS 


800  1000  1200  1400  1600  1800  2000 

INITIAL  POST-SHOCK  TEMPERATURE  (K) 


3.5%  ETHENE  PYROLYSIS,  MODEL 


800  1000  1200  MOO  1600  1800  2000 

INITIAL  POST-SHOCK  TEMPERATURE  (K) 


EMM 


ISoft 
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3.5%  ETHENE  PYROLYSIS,  AUPHATICS 


3.5%  ETHENE  PYROLYSIS,  MODEL 
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3.5%  ETHENE  PYROLYSIS,  AROMATICS 


3.5%  ETHENE  PYROLYSIS,  MODEL 


BENZENE  FORMATION 
3.5%  Ethene 


TEMPERATURE  (K) 
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KWU'U  1/  W  V 


jir^iTi,  »yjy*YFVS  SIX  UV  W.H/'.W  HT  VJCWJPJTIl 

ss 


: 


THE  PYROLYSIS  OF  ACETYLENE  INITIATED  BY  ACETONE 

by 

M.  B.  Colket,  111.,  H.  B.  Palmer*  and  D.  J.  Seery 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 

*  Pennsylvania  State  University 
University  Park,  PA  16802 
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A  detailed,  chemical  kinetic,  radical-chain  mechanism  is  used  to  model  the  pyrolysis  of  accetylene 
near  1000K.  The  assumption  of  a  bimolecular  initiation  process  is  found  to  be  inconsistent  with 
thermochemistry.  Since  experimental  evidence  convinces  the  authors  of  the  presence  of  a  chain 
mechanism,  alternative  sources  of  initiation  are  considered.  Acetone,  a  common  impurity  even  in 
‘purified’  acetylene  was  found  to  dominate  radical  initiation  during  the  pyrolysis  of  acetylene  near 
1000K  despite  concentration  levels  only  0.1%  that  of  acetylene.  Modeling  results  compare  favorably 
with  the  experimental  results  of  Munson  and  Anderson  for  acetylene  decay  and  the  formation  of 
products  vinylacetylene,  benzene,  and  ethylene.  Rate  constants  for  the  following  reactions  are  found 
to  be 


C2H3  +  C2H3 
n  -  C4H5 
n  -  C4H5  +  C2H2 
n  —  C4H5  +  C4H4 


C4H6  k4  =  2.0  x  1013  cc/mole/sec 

C4H4  +  H  k32  =  1.6  x  10“  X  e~{41  4  kc«l/«nole/RT)  8e(.-l 

C6H6  +  H  kS3  =  4.5  x  1012  x  e~(1°  kc,,'mo,e/RT)cc/rnole/sec 

C8H8  +  H  k35  =  1.2  x  1013  x  e~(5  k«'/moi«/RT)  cc/mole/sec 


to  optimize  the  fit  to  the  experimental  data.  A  sensitivity  analysis  shows  that  the  computed  results 
were  most  sensitive  to  rate  constants  for  these  and  a  few  other  reactions. 

By  addition  of  high  temperature  reactions  (involving  ethynyl  radicals),  the  mechanism  was  found  to 
be  consistent  with  experimental  results  using  single-pulse  shock  tubes  and  the  very  high  temperature 
study  of  Frank  and  Just.  Calculations  are  quantitatively  consistent  with  their  experimental  finding 
that  acetone  plays  a  negligible  role  at  high  temperature,  low  concentration  conditions. 

IntiodMCtign 

Speculation  on  the  mechanism  of  acetylene  pyrolysis  has  occurred  for  many  years.  Experiments 
at  low  temperature  (600  -  1200K)  during  the  ‘50’s  and  ‘60’s  resulted  in  a  variety  of  suggestions 
including  radical-chain  processes  and  non-chain  dimerizations  involving  diradical  or  electronically 
excited  intermediates.  During  the  ‘70’s,  a  radical-chain  mechanism  was  preferred1  ~s  and  initiation 
was  usually  believed  to  be  due  to 


|  C,H2  +  C2H2 

. 

* 

> 


— *  n  -  C4H3  +  H 


A-  1 


(1) 


»»  « 


where  11-C4H3  is  the  HCCCHCH  radical.  The  low  A-factor  and  activation  energy  used  to  describe 
the  experimental  data  (e.g.,  A  =  2  x  1017  cc/mole-sec,  E  =  45.9  kcal/mole,  Ref.  1)  were  consistent, 
with  the  type  of  reaction  as  well  as  thermochemistry  widely  in  use  at  the  time.  However,  the  heat  of 
formation  of  n-C4Hs  is  now  believed  to  be  close  to  125  kcal/mole*,  much  higher  than  the  previously 
used  value  of  102  kcal/mole  (equated  to  the  early  calculation  for  the  HCCCCH2  radical  by  Duff 
and  Bauer*).  Thus,  the  endothermicity  of  Reaction  1  is  about  68  kcal/mole  and  not  46  kcal/mole 
as  previously  believed.  Attempts  (at  UTRC)  to  model  acetylene  pyrolysis  with  Ej  e=  68  kcal/mole 
leads  to  an  unacceptably  high  A-factor  for  this  reaction.  This  result  obviously  raises  serious  doubt 
regarding  the  chain  mechanism. 

Furthermore,  recent  results6,7  suggest  that  acetylene  pyrolyzes,  by  first  isomerizing  to  vinylidene 
(HjCC:)  which  then  inserts  into  a  C-H  bond  of  acetylene  to  form  vinylacetylene  directly.  Duran, 
et  al.  interpreted  experimental  data  to  reach  this  conclusion  while  Kiefer  argues  this  sequence  is 
consistent  with  thermodynamics. 

In  our  own  analysis,  we  cannot  disprove  a  non-chain  dimerization;  however,  upon  review  of  our  own 
experimental  data  as  well  as  that  obtained  by  others  we  reach  the  unavoidable  conclusion  that  a 
radical-chain  process  occurs.  In  attempts  to  model  experimental  data,  we  argue  that  thermal  de¬ 
composition  of  acetone,  a  common  impurity  in  acetylene,  is  quite  probably  responsible  for  initiating 
the  radical-chain  in  laboratory  studies  of  the  pyrolysis  of  acetylene. 

In  this  paper,  evidence  supporting  a  radical-chain  will  be  presented,  issues  related  to  acetone- 
initiation  are  discussed,  and  modeling  results  are  presented  and  compared  to  existing  experimental 
data. 

Evidence  for  Chain  Mechanism 

Palmer  and  Cullis’  review  article8  cites  results  such  as  inhibition  by  radical  scavengers,  wall  effects, 
and  existence  of  an  induction  period  as  evidence  that  the  thermal  decomposition  of  acetylene  is  a 
radical  process.  Unfortunately,  these  observations  cannot  distinguish  between  diradical  (non-chain) 
and  monoradical  (chain)  mechanisms. 

Evidence  for  a  chain  arises  from  the  work  of  Callear  and  Smith®  who  investigated  the  reaction 
between  acetylene  and  atomic  hydrogen  in  a  quartz  vessel  at  293K.  Atomic  hydrogen  was  generated 
by  decomposing  H2  using  a  low-pressure  mercury  discharge  lamp.  Major  products  included  ethene, 
1,3-butadiene,  benzene,  and  trans-l,3,5-hexatriene.  To  explain  these  products,  Callear  and  Smith 
suggested  a  reaction  mechanism  similar  to  that  used  in  the  present  study  except  that  in  the  low 
temperature,  hydrogen-rich  investigation,  CsH*,  n-C4Hg,  and  1-CeHr  (1,3,5-hexatrien-l-yl)  radicals 
formed  principally  CjH4,  1,3-C4Hb,  and  1,3,5-CeHg  rather  than  losing  H-atoms  to  form  CjHa, 
C4H4  and  benzene.  Some  benzene  was  still  formed  at  low  temperatures  presumably  because  of  the 
high  stability  of  the  aromatic.  If  a  chain  process  exists  at  room  temperature,  it  should  also  occur 
at  elevated  temperatures. 

At  higher  temperatures,  there  also  is  a  multiplicity  of  products  from  the  thermal  decomposition  of 
acetylene.  Near  1000K,  using  20%  C2H2  in  helium,  Munson  and  Anderson10  found  vinylacetylene 
and  also  benzene,  ethene,  methane  and  hydrogen.  Working  at  even  higher  temperatures,  Ogura5 
detected  the  same  products  during  pyrolysis  in  a  single  pulse  shock  tube  and  also  diacetylene  and 
traces  of  1,3-butadiene,  and  C$  and  Cs-species.  Under  similar  conditions  but  with  a  heated  sampling 
system,  Colket11  detected  the  same  species  and  benzene,  phenylacetylene,  and  traces  of  naphthalene 
as  well.  Although  the  vinylidene  mechanism  apparently  explains  the  overall  decomposition  rates  of 


the  reactant,  it  does  not  describe  the  variety  of  products  that  have  been  observed  experimentally. 
Thus,  contributions  from  a  chain  sequence  are  expected. 

Further  indications  of  a  chain  mechanism  are  provided  by  additional  single-pulse  shock  tube  data 
in  Fig.  1.  This  figure  compares  vinylacetylene  and  benzene  production  from  three  different  gases: 
(a)  2040  ppm  biacetyl  in  argon,  (b)  3.7%  acetylene  in  argon  and  (c)  a  mixture  of  3.55%  acetylene 
and  1500  ppm  biacetyl  in  argon.  The  production  of  C4H4  and  CeHe  is  substantially  larger  from  the 
mixture  than  would  be  expected  from  purely  additive  considerations.  Near  1I00K,  the  mixture  pro¬ 
duces  an  order  of  magnitude  greater  C4H4  and  CeHe  than  are  produced  from  acetylene  or  biacetyl. 
Since  there  is  no  reason  to  expect  that  biacetyl  enhances  the  rate  of  the  vinylidene  mechanism,  this 
enhancement  must  be  a  result  of  a  chain  mechanism  initiated  by  biacetyl  decomposition  first  into 
acetyl,  then  to  methyl  radicals.  Mechanisms  by  which  methyl  radicals  are  converted  to  H-atoms 
(required  for  the  acetylene  chain)  are  discussed  elsewhere  in  this  manuscript.  Methyl  radical  addi¬ 
tion  to  acetylene  and  to  subsequent  adducts  forms  a  variety  of  products,  e.g.,  C3H4,  C4H6,  C4H8, 
etc.,  but  the  concentrations  of  these  species  are  small  relative  to  those  of  C4H4  and  CeHe  There 
appears  to  be  no  facile  method  to  produce  vinylacetylene  and  benzene  at  the  observed  rates  other 
than  from  the  chain  mechanism  related  to  acetylene  pyrolysis. 

Finally,  Ogura3  shock-heated  equimolar  mixtures  of  C2H2/C2D2  diluted  in  argon.  He  found  over 
40%  of  the  vinylacetylene  in  the  form  of  the  C4H3D  and  C4HD3  isotopes  at  temperatures  as  low  as 
1070K.  It  is  highly  unlikely  that  such  significant  scrambling  can  be  accounted  for  by  a  vinylidene 
mechanism.  Again  a  radical-chain  sequence  is  favored. 

In  summary,  (1)  a  radical-chain  mechanism  has  been  observed  at  room  temperature;  (2)  acety¬ 
lene  pyrolysis  produces  a  wide  variety  of  products;  (3)  a  large  enhancement  in  the  production  of 
vinylacetylene  and  benzene  is  observed  when  acetylene  is  a  co-pyrolyzed  with  biacetyl,  and  (4) 
significant  scrambling  was  observed  during  pyrolysis  of  C2H2/C2D2  mixtures.  Since  these  phenom¬ 
ena  are  either  indicative  of  a  chain  mechanism  or  can  only  be  explained  by  a  chain  mechanism, 
this  analysis  strongly  supports  the  existence  of  a  chain  during  acetylene  pyrolysis  at  temperatures 
of  900-1400K.  It  does  not  exclude  possible  contributions  of  a  vinylidene  route  or  other  non-chain 
mechanism,  nor  does  it  provide  information  on  the  relative  contributions  of  the  two  mechanisms  if 
they  coexist.  Nevertheless,  since  recent  modeling  efforts  have  focused  on  a  non-chain  mechanism, 
this  work  focuses  solely  on  chain  processes  to  examine  whether  or  not  acetylene  decomposition  and 
product  formation  can  be  described  adequately. 

Preliminary  Modeling 

For  low  temperature  pyrolysis  (900-1400K),  Tanzawa  and  Gardiner*  invoked  bimolecular  initiation 
followed  by  addition  of  H-atoms  to  C3H2  to  form  vinyl  radicals.  The  proposed  chain  mechanism 
accounted  quite  well  for  the  experimental  results  fo?  acetylene  decomposition.  However,  it  did  not 
describe  the  known  production  of  benzene  and  higher  molecular  weight  species,  which  together 
accounted  for  more  than  50  percent  of  the  decomposed  acetylene  (near  1000K).  Frenklach  and 
coworkers12  extended  this  mechanism  to  describe  formation  of  heavier  species  including  PAH  and 
soot  formation  during  pyrolysis  of  acetylene  at  higher  temperatures.  Colket11,  subsequently,  used 
a  similar  mechanism  for  comparison  to  experimental  profiles  of  species  with  molecular  weights  up 
to  102. 

An  unsettling  feature  of  these  models  is  the  rate  expression  for  the  initiation  step  (Reaction  1).  Using 
a  recent  value4  for  the  heat  of  formation  of  n-C4Hs  of  124  kcal/mo  Reaction  1  is  68  kcal/mole 
endothermic.  However,  the  activation  energy  used  previously  for  this  reaction  is  about  46  kcal/mole, 


significantly  lower  than  its  endothermicity.  Hence,  we  attempted  to  model  experimental  data  using 
Ei  =  68  kcal/mole.  This  value  not  only  is  required  to  satisfy  thermochemistry  but  also  is  attractive 
when  one  performs  a  steady-state  analysis  on  the  simplified  sequence  presented  in  Table  I.  Such  an 
analysis  leads  to 

-  -Cd\H2)  =  2k3A  (^) 4  [C,H,]’  +  2ki[CjHj]1  Eq.  A 

in  agreement  with  the  often  quoted  experimental  observation  of  second  order  decomposition  kinetics. 
(This  sequence  differs  from  that  proposed  by  Tanzawa  and  Gardiner1  in  that  vinyl  recombination 
was  found  to  dominate  radical  termination.  In  the  earlier  modeling1  the  high  stability  assumed  for 
C4H3  led  to  its  important  role  in  termination  and  C4H3  +  H  dominated.  Termination  by  vinyl 
recombination  has  also  been  suggested  by  Ogura8.)  Assuming  a  long  chain,  the  second  term  in  Eq. 
A  can  be  ignored  and  the  overall  activation  energy  is  Eov  =  E3a  +  (Ei-E4)/2.  Assuming  E3A  =  4, 
E<  =  0  kcal/mole  and  taking  Ei  =  68,  then  Eov  =  38  kcal/mole  which  is  in  excellent  agreement  with 
measured  values.  The  ‘old’  value  for  Ej  of  46  kcal/mole  gives  a  steady-state  E0v  of  27  kcal/mole, 
substantially  lower  than  experimental  results. 

Despite  the  attractiveness  of  the  ‘updated’  Ei,  our  attempts  to  model  experimental  results  were 
not  fully  satisfying  Quite  reasonable  agreement  could  be  obtained  with  both  the  decay  profiles  of 
acetylene  and  product  formation  as  measured  by  Munson  and  Anderson10.  The  mechanism  used 
was  similar  to  that  used  by  Colket11  except  that 

kjigh  =6x  10He(-68000e‘,/mole/RT)cc/mole/8ec 

and 

k<  =  2  x  1012cc/mole/sec. 

The  first  expression  is  denoted  k^'*h  since  it  is  higher  than  the  maximum  value  one  would  expect. 
The  A-factor  is  obviously  much  too  high  for  a  bimolecular  initiation  of  this  type.  The  value  of 
k«,  for  vinyl-vinyl  recombination,  is  more  than  five  times  lower  than  that  typically  used  for  this 
reaction.  Consequently,  these  results  must  be  viewed  with  caution. 

Indications  of  Acetone  Initiation 

With  concern  about  the  above  results  yet  with  the  evidence  for  the  existence  of  a  chain,  possible 
initiation  steps  were  re-examined.  Of  the  likely  possibilities 

C2H2  — *  C2H  +  H  -  128kcal/moie 
C2H2  +  C2H2  —  C:H  +  C2H3  -  87kcal/mole 
C2H2  +  C2H2  — *  i  -  C4H3  +  H  -  60kcal/mole 

only  the  third  has  a  sufficiently  low  endothermicity,  yet  it  also  would  be  expected  to  have  a  low 
A-factor  in  order  to  form  the  H2CCCCH  radical.  This  reaction  remains  a  possibility  and  may 
contribute.  Further  consideration  of  the  initiation  problem  led  us  to  consider  impurities  as  potential 
initiators.  Acetone,  stored  in  all  acetylene  bottles  as  an  inhibitor,  is  a  well-known  impurity  of 
acetylene  and  coelutes  from  tanks  in  concentrations  of  one  to  20%  that  of  acetylene18.  When 
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precautions  are  taken  to  reduce  the  concentration  of  impurities,  acetone  concentrations  are  typically 
decreased  to  approximately  0.1%  of  the  acetylene  concentration.  Assuming  that  a  radical  initiation 
rate  of  k j'gh  (CjHjj2  is  needed  to  produce  sufficient  radicals  for  acetylene  pyrolysis,  then  this  rate 
can  be  compared  to  the  known14  initiation  by  acetone. 


CH3COCH3 

CHsCO 


CH3  +  CH3CO  (5) 
CHS  +  CO  (6) 


Assuming  [acetone]/[C2H2]  =  0.001,  one  obtains 


rfe°n^  =  0.0004  x  e(-4000  «'/”°'*/"T)[C2H2rl 

k  1 1 C  2  H  2  ] 


with  Szwarc  and  Taylor's  expression14  for  ks  (=2.4  x  1015  exp(-72000  cal/mole/RT)  sec-1)  At 
20%  C2H2  and  1000K  (Munson  and  Anderson’s  conditions),  this  ratio  is  22  while  at  [C2H2]  =  2.35 
x  10~5  mole/cc  and  1200K  (Ogura),  the  ratio  is  3.2.  Since  this  ratio  is  greater  than  one  at  both 
experimental  conditions,  radical  initiation  by  thermal  decomposition  of  acetone  dominates  over  that 
of  acetylene 

Considering  that  k  J‘8h  is  substantially  higher  than  can  be  expected  from  thermodynamic  arguments 
and  that  the  radical  initiation  rate  from  acetone  is  higher  than  the  required  rate  of  initiation,  it 
seems  quite  reasonable  to  expect  that  acetone  plays  a  role  in  the  initiation  process.  For  reference,  it 
is  instructive  to  examine  the  ratio  k^acetone] /k  1  [CsH*)2  when  kj  is  assigned  the  value  determined 
by  Tanzawa  and  Gardiner1.  This  ratio  is  then  0.13  and  0.12  for  the  previously  discussed  conditions. 
These  values,  in  all  probability,  led  to  early  conclusions  that  acetone  initiation  contributed  negligibly 
to  acetylene  pyrolysis.  This  major  difference  in  conclusions  suggests  a  major  difference  in  the  two 
values  of  k!  and  reasons  for  this  difference  need  to  be  examined. 

The  absolute  value  of  kj  determined  by  Tanzawa  and  Gardiner  is  nearly  two  orders  of  magnitude 
higher  than  that  found  in  the  present  analysis.  The  principal  reason  for  their  higher  value  is  be¬ 
lieved  to  be  a  result  of  the  low  AH /  used  for  the  n-C4Hj  radical  (~  20  kcal/mole  lower  than  that 
used  in  the  present  analysis).  The  resultant  high  stability  of  this  radical  as  assumed  in  their  model 
led  to  high  rates  of  radical  termination  via  n-C4Hs  +  H  +  M.  This  step  played  only  a  minor  role  in 
the  present  study;  rather  vinyl-vinyl  recombination  dominated  termination.  Thus  radical  termina¬ 
tion  was  overestimated  substantially  in  the  earlier  studies  and  the  initiation  rate  was  increased  to 
compensate. 

The  important  conclusion  from  this  cursory  analysis  is  that  acetone  decomposition  is  fast  enough 
to  influence  the  radical  concentrations  during  acetylene  pyrolysis. 

Issues  Raised  by  Acetone  Initiation 

Several  critical  issues  must  be  raised  if  initiation  by  acetone  contributes.  First  of  all,  a  rapid 
mechanism  for  conversion  of  methyl  radicals  to  H-atoms  must  be  present;  secondly,  the  ‘well-known’ 
second  order  behavior  of  acetylene  must  be  addressed,  specifically  to  see  if  this  reaction  order  can 
be  predicted  and  if  not,  why  not;  and  thirdly,  if  acetone  initiation  dominates,  then  why  is  there 
such  good  experimental  agreement  (see  Ref.  8)  on  acetylene  decomposition  from  study  to  study 
when  initial  concentrations  of  acetone  can  be  expected  to  vary. 
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It  is  readily  shown  that  methyl  radicals  can  be  rapidly  converted  to  H-atoms.  A  mechanism  for 
conversion  is  shown  in  Table  II.  Rate  constants  were  obtained  from  the  literature  where  possible. 
The  principal  fate  of  methyl  near  1000K  will  be  its  addition  to  acetylene  (R7).  According  to  a  recent 
QRRK  analysis15  the  resultant  adduct  should  stabilize  under  these  conditions.  Using  rate  constants 
(in  Table  II  and  IV)  and  thermodynamic  parameters  nearly  identical  to  those  given  by  Dean  and 
Westmoreland15  (see  Table  VI)  the  rate  of  CH3CHCH  isomerization  to  ally  I  radical  relative  to  its 
decomposition  to  methyl  acetylene  and  H  is  approximately  five.  Due  to  the  stability  of  allyl,  its 
decomposition  is  slow  and  principally  adds  to  acetylene  to  form  C5H7  (R9).  This  linear  radical  can 
then  either  decompose  to  reactants  (R(-8))  or  form  cyclopentadiene  plus  atomic  hydrogen  (RIO). 
Thus,  methyl  radicals  produced  from  acetone  decomposition  will  be  rapidly  converted  to  H-atoms 
Since  the  k[acetone]/ki[CjHj]2  ratios  as  previously  calculated  are  significantly  greater  than  one, 
production  of  H-atoms  via  acetone  decomposition  appears  to  be  sufficiently  fast  to  initiate  the 
thermal  decomposition  of  acetylene. 

At  higher  temperatures,  methyl  radical  addition  to  acetylene  must  compete  with  methyl  recombina¬ 
tion  (R17),  since  the  rapid  decomposition  of  acetone  quickly  produces  high  concentrations  of  methyl. 
Also,  methyl  can  extract  hydrogen  from  acetylene  (R16),  since  at  high  temperatures  thermal  energy 
is  available  to  overcome  the  substantial  endothermicity  (~  25  kcal/mole)  of  this  reaction.  In  either 
high-temperature  case,  a  variety  of  alternative  sequences  will  again  provide  H-atoms 

The  order,  n,  for  the  overall  reaction  as  defined  by 

J  =  kov(C2H2]" 

is  often  cited  as  two.  A  steady-state  analysis  using  reactions  in  Table  I  and  Reactions  5-10  in  Table 
II  leads  to 


-~^t~  =  [CjHj] [acetone]* 


or  assuming  [acetone]/[C2H2)  =  0.001,  then 


-d|CdtH2]  =0.0894  ksA^HCjH,]1 


Eq.  B 


This  steady-state  result  that  n=3/2  appears  to  be  a  major  drawback  of  the  present  proposal  re¬ 
garding  initiation  by  acetone.  A  review  of  the  literature  was  performed  in  an  attempt  to  find 
references  in  which  the  reaction  order  was  determined  or  which  provide  data  from  which  the  order 
could  be  computed.  A  summary  of  reaction  orders  is  presented  in  Table  III.  Although  there  is 
substantial  variation  in  the  reaction  order,  most  values  fall  in  the  range  of  1.5  to  2.0.  The  low 
temperature  data  from  Silc.ocks16  as  well  as  the  flow  reactor  data  from  Palmer  and  Dormiah17  were 
originally  interpreted  assuming  simultaneous  homogeneous  and  heterogeneous  processes.  Silcocks 
found  the  reaction  orders  to  be  two  and  one,  respectively,  while  Palmer  and  Dormish  found  that 
this  assumption  was  consistent  with  their  data.  An  attempt  has  been  made  to  reinterpret  the  lat¬ 
ter  experiments,  with  a  resultant  overall  order  ranging  from  1.2  to  1.7.  Munson  and  Anderson’s 
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low  temperature  experiments10  exhibit  significant  scatter  in  plots  of  log(rate)  vs.  log[CjHj]  but 
the  order  appears  to  be  1.7  ±  0.2.  Their  higher  temperature  experiments  clearly  give  n  =  1.5  ± 
0.1.  The  results  from  Colket11  of  n  =  1.6  ±  0.2  are  from  comparison  of  single-pulse  shock  tube 
experiments  at  4.9%  (unpublished),  3.5%  (Ref  11)  and  0.2%  acetylene  (unpublished)  The  order 
obtained  by  Ogura3  is  substantially  higher  than  from  other  experiments.  Results  from  Cullis  and 
Franklin18  at  983K  give  beautiful  agreement  with  the  assumption  of  a  second  order  reaction.  If 
indeed  a  chain  mechanism  controls  the  reaction,  however,  then  Cullis  and  Franklin's  data  needs 
reinterpretation.  They  pyrolyzed  pure  acetylene  at  total  pressures  ranging  from  96  to  402  torr. 
Under  these  conditions,  some  critical  reactions,  particularly  H  +  CjHj  «-*  CjH3  can  be  expected  to 
exhibit  pressure  dependence  (see  Warnatz19)  and  will  lead  to  a  net  decrease  in  overall  reaction  rate 
with  decreasing  pressure.  If  pressure  dependence  is  included,  a  preliminary  analysis  indicates  that 
Cullis  and  Franklin’s  data  are  no  longer  inconsistent  with  the  assumption  of  an  overall  reaction 
order  of  1.5.  Results  from  other  authors  listed  in  Table  III  do  not  support  the  ‘well-known’  reaction 
order  of  two. 

Thus,  with  the  exception  of  the  very  low  temperature  data  of  Silcocks  and  the  value  from  Ogura,  the 
assumption  of  an  overall  reaction  order  of  1.5  as  suggested  by  acetone  initiation  is  not  unreasonable. 
In  fact,  the  experimental  data  seem  to  support  a  value  of  1.5  as  well  as  or  perhaps  better  than  the 
order  2. 

If  acetone  is  the  initiator,  then  why  should  experiments  performed  in  many  different  laboratories, 
presumably  with  different  initial  acetone  concentrations,  all  result  in  similar  overall  rate  constants? 
This  is  a  question  which  we  could  only  partially  address.  A  variety  of  purification  techniques  have 
been  used.  Often  the  level  of  impurities  was  not  determined  or  not  determined  with  sufficient 
accuracy.  Based  on  our  own  experience,  the  work  of  Hamins,  et  a!,s  and  reported  purity  levels,  we 
estimate  that  typical  purified  samples  of  acetylene  contain  about  0.1%  acetone,  with  a  total  range 
of  0.05%  to  0.2%.  This  variation  of  a  factor  of  four  leads  to  a  factor  of  two  variation  in  overall 
decomposition  rates  (due  to  the  half-order  dependence  on  acetone  concentration).  This  variation 
in  overall  rates  is  within  the  measured  experimental  uncertainty. 

Development  of  Kinetic  Model 

The  kinetic  model  used  in  this  analysis  is  based  on  (1)  the  initiation  sequence  as  indicated  in 
Table  II  largely  using  published  rate  constants  and  (2)  an  acetylene  pyrolysis  model  as  refined  by 
Colket11  to  describe  results  obtained  in  a  single-pulse  shock  tube.  Model  results  (using  a  program 
based  on  CHEMK1N20,  LSODE21  and  isothermal,  isobaric  assumptions)  were  compared  with  the 
flow  reactor  data  of  Munson  and  Anderson10.  The  mechanism  and  rate  constants  used  to  model 
the  flow  reactor  data  are  shown  in  Tables  II  and  IV.  Bimolecular  initiation  by  acetylene  (Rl)  was 
included  in  the  model  using  a  rate  constant  which  can  be  expected  to  be  an  upper  limit.  In  the 
modeling  of  the  Munson  and  Anderson  data,  this  reaction  always  had  a  negligible  role.  A  few  rate 
expressions  have  been  modified  in  order  to  provide  better  fits  to  the  flow  reactor  data.  Selection  of 
rate  constants  for  modification  was  accomplished  using  a  combination  of  a  reaction  pathway  analysis 
(with  computer  programs  developed  at  UTRC)  and  a  sensitivity  analysis  using  CHEMSEN22.  For 
each  rate  expression  modified,  attempts  were  made  to  maintain  consistency  with  the  SPST  modeling 
results  by  adjusting  both  the  A-factor  and  the  activation  energy.  In  most  cases  this  was  possible 
since  the  flow  reactor  data  were  most  dependent  on  kinetics  over  the  range  900-1000K,  whereas  for 
the  SPST  results,  the  range  was  1200-1500K.  In  the  case  of  the  competitive  sinks  for  the  n-C<H& 
radical, 
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this  was  not  possible,  but  an  attempt  was  made  instead  to  maintain  the  ratio  of  rate  constants 
as  determined  in  modeling  the  SPST  data.  Reaction  33  has  been  simplified  from  the  previous 
3-step  reversible  sequence  involving  intermediate  formation  of  £-CeHr  and  c-CeHr,  since  a  QRRK 
analysis  by  Westmoreland23  indicates  that  at  the  conditions  of  the  present  study,  the  reaction 
rapidly  proceeds  as 

c2h2  +  n  -  c4h5  —  t  -  c6h;  -  c  -  c6h;  c6h6  +  h 


Modeling  results  near  1000K  using  the  irreversible  one-step  or  the  three-step  reversible  sequence 
were  virtually  indistinguishable 

Other  differences  between  the  present  and  previous  model  include  rate  expressions  for  H-atom 
extraction  by  phenyl  or  vinyl  from  vinylacetyiene.  These  were  equated  to  the  expression  for  vinyl 
plus  1,3-butadiene  suggested  by  Kiefer,  et  al  24 .  Model  results  showed  little  sensitivity  to  rate 
constants  for  these  reactions.  Vinyl-vinyl  recombination  (R4)  was  taken  to  be  2  x  1013  cc/mole/sec 
based  on  a  fit  to  the  Munson  and  Anderson  data.  Rate  constants  for  other  alkyl  radical  terminations 
were  assumed  to  be  equal.  The  selected  value  was  5  x  1012  cc/mole-sec.  After  Frenklach,  et  al  ’*, 
the  two  isomers,  n-C4H3  and  i-C4H3  were  included;  however,  there  is  no  ready  sink  for  i-C4H3  at 
the  low  temperatures  of  the  Munson  and  Anderson  study.  To  resolve  this  problem,  two  reactions 
were  included. 

n  -  C4H3  ♦-*  i  -  C4H3  (61) 
i  -  C4H3  +  C4H4  -»  CsHe  +  H  (66) 

The  first  of  these  was  found  to  be  negligible  near  1000K.  The  second  was  the  dominant  sink  for 
i-C4H3  while  affecting  phenylacetylene  slightly.  The  first  reaction  was  assigned  a  rate  constant 
similar  to  a  value25  for  the  i-C3H?  «-*  n-C3H7  isomerization,  although  an  activation  energy  close 
to  the  strain  energy  in  cyclopropene  was  used.  The  second  reaction  has  been  proposed™  as  an 
alternative  mechanism  for  forming  phenyl  acetylene.  It  is  a  multistep  process  requiring  a  1,1  11- 
atom  shift,  cyclization,  and  loss  of  an  H-atom  and  was  assigned  a  rate  expression  identical  to  that 
of  Reaction  4  which  requires  a  similar  transformation.  Neglecting  the  back  reaction  is  justified 
since  the  reverse  process  is  slow  at  these  temperatures.  A  QRRK  analysis*7  has  been  performed 
for  Reaction  66  and  indicates  that  near  1000K  approximately  50%  of  the  C4H4-i-C4H3  adducts 
will  go  directly  to  products  prior  to  stabilization  to  intermediates.  Therefore  elimination  of  the 
intermediate  steps  is  reasonable. 

A  third  possible  sink  for  i-C4H3  not  included  in  the  modeling  is  the  addition  of  i-C4H3  to  acetylene, 
cyclization  to  form  a  five-membered  ring  and  then  isomerization  to  phenyl. 

Third  body  effects  have  been  assumed  to  be  small  for  this  analysis.  The  greatest  deviation  from 
this  assumption  will  occur  for  the  H-atom  addition  to  acetylene.  Fall-off  curves  given  by  Warnatz*9 
indicate  that  the  k/koo  for  this  reaction  is  approximately  0.5  at  one  atmosphere  of  argon  and  1000K. 


In  the  present  study,  20%  of  the  gas  is  acetylene,  a  more  efficient  third  body,  suggesting  the  k/koo 
is  larger  than  0.5;  therefore,  neglecting  pressure  effects  is  believed  to  result  in  only  a  small  error  A 
reduction  of  k2  by  a  factor  of  two  results  in  a  reduction  of  the  overall  decomposition  rate  of  about 
5-10%. 

Some  modeling  was  also  performed  for  comparison  to  shock  tube  results.  For  these  calculations, 
the  kinetic  mechanism  was  expanded  to  include  ethynyl  and  related  radicals  (see  Table  V). 

Thermodynamic  values  are  reported  in  Table  VI  and  are  the  same  as  those  used  previously11,  except 
for  those  of  the  C4H3  isomers,  phenyl  radical  and  the  Cg-species  which  were  obtained  from  Stein4. 
The  latter  are  identical  to  those  used  in  the  study  by  Frenklach,  et  al11. 

Kinetic  Model  vs.  Theory 

In  Figs.  2  to  5,  calculations  of  acetylene,  benzene,  vinylacetylene,  and  ethene  are  compared  to  ex¬ 
perimental  profiles  as  measured  by  Munson  and  Anderson10  in  a  flow  reactor.  For  these  calculations 
it  was  assumed  that  the  acetone/acetylene  ratio  was  0.001.  Overall  the  agreement  is  quite  satisfac¬ 
tory.  Benzene  is  overpredicted  at  long  times  for  the  highest  temperature  run,  presumably  because 
of  the  conversion  of  benzene  to  higher  molecular  weight  species.  Vinylacetylene  is  underpredicted 
at  the  low  temperatures.  Hydrogen,  not  shown,  is  underpredicted  by  about  a  factor  of  five  to  ten 
Although  substantial  methane  formation  was  observed,  a  model  for  methane  formation  was  not 
included.  Methane  produced  from  methyl  directly  or  via  odd-carbon  species  derived  from  acetone 
could  only  account  for  a  very  small  portion  of  the  methane  observed.  Heterogeneous  processes 
seem  the  most  likely  source  of  methane  although  gas-phase  reactions  are  possible.  For  example, 
consider  addition  of  the  i-C«H5  radical  (H2CCHCCH2)  to  acetylene  followed  by  cyclization  to  a 
five-membered  ring  and  isomerization  to  methylcyclopentadienyl.  Subsequent  reactions  involving 
methylcyclopentadienvl  may  then  form  methyl  radicals  and/or  methane.  The  i-C^Hs  radical  is  gen¬ 
erally  ignored  in  kinetic  modeling  but  may  play  a  significant  role  since  it  should  be  formed  readily 
via  H-addition  to  vinylacetylene. 

The  formation  of  polyaromatic  hydrocarbons  is  not  included  in  the  present  model.  Modeling  results 
for  acetylene,  benzene,  and  vinylacetylene18  fit  the  experimental  data  almost  perfectly  for  the 
Munson  and  Anderson  data  of  1073K  when  formation  of  poiyaromatic  hydrocarbons  are  included. 
There  is  relatively  little  change  to  the  model  predictions  for  the  lower  temperature  runs. 

Model  calculations  have  also  been  compared  to  data  obtained  in  a  single-pulse  shock  tube  by 
Colket11,M  at  3.5%  and  4.9%  acetylene  in  five  and  eight  atmospheres  of  argon,  respectively.  Ex¬ 
perimental  procedures  for  the  unpublished  4.9%  data  were  similar  to  those  previously  described11. 
The  acetone  concentrations  were  equated  to  the  measured  values,  that  is,  0.1%  and  0.2%  the  con¬ 
centration  of  acetylene  for  the  3.5%  and  4.9%  data,,  respectively.  Calculations  were  performed 
using  a  shock  tube  code  modified  to  account  for  quenching  processes.  Again  satisfactory  agreement 
(comparable  to  that  achieved  in  Ref.  11)  is  obtained  for  most  species,  including  the  minor  products, 
styrene  and  ethylene.  The  main  descrepancy  is  that  the  model  underpredicts  product  formation 
at  low  temperatures.  Of  particular  interest  is  Fig.  6  in  which  only  odd-carbon  products  derivable 
from  the  impurity  acetone  are  shown  for  the  4.9%  series  of  runs.  Also  plotted  is  the  sum  of  the 
concentrations  of  acetone  plus  half  of  all  species  containing  an  odd  number  of  carbons.  Acetone 
produces  two  methyl  radicals,  each  of  which  could  produce  an  odd-carbon  species.  The  fact  that 
this  sum  is  essentially  constant  over  the  range  of  temperature  of  experiments  lends  strong  support 
to  a  major  thesis  of  this  work;  i.e.,  acetone  produces  methyl  radicals  which  convert  to  H-atoms  via 
addition  to  acetylenes,  whereas  ethane  formation  via  methyl  recombination  is  small  compared  to 
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methyl  addition  to  acetylene.  Model  predictions  are  not  shown  in  Fig  6  since  prediction  of  the 
odd-carbon  products  is  generally  poor.  This  limitation  is  probably  a  direct  result  of  the  fact  that 
there  is  very  little  experimental  data  on  these  species,  and  many  of  the  rate  constants  in  Table  II 
were  estimated. 

Ogura  also  obtained  single-pulse  shock  tube  data  at  5%  acetylene  in  about  two  atmospheres  of  argon 
Model  predictions  are  compared  to  the  Ogura3  data  in  Fig.  7.  for  vinylacetylene,  diacetylene  and 
hydrogen.  Again  comparison  between  experiments  and  theory  seem  quite  satisfactory.  For  the 
model  calculations,  dwell  times  of  one  millisecond  were  assumed  and  quenching  rates  comparable 
to  those  measured  by  Colket11  were  used.  The  acetone/acetylene  ratio  was  assumed  to  be  0  001 
Attempts  to  fit  Ogura’s  experimental  profiles  of  methane,  methylacetylene,  allene,  and  ethylene  for 
experiments  with  10%  acetylene  in  argon  were  also  made.  The  agreement  was  good  for  methane  and 
satisfactory  for  ethylene  except  at  low  temperatures  when  the  model  significantly  underpredicted 
the  experimental  data  The  model  overpredicts  the  Cj-hydrocarbons  by  about  a  factor  of  two 
to  three.  This  difference  does  not  seem  unreasonable  considering  the  general  lack  of  knowledge 
regarding  pyrolysis  of  odd-carbon,  acetylenic  species. 

Vinylidene  vs.  Chain 

Perhaps  the  most  significant  pieces  of  experimental  evidence  against  a  radical  chain-mechanism  are 
the  results  by  Hou  and  Anderson39  and  more  recently  by  Duran,  et  al6.  The  earlier  study  involved  a 
molecular  beam,  mass  spectrometric  analysis  of  acetylene  pyrolysis  products  from  a  flow  reactor  near 
1000K.  The  usefulness  of  their  facility  was  ‘proved’  by  its  ability  to  detect  methyl  radicals  during 
pyrolysis  of  tetramethyl  lead  and  di-t-butylperoxide.  They  detected  no  radicals  during  acetylene 
pyrolysis.  Calculations  for  their  experimental  conditions  indicate  that  during  acetylene  pyrolysis 
their  total  radical  concentrations  always  were  substantially  less  than  one  part  per  million.  It  appears 
unlikely  that  their  facility  had  the  sensitivities  to  observe  species  at  such  low  concentrations,  whereas 
methyl  radicals  readily  produced  from  the  test  compounds  would  be  expected  to  achieve  substantial 
concentrations.  We  think  that  these  early  experiments  do  not  disprove  a  radical-chain  mechanism. 

Recent  pyrolysis  experiments  by  Duran,  et  al6.  were  also  performed  in  a  flow  reactor  with  mass 
spectrometric  results.  Additions  of  benzene  and  toluene  showed  respectively,  no  inhibition  and  an 
acceleration  The  acceleration,  unexplained  by  the  vinylidene  mechanism,  is  easily  explained  via 
chain  processes  since  initiation  via  toluene  (CtHs  — *  C7H7  +  H)  using  any  of  the  recently  proposed 
rate  expressions  (see  Ref.  30-31  and  citations  contained  therein)  can  be  shown  to  be  substantially 
faster  than  initiation  by  acetone.  In  addition,  Duran  observed  a  variety  of  vinyl  benzenes  when 
benzene  and  acetylene  were  copyrolyzed.  They  interpreted  these  results  as  evidence  for  insertion 
of  the  carbene  (vinylidene)  into  the  benzene  ring.  Alternatively,  these  vinylbenzenes  may  also  be 
the  result  of  the  displacement  (C3H3  4-  CsHs  — *  H  +  CeHsCjHs),  which  is  similar  to  the  known” 
methyl-substitution  in  the  benzene-toluene  system.  Copyrolysis  of  toluene  and  acetylene  produced 
principally  a  species  at  mass  118,  interpreted  to  be  methyl  styrene.  Instead  this  species  may  be 
phenyl  propene  formed  by  acetylene  addition  to  the  stable  benzyl  radical.  Consequently,  we  question 
their  conclusion  of  proof  for  the  vinylidene  mechanism. 


Overall  Reaction  Order  and  Activation  Energy 


also  determined  from  plots  of  predicted  log  (d[C2H2]/dt)  vs.  [C2H2]  as  a  funtion  of  time  over  the 
temperature  range  873  to  1073  K  These  calculations  resulted  in  a  reaction  order  of  1.43,  1  66,  and 
1.67  at  973,  1023,  and  1073K  respectively.  These  values  are  consistent  with  the  overall  order  of  1.5 
obtained  from  the  steady-state  analysis  Assuming  this  latter  value,  an  overall  rate  expression  was 
determined  using  the  model  prediction  for  Munson  and  Anderson’s  conditions.  The  rate  expression 
determined  was: 


d[C2H2l 


,  )  -  101,e-(40600  ca,/moU/RT){C2H2]1  8  (cc/molej^sec"1 

dt  /overall 


at  873  -  1073  K  This  expression  can  be  compared  to  that  obtained  from  the  simplified  steady-state 
expression  (sss),  Eq.  B  Equating  k3,t  to  k3  and  using  other  rates  from  Tables  II  and  IV,  this 
simpilied  expression  becomes 


d[C2H 


=  3.4  x  1o"e-(40000  c*!/mcl*/RT)(C2H2]1  5  (cc/mole^sec"1 

dt  '  overall  1 


The  lower  A-factor  from  the  detailed  model  is  due  to  a  combination  of  the  reverse  of  R3,  alternative 
termination  reactions  (e  g  ,  R54),  and  the  fact  that  some  acetone  decomposes  via  a  chain  process 
rather  than  R5  (see  Reactions  21-23  in  Table  II). 

Higher  Temperature  Modeling 

Preliminary  modeling  has  been  performed  for  shock  tube  conditions  at  higher  temperatures.  The 
scenario  regarding  initiation,  termination,  fate  of  acetone  and  methyl  radicals  changes  somewhat. 
The  changes  include  a  greater  fraction  of  acetone  decomposing  via  radical  attack.  This  leads  to  the 
production  of  ketene  which  also  is  an  effective  radical  initiator.  The  resulting  sequence, 


H  +  CH3COCH3  «-*  CH2COCH3  +  H2  (22a) 


CH2COCH3  ~  CH2CO  +  CH3  (22b) 
CHrCO  -  CHj  +  CO  (24) 

CH2  +  C2H2  *-*  CjHj  +  H  (28) 

C3H3  4-  C3H3  <-*  C«He  (29) 


leads  to  the  net  production  of  only  one  methyl  radical  (the  H-atom  is  regenerated)  compared 
to  the  two  radicals  from  initiation  directly  through  acetone  decomposition  (Reaction  5).  This 
sequence  is  very  dependent  on  the  uncertain  high  temperature  chemistry  of  methylene  (CH2)  and 
propargyl  (C3H3)  radicals.  Besides  recombination,  the  latter  can  add  to  acetylene  to  form  a  C$H$ 
radical  whose  stability  is  enhanced  by  resonance.  Other  potential  complicating  features  are  that 
aome  species  formed  as  by-products  of  acetone  decomposition,  e.g.  cyclopentadienes,  have  weak 
C-H  bonds  because  of  the  resonantly-enhanced  stability  of  the  hydrocarbon  radical  (in  this  case 
cyclopentadienyl).  Consequently,  even  though  acetone  may  be  destroyed  early  during  reaction,  some 
of  the  products  will  continue  to  provide  radicals  and  maintain  the  decomposition  of  acetylene.  At 
high  temperatures,  vinyl  recombination  no  longer  dominates  radical  termination.  Methyl  radicals 
and  radical  adducts  formed  by  addition  to  acetylene  contribute  to  chain  termination. 


At  very  high  temperatures  (>2000K),  a  significantly  different  story  emerges.  Acetone  decomposes 
within  a  few  nanoseconds  to  methyl  radicals  which  rapidly  equilibrate  with  ethane  in  the  recombi¬ 
nation  reaction.  Eventually,  methyl  is  converted  to  two  H-atoms  via  the  sequence  Reaction  18,20, 
but  the  amount  of  H-atoms  is  small  compared  to  H-atom  production  from  pure  acetylene  pyrolysis 
as  initiated  by  cleavage  of  the  CjH-H  bond33.  In  fact  our  calculations  are  consistent  w  ith  the  results 
of  Frank  and  Just33  who  found  that  ‘purified’  and  unpurified  acetylene  produced  the  same  H-atom 
profiles  at  50ppm  acetylene  and  2400K.  In  contrast,  near  2000K  they  found  that  high  concentration 
runs  require  purification.  Again  this  result  is  identical  to  our  calculations.  In  fact,  our  model  pre¬ 
dicts  that  the  accelerated  rate  they  observe  at  1800-2000K  for  the  initiation  rate  can  be  explained 
by  the  presence  of  small  amounts  of  acetone  (—0.1%  in  acetylene) 

Sensitivity  and  Reaction  Pathway  Analyses 

A  sensitivity  analysis  has  been  performed  using  CHEMSEN  for  the  isobaric,  isothermal  conditions 
of  Munson  and  Anderson.  Results  of  these  calculations  are  reproduced  in  Table  VII  for  973  K 
after  10%  decomposition  for  the  six  reactions  most  sensitive  to  the  computed  results  for  each  of 
the  five  species  listed  in  Table  VII.  Rate  constants  for  several  of  these  reactions  were  determined  in 
this  work  Also,  a  reaction  pathway  analysis  was  performed  using  programs  developed  at  UTRC 
Reactant  decay  and  product  formation  as  determined  using  this  analysis  are  described  below 

Acetylene  was  calculated  to  decay  principally  by  three  reactions: 

H  +  C2H2  ~  C2H3  (2) 

C2H3  +  C2H2  «-*  n  —  C4H5  (3) 

n  -  C«Hs  +  C2H2  -4  C6H6  +  H  (33) 

The  net  rates  of  the  first  two  are  essentially  equal  whereas  the  rate  of  the  third  is  slower  because 
of  the  decay  of  n-C^Hs  to  vinylacetylene.  Vinylacetylene  is  formed  almost  entirely  by  H-atom 
elimination  from  n-C<H$  (Reaction  32).  Its  decay  is  due  predominantly  to: 


n  —  C«H&  +  C4H4  «-»  CgHg  +  H  (35) 


which  inturn  forms  styrene,  although  H-atom  abstraction  from  vinylacetylene  by  phenyl  or  vinyl 
(R41-42,  64-65)  also  contributes  to  vinylacetylene  decay.  The  formation  of  benzene  is  from  n-C«H$ 
addition  to  acetylene  (R33).  As  the  reaction  proceeds,  additional  benzene  is  produced  by  phenyl 
(from  R38)  abstracting  H-atoms  from  molecular  hydrogen,  vinylacetylene  or  ethene.  Ethene  arises 
mainly  from  H-atom  attack  on  styrene  (R38)  although  some  is  produced  by  abstraction  of  H-atoms 
from  hydrogen  and  vinylacetylene  by  vinyl  radicals  (R47  and  R65).  According  to  the  proposed 
model,  which  underpredicts  hydrogen  by  a  factor  of  five  to  ten,  hydrogen  is  formed  by  H-atom 
attack  on  styrene  (R36). 


In  order  to  fit  the  Munson  and  Anderson  data,  values  of  four  rate  constants  have  been  modified 
from  values  previously  proposed.  The  related  reactions  are 


C2 1 1 3  i"  CjHj  <— *  C4H6  (4) 

n  -  C<H<  —  C<H<  +  H  (32) 

n  —  C<H5  -*■  CjHj  — •  CeHe  +  H  (33) 

n  -  C*Hs  ■+■  C4H<  —  CgHg  +  H  (35) 


k4,  the  rate  constant  for  vinyl-vinyl  recombination,  was  selected  to  be  2  x  1013  cc/mole^ec  in  this 
work.  This  value  although  slightly  high  is  consistent  with  rate  coefficients  for  recombination  of 
other  low  molecular  weight  hydrocarbons  (see  Ref  34).  The  value  used  here  is  about  50%  larger 
than  the  expression  suggested  by  Ebert,  et  al35.  The  rate  constant  for  R32  was  based  on  the 
thermochemical  estimate  by  Weissman  and  Bensonse  but  the  A-factor  was  increased  by  60%  to 
enhance  the  fit  to  experimental  data  Equilibrium  calculations  give  a  reverse  rate  constant  of 
about  8  x  1012  cc/mole/sec  at  1000K.  This  value  is  about  a  factor  of  five  higher  than  the  H  + 
C2H2  reaction.  Errors  in  the  thermochemistry  and  possible  contributions  from  reactions  involving 
vinylidene  may  be  the  cause  of  the  high  rates  for  R32  and  R(-32)  as  required  by  this  analysis 
At  1000K,  the  irreversible  rate  constant,  kg3,  is  about  a  factor  of  four36,  two1123,  and  one  and 
one-half37  lower  than  values  previously  proposed.  The  major  difference  of  this  work  is  the  high 
activation  energy  (lOkcal/mole)  compared  to  other  evaluations:  3.737,  4 ,9(from  a  fit  to  reported23 
calculations,  6.930',  and  911.  Due  to  this  difference,  the  extrapolated  value  from  the  present  work 
is  essentially  identical  to  other  values23,37  near  1200-1300K.  kgs  =  1.1  x  1013e"BOO°/BT  as  used  in 
this  study  is  substantially  lower  than  the  high  value  suggested  by  Colket11  but  is  about  four  times 
larger  than  the  value  reported  by  Cole,  et  al.37.  As  discussed  previously,  the  i-C4Hs  radical  may 
play  a  role  or  C4H4-C4H4  reactions  may  contribute. 


Rate  constants  as  determined  in  this  work  are  generally  within  a  factor  of  four  of  other  values 
recently  proposed  and  experimental  data  for  these  reactions  is  extremely  limited:  usually,  rate 
expressions  have  been  estimated  or  determined  from  detailed  modeling.  Consequently,  we  believe 
that  the  evaluations  of  rate  constants  in  this  study  are  very  reasonable  and  lend  strong  support  for 
the  general  thesis  of  this  paper  regarding  the  importance  of  acetone. 


Uncertainties  in  the  analysis  include  possible  contributions  due  to  reactions  involving  vinylidene  and 
related  species.  As  discussed  previously,  we  are  unable  to  disprove  the  existence  of  such  mechanisms. 
Due  to  the  generally  good  agreement  between  the  modeling  (without  invoking  such  mechanisms) 
and  experimental  data,  we  believe  that  contributions  due  to  these  reactions  are  negligible;  instead, 
the  chain  median  -m  seems  to  dominate  in  most  cases.  It  is  possible  that  the  vinylidene  mechanism 
contributes  to  vinylacetylene  production  at  873K  (see  Fig.  3)  when  the  chain  underpredicts  the 
experimental  data.  Other  uncertainties  are  related  to  the  lack  of  a  model  for  methane  formation 
at  low  temperatures.  At  973K,  Munson  and  Anderson  found  that  about  1.2%  of  the  initial  carbon 
was  converted  to  methane.  Often  it  is  suggested  that  this  methane  formation  is  catalytic.  An 
alternative  and  speculative  mechanism,  described  in  this  paper  but  not  included  in  the  model 
involves  the  intermediate  formation  of  the  i-C^Hs  radical.  It  may  be  formed  by  H-atom  addition  to 
vinylacetylene  or  isomerization  of  the  excited  h-C4Hb  radical  after  its  formation  by  vinyl  addition 
to  acetylene.  The  formation  of  excited  radicals  via  radical  additions  to  acetylenes  is  another  area  of 
uncertainty  in  the  proposed  model.  A  limited  number  of  QRRK15  calculations  have  been  performed 
and  indicate  that  these  effects  should  be  included  at  elevated  temperatures.  An  example,  is  the 
system 


C:H?  -  CjH2  -  n  -  C4H5 
n  -  *—  n  -  C4H5 

n  -  C4H5  ~  C4H4  +  H 
n  -  C4H5  -  i  -  C4H5 


but  for  simplicity,  these  effects  have  not  been  included  in  the  present  analysis. 


Significant  uncertainties  exist  in  the  kinetics  and  mechanisms  related  to  the  odd-carbon,  acetylenic 
species.  This  belief  is  based  on  the  facts  that  very  little  experimental  data  or  estimates  on  such 
species  exist ;  large  uncertainties  exist  in  the  thermodynamics  of  the  odd-carbon  radicals;  and  predic¬ 
tions  of  profiles  of  odd-carbon  species,  although  qualitative  are  not  quantitatively  accurate.  Despite 
these  drawbacks,  the  major  conclusion  of  this  paper,  that  is  that  acetylene  pyrolysis  is  initiated  by 
a  small  amount  of  acetone,  should  not  be  altered. 


Conclusions 


The  major  conclusion  of  this  study  is  that  the  pyrolysis  of  pure  acetylene  has  not  been  investigated 
at  temperatures  below  2000K  The  very  high  C-H  bond  strengths  in  acetylene  prevent  a  rapid 
initiation  of  radicals;  instead  very  low  concentrations  of  impurities  are  the  source  of  radicals  which 
initiate  and  sustain  a  chain  decomposition.  Acetone,  a  contaminant  present  in  kinetic  studies 
of  acetylene  pyrolysis  has  been  shown  to  dominate  radical  initiation  The  overall  reaction  rate 
determined  both  from  a  simple  steady-state  analysis  and  detailed  chemical  kinetic  modeling  has  an 
activation  energy  near  40  kcal/mole  and  a  reaction  order  near  1.5.  The  overall  activation  energy 
is  consistent  with  previous  experiments  while  the  reacion  order  is  lower  than  the  usually  assumed 
value  of  two  A  detailed  review  of  the  literature  shows  in  fact  minimal  support  for  a  reaction  order 
of  two  while  experiments  are  reasonably  consistent  with  an  order  of  1.5.  A  few  rate  constants  have 
been  selected  to  provide  a  fit  to  the  Munson  and  Anderson  data.  Values  of  these  rate  constants  are 
in  general  consistent  with  values  previously  published.  The  modeling  makes  no  attempt  to  include 
reactions  involving  the  vinylidene  radical.  Although  there  is  no  proof  against  mechanisms  involving 
this  carbene,  there  appears  now  to  be  ample  evidence  in  support  of  a  chain  mechanism;  firstly, 
experimental  evidence  from  a  wide  variety  of  sources  and  conditions  support  the  presence  of  a  chain, 
and  secondly,  a  mechanism  initiated  by  thermal  decomposition  of  acetone  adequately  describes 
experimental  results.  Vinylidene-related  mechanisms  may  contribute  to  acetylene  pyrolysis  but  we 
believe  that  a  chain  mechanism  dominates  all  acetylene  pyrolysis  above  900K. 
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TABLE  II 

PROPOSED  SET  OF  REACTIONS  AND  RATE  COEFFICIENTS 
FOR  INITIATION  BY  ACETONE 
logic  =  logA  +  nlogT  -  E/R/T/2.303  * 


Reactions 


Forward 
Rate  Constant 


Reverse 
Rate  Constant 


logA 

n 

E 

logA 

n 

E 

Ref 

5 

CH3COCH3=CH3CO+CH3 

14.38 

0.0 

72.0 

11.55 

0.0 

-5.6 

14 

6 

CH3CO=CH3+CO 

12.48 

0.0 

16.7 

10.68 

0.0 

4.6 

a 

7 

CH3+C2H2=CH3CHCH 

11.79 

0.0 

7.7 

12.87 

0.0 

33.9 

38 

8 

CH3CHCH=C3H5 

13.15 

0.0 

36.0 

14.10 

0.0 

57.4 

15 

9 

C3H5+C2H2=C5H7 

12.00 

0.0 

8.0 

13.11 

0.0 

20.1 

est 

10 

C5H7-C-C5H6+H 

10.30 

0.0 

5.0 

0.00 

0.0 

0.0 

est 

11 

H+C3H4=CH3CHCH 

12.76 

0.0 

3.1 

12.52 

0.0 

38.2 

39 

12 

H+ALLENE=C3H5 

12.60 

0.0 

2.7 

13.13 

0.0 

60.9 

39 

13 

2C3H5=C3H6+C3H4 

12.70 

0.0 

0.0 

13.20 

0.0 

30.5 

est 

14 

C3H5+C2H3S=C5H8 

12.70 

0.0 

0.0 

15.62 

0.0 

82.5 

est 

15 

C3H5+H=C3H6 

13.60 

0.0 

0.0 

14.81 

0.0 

87.0 

est 

16 

CH3+C2H2=CH4+C2H 

12.40 

0.0 

35.0 

13.11 

0.0 

15.1 

est 

17 

2CH3=C2H6 

14.38 

-.4 

0.0 

16.31 

0.0 

87.1 

a 

18 

2CH3=C2H5+H 

14.90 

0.0 

26.5 

16.92 

0.0 

16.3 

a 

19 

C2H6+CH3=C2H5+CH4 

-0.26 

4.0 

8.3 

14.71 

0.0 

25.1 

a 

20 

C2H5=C2H4+H 

13.30 

0.0 

39.7 

12.50 

0.0 

1.9 

a 

21 

CH3+CH3COCH3- 

11.60 

0.0 

9.7 

0.00 

0.0 

0.0 

40 

CH4+CH2CO+CH3 

22 

H+CH3COCH3- 

13.28 

0.0 

6.4 

0.00 

0.0 

0.0 

39 

H2+CH2CO+CH3 

23 

C2H3+CH3COCH3- 

12.48 

0.0 

6.4 

0.00 

0.0 

0.0 

est 

C2H4+CH2CO+CH3 

24 

CH2CO=CH2+CO 

14.00 

0.0 

71.0 

11.65 

0.0 

-4.4 

c 

25 

CH3+H=CH2+H2 

14.86 

0.0 

15.1 

14.09 

0.0 

10.1 

41 

26 

CH3+CH2=C2H4+H 

13.30 

0.0 

0.0 

15.89 

0.0 

61.3 

41 

27 

2CH2=C2H2+H2 

13.30 

0.0 

0.0 

15.65 

0.0 

128.7 

41 

28 

CH2+C2H2=C3H3+H 

12.26 

0.0 

0.0 

13.36 

0.0 

11.5 

a 

29 

2C3H3-C6H6 

12.70 

0.0 

0.0 

0.00 

0.0 

0.0 

est 

30 

CH2CO+H=CH3+CO 

12.85 

0.0 

3.0 

11.87 

0.0 

37.0 

a 

31 

CH2CO+CH3=C2H5+CO 

12.30 

0.0 

3.0 

13.34 

0.0 

26.7 

est 

•  NOTES:  Units  for  A:  cc, moles, sec. 

Units  for  E:  kcal/mole. 

represents  forward  and  reverse  directions  included  in  model 
represents  forward  direction  only  included  in  model. 

PW  indicates  rate  evaluted  from  the  present  work, 
a  indicates  rate  constants  from  Ref.  19. 

c  rate  constant  a  factor  of  three  lower  from  that  reported  in 
Ref.  19. 

est  estimate  based  on  thermodynamics  and/or  analogous  reactions 


A-  19 


TABLE  III 


Experimental  Overall  Reaction  Orders 


Apparatus 

Temp  (K) 

Order  (n) 

Silcocks10 

static  reactor 

625-745 

l(het  ),2(hom.) 

Cullis  and  Franklin18 

static  reactor 

983 

2 

Munson  and  Anderson10 

flow  reactor 

873,923 

973-1073 

1.7±0.2 

1.5±0.1 

Palmer  and  Donnish17 

flow  reactor 

1333-1528 

l(het.),2(hom  ) 
[1.2  to  1.7]* 

Ogura3 

SPST 

1000-1670 

(-2  3) 

Colket 28 

SPST 

1100-1400 

1  6±0.2 

Aten  and  Greene42 

shock  tube 

1400-2500 

1.72±0.1 

Cundall,  et  al43 

shock  tube 

1500-2500 

r* 

Gay,  et  al44 

shock  tube 

1600-2400 

not  2 

Towell  and  Martin45 

flow  reactor 

1220-1400 

1.5“ 

reevaluated  in  present  study 
overall  order  for  formation  of  products 


TABLE  IV 

PROPOSED  SET  OF  REACTIONS  AND  RATE  COEFFICIENTS 
FOR  ACETYLENE  PYROLYSIS  AT  LOW  AND  INTERNED I ATE  TEKPERTURES 
log*  =  logA  +  nlogT  -  E/R/T/2.303  * 

Forward  Reverse 

Reactions  Rate  Constant  Rate  Constant 


logA 

n 

E 

logA 

n 

E 

Ref 

1 

2C2H2=n-C4H3+H 

13.00 

0.0 

67.0 

12.48 

0.0 

-2.0 

est 

2 

H+C2H2=C2H3 

12.74 

0.0 

2.5 

13.01 

0.0 

43.7 

49 

3 

C2H3+C2H2=n-C4H5 

12.04 

0.0 

4.0 

13.88 

0.0 

38.4 

11 

4 

2C2H3=C4H6 

13.30 

0.0 

0.0 

16.71 

0.0 

101.1 

pw 

32 

n-C4H5=C4H4+H 

14.20 

0.0 

41.4 

13.60 

0.0 

3.0 

pw 

33 

n-C4H5+C2H2-C6H6+H 

12.65 

0.0 

10.0 

0.00 

0.0 

0.0 

pw 

34 

C2H3+C4H4-C6H6+H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

b 

35 

C4H4+n-C4H5-C8H8+H 

13.08 

0.0 

5.0 

0.00 

0.0 

0.0 

pw 

36 

C8H8+H=C8H7  +H2 

14.60 

0.0 

7.0 

13.65 

0.0 

12.0 

11 

37 

C8H7  =C8H6+H 

12.85 

0.0 

45.0 

12.95 

0.0 

0.7 

est 

38 

C2H4+C6H5=C8H8+H 

11.27 

0.0 

2.1 

12.38 

0.0 

1.6 

d 

39 

C6H6+H=C6H5+H2 

14.40 

0.0 

16.0 

12.69 

0.0 

9.8 

a 

40 

C2H2+C6H5=C8H6+H 

12.00 

0.0 

4.0 

13.88 

0.0 

6.0 

b 

41 

C2H3+C4H4=C2H4+n-C4H3 

13.30 

0.0 

14.5 

13.23 

0.0 

13.4 

24 

42 

C4H4+C6H5=C6H6+n-C4H3 

13.30 

0.0 

14.5 

13.60 

0.0 

18.9 

24 

43 

L-C6H5=n-C4H3+C2H2 

14.00 

0.0 

36.0 

11.69 

0.0 

1.0 

11 

44 

C6H5=L-C6H5 

13.54 

0.0 

65.0 

9.91 

0.0 

1.1 

11 

45 

C2H3+C6H6=C2H4+C6H5 

13.30 

0.0 

14.5 

12.94 

0.0 

9.1 

e 

46 

C2H3+H=C2H4 

13.60 

0.0 

0.0 

15.55 

0.0 

105.1 

est 

47 

C2H4+H=C2H3+H2 

14.84 

0.0 

14.5 

13.49 

0.0 

13.8 

a 

48 

C2H3+H=H2+C2H2 

13.00 

0.0 

0.0 

13.34 

0.0 

63.1 

a 

49 

C4H4+H=n-C4H3 +H2 

13.90 

0.0 

14.5 

12.49 

0.0 

12.7 

a 

50 

C6H6=C6H5+H 

15.70 

0.0 

107.9 

13.39 

0.0 

-2.6 

47 

51 

C6H5+C6H6=C12H10+H 

11.80 

0.0 

11.0 

13.04 

0.0 

8.7 

b 

52 

2n-C4H5=C8H10 

12.70 

0.0 

0.0 

16.92 

0.0 

106.9 

est 

53 

n-C4H5+H=C4H6 

13.60 

0.0 

0.0 

15.44 

0.0 

107.9 

est 

54 

n-C4H5+C2H3=C6H8 

12.70 

0.0 

0.0 

16.33 

0.0 

104.0 

est 

55 

C4H6+H=n-C4H5+H2 

14.00 

0.0 

14.5 

12.76 

0.0 

10.9 

24 

56 

C4H6+C2H3«n-C4H5+C2H4 

13.30 

0.0 

14.5 

13.41 

0.0 

11.7 

e 

57 

C6H5+C2H3-C8H8 

12.70 

0.0 

0.0 

15.76 

0.0 

104.5 

est 

58 

i-C4H3+H2«C2H2+C2H3 

10.70 

0.0 

20.0 

10.58 

0.0 

17.8 

11 

59 

i-C4H3*C4H2+H 

12.00 

0.0 

49.0 

12.05 

0.0 

0.6 

11 

60 

n-C4H3«C4H2+H 

12.60 

0.0 

40.0 

12.95 

0.0 

-0.4 

11 

61 

n-C4H3*i-C4H3 

13.00 

0.0 

53.0 

13.30 

0.0 

61.0 

est 

62 

i-C4H3+H»C4H2+H2 

13.00 

0.0 

0.0 

13.65 

0.0 

56.0 

11 

63 

C4H4+H»i-C4H3+H2 

14.49 

0.0 

14.5 

13.38 

0.0 

20.7 

11 

64 

C6H5+C4H4=C6H6+i-C4H3 

13.30 

0.0 

14.5 

13.90 

0.0 

26.9 

e 

65 

C2H3+C4H4«C2H4+i-C4H3 

13.30 

0.0 

14.5 

13.54 

0.0 

21.4 

e 

66 

C4H4+i-C4H3-C8H6+H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

est 

67 

H2+M-2H+M 

12.35 

-.5 

92.5 

9.98 

0.0 

-13.1 

a 

•  MOTES:  see  notes  for  Table  II 

b  see  citations  in  Ref.  11.  (Reverse  rate  constants  nay  differ 
slightly  due  to  differences  in  the  thermodynamics  and  the 
temperature  range  of  the  fit  for  the  reverse  rate.) 
d  rate  constant  divided  by  2  from  that  cited  in  Ref.  11. 
e  see  text  . 


TABLE  V 

PROPOSED  SET  OF  REACTIONS  AND  RATE  COEFFICIENTS 
FOR  HIGH  TEMPERATURE  ACETYLENE  PYROLYSIS 
logic  =  logA  +  nlogT  -  E/R/T/2.303  • 

Forward  Reverse 


Reactions 

Rate  Constant 

Rate 

Constant 

log  A 

n 

E 

logA 

n 

E 

Ref 

68 

C2H2+H=C2H+H+M 

16.62 

0.0 

107.0 

15.21 

0.0 

-18.2 

b 

69 

C2H+C4H4=C2H2+i-C4H3 

13.60 

0.0 

0.0 

13.29 

0.0 

27.0 

b 

70 

n-C4H3=C2H2+C2H 

14.30 

0.0 

57.0 

13.42 

0.0 

0.8 

11 

71 

n-C4H3+H=C4H2+H2 

12.00 

0.0 

0.0 

12.95 

0.0 

64.0 

11 

72 

i-C4H3+H=C4H4 

13.78 

0.0 

0.0 

15.49 

0.0 

98.2 

est 

73 

C2H+H2=H+C2H2 

12.85 

0.0 

0.0 

13.65 

0.0 

20.8 

b 

74 

C2H+C2H2=C4H2+H 

13.60 

0.0 

0.0 

14.84 

0.0 

15.8 

b 

75 

C4H2=C4H+H 

14.89 

0.0 

120.0 

13.25 

0.0 

1.9 

b 

76 

C2H+C4H2=C6H2+H 

13.60 

0.0 

0.0 

15.04 

0.0 

15.5 

b 

77 

C4H+C2H2=C6H2+H 

13.30 

0.0 

0.0 

14.97 

0.0 

8.5 

11 

78 

C6H2=C6H+H 

14.89 

0.0 

120.0 

13.02 

0.0 

5.2 

b 

79 

C4H+H2=H+C4H2 

13.30 

0.0 

0.0 

14.34 

0.0 

13.8 

b 

80 

C6H+H2=H+C6H2 

13.30 

0.0 

0.0 

14.56 

0.0 

10.5 

b 

81 

C2H+C6H6=C6H5+C2H2 

13.30 

0.0 

0.0 

12.40 

0.0 

14.7 

b 

82 

C4H+C6H6=C6H5+C4H2 

13.30 

0.0 

0.0 

12.63 

0.0 

7.6 

b 

83 

C2H3+C4H2=C4H4+C2H 

13.48 

0.0 

23.0 

13.47 

0.0 

3.1 

11 

84 

C4H4=C2H+C2H3 

16.00 

0.0 

105.0 

13.37 

0.0 

-16.2 

11 

•  NOTES:  Units  for  A:  cc, moles, sec. 

Units  for  E:  kcal/mole. 

«=«  represents  forward  and  reverse  directions  included  in  model. 

represents  forward  direction  only  included  in  model. 

PW  indicates  rate  evaluted  from  the  present  work, 
b  see  footnote  b  in  Table  IV. 

est  estimate  based  on  thermodynamics  and/or  analogous  reactions. 
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TABLE  VI 

SELECTED  THERMODYNAMICS  AT  300K 


SYMBOL 

NAME/structure 

S° 

eu 

AHf 

kcal/mole 

CuHio 

biphenyl 

93.6 

43.6 

CgHio 

octatetraene 

90.4 

53.2 

c8h8 

styrene 

82.6 

35.3 

c8h7 

CcH5CCH2 

83.8 

83.2 

c8h6 

phenylacetylene 

76.4 

75.2 

c6h8 

hexatriene 

79.4 

39.6 

c6h6 

benzene 

65.2 

19.8 

c6h5 

phenyl 

69.1 

78.2 

c6h2 

triacetylene 

71.1 

169.7 

c6h 

CcH 

74.3 

233.2 

c4h6 

1,3-butadiene 

66.6 

26.1 

n-C4H5 

H2CHCHCH 

69.1 

82.5 

c4h4 

vinylacetylene 

66.1 

69.4 

i-C4H3 

HCCCCH2 

68.1 

116.1 

n-C4H3 

HCCCHCH 

69.4 

124.1 

c4h2 

diacetylene 

59.9 

111.7 

c4h 

C4H 

62.8 

179.0 

c2h3 

vinyl 

54.5 

65.7 

c2h 

ethynyl 

49.6 

128.5 

c5h8 

1,4-pentadiene 

76.6 

18.7 

CsH7 

1,4-pentadien-l-yl 

79.3 

78.2 

c-C8H6 

cyclopentadiene 

64.6 

32.0 

CH3CO 

acetyl 

62.0 

-6.0 

c2h5 

ethyl 

54.6 

26.4 

CH3CHCH 

CH3CHCH 

67.5 

61.6 

c3H5 

ally! 

62  I 

39  4 

c3h4 

methylacetylene 

59.3 

44.4 

c3h3 

propargyl 

60.0 

81.5 

ch2 

methylene 

46.4 

92.4 

ch2co 

ketene 

57.8 

-12.4 

note:  Thermodynamics  of  species  not  listed  in  this  table  are  essentially 
identical  to  those  listed  in  Ref.  25. 


A-  23 


L. 


Normalized  Sensitivity  Coefficients 


dy 

dk 


x 


k* 

y 


2)  H  4  CjH2  -  C2H3 

3)  C2H3  4  C2H2  ♦— *  n-C^Hs 

4)  C2H3  4  C2H3  ♦— *  C«Hc' 

5)  ch3coch3  -  ch3  +  ch3co 

23)  C2H3  4  CH3COCH3  —  C2H<  4  CH2CO  4 

32)  n-C^Hs  *— *  C4H4  4  H 

33)  n-C<Hs  4  C2H2  — *  CeHe  +  H 

34)  C2H3  +  C4H4  -  C6Hg  4  H 

35)  n-C<H5  4  C4H4  -  CfH8  4  H 

36)  C*Hg  4  H  «-*  CgHr  4  H2 

4l)  C2H3  4  C4H4  *—*  C2H<  4  n-C4H3 
45)  C2H3  4  CgHo  <-*  C2H4  4  CgHs 


c2h2 

C2H4 

C4H4 

CeH6 

h2 

-0.014 

0.06 

0.13 

0.09 

-0.52 

-0.030 

0.10 

0.18 

0.21 

0.61 

0.025 

-0.42 

-0.11 

-0.23 

-0.55 

-0.050 

0.85 

0  22 

0  47 

1.11 

- 

0.20 

-0.01 

-0.01 

-0.02 

-0.020 

0.35 

0.68 

-0.14 

0.64 

-0  035 

0.13 

-0.13 

0.73 

-0.04 

- 

-0.02 

-0.12 

0  10 

-0.01 

- 

-0.01 

-0.30 

-0.10 

0.54 

- 

-0.08 

-0.02 

-0.04 

0.54 

- 

0.18 

- 

- 

- 

- 

0.22 

- 

- 

- 

*  As  calculated  by  CHEMSEN  at  10%  decomposition  of  acetylene,  973K,  20%  acetylene  in  one 
atmosphere  of  argon. 


LIST  OF  FIGURES 


1.  Comparison  of  Production  of  Vinylacetylene  and  Benzene  in  a  Single-Pulse  Shock  Tube.  3.7% 
acetylene,  2040ppm  biacetyl,  3.55%  acetylene/ I500ppm  biacetyl.  Squares  -  vinylacetylene, 
triangles  -  benzene.  Total  pressure  ~  eight  atmospheres,  dwell  times  ~  500-700  microseconds. 

2.  Comparison  of  Data  and  Model  Predictions  of  Acetylene  Decay.  Symbols  and  aolid  lines  are  data 
and  best  fits  from  Munson  and  Anderson.  Dotted  lines  are  model  predictions  from  the  present  study. 

3.  Comparison  of  Data  and  Model  Predictions  of  Vinylacetylene  Formation  and  Decay  (see  caption 
of  Fig.  2.) 

4.  Comparison  of  Data  and  Model  Predictions  of  Benzene  Formation,  (see  caption  of  Fig  2.) 

5.  Comparison  of  Data  and  Model  Predictions  of  Ethylene  Formation,  (see  caption  of  Fig  2  ) 

6  Single-pulse  shock  tube  data  for  4.9%  acetylene,  lOOppm  acetone  in  eight  atmospheres  of  argon 
Only  odd-carbon  species  are  shown  including  acetone,  allene  plus  methyl  acetylene,  four  Cj-species 
(one  of  two  dominant  peaks  idendified  as  cylcopentadiene),  and  toluene  Top  line  is  molar  sum  of 
the  acetone  concentration  and  one-half  of  all  the  odd-carbon  species. 

7.  Comparison  of  Data  and  Model  for  Vinylacetylene,  Diacetylene,  and  Hydrogen  for  Pyrolysis  of  Five 
Percent  Acetylene.  Single-pulse  shock  tube  data  is  from  Ogura.  Solid  Lines  are  model  predictions. 
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SOME  THOUGHTS  ON  MODELING  PRE-PARTICLE  CHEMISTRY 

by 

Meredith  B.  Colket,  III 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


Within  the  last  decade,  pre-particle  kinetics  has  been  identified 1-3  as  an  important  step  in  soot 
formation.  This  chemistry  limits  particle  inception  rates  which  in  turn  limits  the  total  number  of 
soot  particles.  Since  the  number  of  particles  limit  the  total  surface  area  available  for  surface  growth 
and  since  particle  sizes  are  typically  limited  to  300  to  500A,  the  initial  chemistry  leading  to  soot 
inception  limits  both  the  rate  of  growth  (i.e.,  mass  addition)  and  the  total  amount  of  soot  formed  in 
a  given  environment.  This  very  important  conclusion  has  led  to  a  substantial  effort  directed  towards 
understanding  pre-particle  chemistry  and,  more  generally,  the  pyrolysis  of  hydrocarbons.  Chemistry 
involving  neutral  species  has  been  a  focus  of  this  research  since  recent  results4-5  indicate  that  ions 
play  at  most  a  minor  role.  This  approach6  seems  quite  reasonable  until  such  time  as  non-ionic 
reactions  can  be  shown  to  be  insufficient  for  prediction  of  particle  inception  rates. 

Research  in  pre-particle  kinetics  has  resulted  in  the  identification  of  several  possible  mechanisms 
for  the  formation  of  a  single  aromatic  ring  These  reactions  involve  addition  of  radical  intermediates 
to  acetylenic  (and  sometimes  olefinic)  compounds.  Many  of  these  are  summarized  in  Table  I  where 
only  initial  reactants  and  forward  rate  constants  are  reported.  Typically,  these  reactions  include  the 
formation  of  a  linear  adduct  followed  by  cyclization,  and  loss  (or  addition)  of  an  H-atom.  Energetics 
of  some  of  these  reactions  have  been  described  recently T . 

Table  I  presents  a  wide  range  of  possibilities  for  ring  formation.  Generally,  for  a  given  fuel  and 
temperature  range,  a  particular  sequence  will  dominate.  For  example,  during  pyrolysis  of  allene, 
the  Cs  reactions  are  most  important;  in  1,3-butadiene  flames,  reactions  involving  the  butadienyl 
radical  dominate;  and  during  acetylene  pyrolysis,  steps  involving  radical  addition  to  acetylene  are  the 
fastest.  Temperature  is  also  important.  At  low  temperatures  (<1500K)  during  acetylene  pyrolysis8, 
rings  are  formed  principally  via  n-  CtHs  addition  to  acetylene  yet  at  higher  temperatures  T1-C4H3 
addition  to  acetylene  dominates.  The  temperature  dependence  of  mechanistic  pathways  to  rings  is 
simply  explained  by  the  temperature  dependence  of  radical  concentrations.  At  low  temperatures,  n- 
C4H&  is  formed  rapidly  via  H-atom  addition  to  acetylene  to  form  vinyl  which  then  adds  to  acetylene; 
whereas  at  high  temperatures,  the  favored  H-atom  abstration  from  acetylene  forms  ethenyl  {CjH) 
which  then  adds  to  acetylene  thereby  forming  n  -  C4H1.  In  a  similar  matter,  fuel  dependence  of 
mechanisms  is  explained  by  the  principal  species  present  during  pyrolysis  or  oxidation.  The  net 
conclusion  is  that  there  is  no  unique  pathway  to  ring  formation. 

This  conclusion  is  discouraging  when  considering  that  formation  of  the  first  ring  is  only  the  first 
step  in  a  comprehensive  model  for  soot  formation.  The  potential  complexity  of  such  a  model  could 
be  overwhelming  Consider,  as  an  example,  the  first  attempt  at  detailed  modeling9  which  included 
619  reversible  reactions  to  describe  soot  formation  during  pyrolysis  of  acetylene.  In  subsequent 
publications10-11,  the  reaction  complexity  has  increased  further  in  order  to  include  additional  fuels 
and  oxidation  steps.  Hopefully,  there  is  a  way  to  simplify  such  models.  In  fact,  there  is  some 
evidence  supporting  such  an  expectation  in  flames  or  during  oxidation.  During  acetylene  oxidation 
in  a  shock  tube10,  the  dominant  sequence  for  ring  formation  was  found  to  be 
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H  +  Cji/j  <— >  C%Hs 
C%H $  -t-Cjffj  «-»  n  —  C4//5 
n  —  C4//5  +  Cjf/j  *-*  t  —  C6Ht 
£  —  C(,Hy  *—>  c  —  CqHj 
c  -  CqHi  ~  Cetf6  +  H 

Recent  measurements12  on  a  methane/air  diffusion  flame  are  consistent  with  this  conclusion. 
Other  encouraging  results  have  been  obtained  on  premixed  flames.  For  example,  by  controlling 
temperature  while  using  a  wide  range  of  fuels,  sooting  points  were  correlated  independent  of  fuel 
structure18.  Thus,  it  is  probable  that,  in  each  of  these  flames  with  different  fuels,  soot  is  formed  by 
the  same  mechanism.  Similar  conclusions  were  obtained14  when  soot  formed  in  a  lightly-sooting, 
ethylene  premixed  flame  was  found  to  be  essentially  identical  to  that  formed  in  a  flame  in  which 
a  portion  of  the  ethylene  fuel  was  substituted  with  toluene.  Toluene  is  normally  considered  to  en¬ 
hance  significantly  soot-forming  characteristics.  In  diffusion  flames  and  in  heavily  sooting  premixed 
flames,  soot  formation  is  dependent  on  the  structure  of  the  fuel.  However,  such  fuel-dependent 
soot  formation  may  simply  be  an  enhancement  of  the  underlying  fuel-independent  Boot  formation 
mechanism.  This  acceleration  may  arise  principally  from  an  enhancement  in  the  initial  rate  of  ring 
formation  while  mechanisms  related  to  the  growth  of  PAH’s,  soot  inception,  or  soot  growth,  etc., 
remain  unchanged.  There  is  some  encouragement  to  assume  this  optimistic  viewpoint.  Recently, 
sooting  tendencies  for  several  fuels  have  been  correlated15  with  benzene  production  as  measured 
using  a  time-of-flight  mass  spectrometer  coupled  to  a  shock  tube.  Recent  measurements  in  a  single- 
pulse  shock  tube  are  consistent  with  this  observation.  Thus,  the  influence  of  fuel  structure  on  soot 
formation  may  lie  (in  the  first  order)  in  the  propensity  of  a  given  fuel  to  initiate  the  formation  of 
aromatic  rings.  If  this  is  true,  then  fuel-type  effects  can  be  simulated  in  a  soot  model  by  includ¬ 
ing  reactions  (specific  to  the  particular  fuel)  for  the  formation  of  aromatics  with  one  (or  possibly 
two)  ring(s).  The  remainder  of  the  model  would  be  an  ’underlying’  mechanism  which  is  similar  to 
that  occurring  in  other  flame  or  pyrolytic  systems.  This  approach  is  similar  to  that  utilized  in  the 
extension  of  the  acetylene  model0  to  other  fuels  and  to  oxidation10-11. 

Conversion  of  identifiable  hydrocarbon  species  into  nascent  soot  particles  is  a  research  topic 
of  which  very  little  is  presently  known.  Due  to  our  inability  to  measure  dynamically  and  quan¬ 
titatively  species/par  tides  in  size  ranges  of  5-20A,  it  seems  unlikely  that  speculative  models  can 
be  confirmed  experimentally,  at  least  for  the  near  future.  Consequently,  we  will  have  to  rely  on 
thought  experiments  and  computer  modeling  to  test  out  various  theories  of  soot  inception.  Recent 
proposals  include  condensation  of  high  boiling  point  compounds16  and  formation  of  concentric,  im¬ 
perfect  spheres  consisting  of  carbon  atoms  bonded  in  pentagonal  and  hexagonal  configurations17. 
Some  authors10  have  avoided  the  uncertainty  of  nucleation  by  simply  allowing  large  poly-nuclear 
aromatics  to  continue  growing  by  sequential  addition  of  acetylene  (and  H-atom  abstractions).  In 
this  model  there  is  no  upper  limit  to  size.  Other  authors*'18  have  suggested  that  nucleation  occurs 
by  collision  of  two  (or  three)  pre-particle  species,  thus  forming  a  three-dimensional  aggregate.  Data 
on  particle  inception  rates18  have  then  been  used  to  estimate  that  these  pre-particle  species  are 
relatively  small  polycyclic  aromatics  comparable  in  size  to  pyrene.  This  proposal  although  yet  to  be 
proven  is  obviously  attractive  due  to  the  potential  of  simplifying  an  otherwise  extremely  complex 
system. 

Is  it,  or  will  it  be  possible  to  construct  a  comprehensive  model  for  pre-particle  kinetics  and 
soot  formation?  Based  on  the  above  analysis,  the  author  has  increased  hope  that  such  a  task  can 
be  accomplished.  In  summary,  the  reasons  for  this  optimism  include: 


(1)  the  probable  existence  of  an  ’underlying’  mechanism  for  pre-particle  chemistry  which  for  flames 
is  independent  of  fuel; 

(2)  the  hope  that  fuel-type  dependence  may  be  limited  to  enhancement/suppression  of  the  forma¬ 
tion  of  the  first  one  or  two  rings; 

(3)  the  fact  that  a  mechanism  has  been  proposed10  which  offers  a  starting  or  reference  point  for 
revised  models  of  pre-particle  chemistry;  and 

(4)  the  existence  of  several  reasonable  proposals  for  soot  inception,  i.e. ,  conversion  of  gas-phase 
species  into  solid  particles, 

Much  remains  to  be  done  including  obtaining  new  data  on  thermochemistry,  kinetics,  and 
mechanisms.  In  addition,  methods  for  conversion  of  a  detailed  chemical  kinetics  model  and  a  model 
describing  growth/  dynamics/  agglomeration/  and  oxidation  into  a  form  usable  to  designers  of 
engines  must  still  be  determined 
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TABLE  I 

Recent  Proposals  for  the  Formation  of  the  First  Aromatic  Ring 


Initial 

Forward  Rate  Constant1 

Reference 

Reactants 

A  _ 

E 

CH3  +  RH... 

- 

- 

Benson 

n  —  C+Hz  +  C3H3 

5.0*10n 

0.0 

Colket(1987) 

1.0*1013 

0.0 

Frenklach,  et  al.(1985) 

7.5*10n 

5.4 

Westmoreland 1  ( 1 986) 

n  —  CtHs  -t-  Ci  H? 

6.5*101J 

9.0 

Colket(1987) 

3.2*10u 

3.7 

Cole,  et  al.(1983) 

6.7*10n 

4.9 

Westmoreland J  ( 1 986) 

1.0*1013 

0.0 

Frenklach,  et  al.(1985) 

4.0*10lJ 

6.9 

Weissman  and  Benson  (1984) 

n  —  C\H$  +  C4H4 

3.2*10u 

0.6 

Cole,  et  a). (1983) 

7.9*1013 

3.0 

Colket(1987) 

n  —  CiHs  -f  C3H4 

3.2*10" 

3.7 

Cole,  et  al  (1983) 

n  —  CiHf,  +  C4H2 

3.2*10" 

1.8 

Cole,  et  al.(1983) 

CiHs  +  CiHe 

3.2*10" 

3.0 

Weissman  and  Benson  (1984) 

C1H3  +  C4H4 

4.0*10" 

0.0 

Colket(1987) 

C3Hz  +  C3H4 

2.2*10" 

2.0 

Wu  and  Kern(1987) 

Cstfs  +  Cjffs 

3.0*10" 

0.0 

Wu  and  Kern(1987) 

1.  k=Aexp(-E/RT) 

cm3/moles/sec  (E 

in  kcal/mole) 

2.  Fit  for  one  atmosphere,  500<T<1500K 


Kinetic  Mechanist  for  Pvrolysis 
of  Acetylene  Near  1000K 

by  M.  B.  Colket,  III 

United  Technologies  Research  Center,  East  Hartford,  CT  06108 

INTRODUCTION 

Recently  a  mechanism  for  acetylene  pyrolysis  was  proposed  to  describe 
experimental  data  obtained  from  a  single-pulse  shock  tube  (1)  over  the 
temperature  range  of  1100  to  2000K.  The  kinetic  mechanism  is  similar  to  those 
proposed  previously  (2,3)  and  includes  mechanisms  to  describe  formation  of 
vinylacetvlene ,  benzene,  and  phenylacetylene .  Subsequently,  a  subset  of  this 
mechanism  was  modified  and  used  (4)  to  match  decay  profiles  of  acetylene  as 
well  as  product  formation  for  pyrolysis  data  obtained  in  a  flow  reactor  (5)  at 
temperatures  of  873  to  1173K. 

It  is  the  objective  of  this  study:  (a)  to  extend  the  mechanism  (4) 
include  formation  of  higher  molecular  weight  species,  with  a  detailed 
discussion  of  growth  from  one  to  two  rings;  (b)  to  discuss  the  possible 
of  the  'odd'  radical,  i-C^H^;  (c)  to  discuss  some  uncertainties  with 

modeling  PAH  formation. 

Description  of  Model 

The  kinetic  model  used  in  this  study  is  listed  in  Table  I  and  contains  89 
reactions  and  43  species.  CHEMKIN  and  LSODE  were  used  for  integration  of  the 
rate  equations.  For  modeling  of  the  flow  reactor  data,  constant  temperature 
and  pressure  conditions  were  imposed  whereas,  for  the  shock  tube  data,  a  shock 
tube  code  modified  (1)  to  simulate  quenching  in  a  single-pulse  shock  tube  was 
utilized. 

Thermodynamic  parameters  (except  for  those  of  vinyl  radicals)  are  the 
same  as  those  used  and  reported  previously  (1)  for  the  lower  molecular  weight 
species.  Parameters  for  species  with  molecular  weights  above  100  AMU  were 
obtained  from  Stein  (6)  and  are  believed  to  be  identical  to  those  used  in 
Frenklach,  et.  al .  (3)  Reference  data  for  vinyl  radicals  and  heavier  species 
are  reported  in  Table  II.  Stein's  thermodynamic  data  for  vinyl  radicals  were 
adopted  since  his  heat  of  formation  (68.4  kcal/mole)  is  higher  than  that  used 
by  Colket  (1)  (65.7  kcal/mole)  and  is  closer  to  some  recent  determinations. 

In  addition,  Stein's  reference  value  for  entropy  is  higher  (by  2  eu)  than  that 
used  by  Colket.  Fortunately,  the  differences  in  values  are  nearly  offset  when 
determining  equilibrium  constants  (In  K  -  -  (AH-TAS)/RT)  and  changes  in 
chemical  kinetic  modeling  results  are  minimal. 

The  chemical  kinetic  model  differs  from  that  used  previously  (4)  in  that 

C„H,  C, H„ ,  C.H,  i-C.H,,  C.H.  and  C,H  and  associated 
2  _  .  4  2  4  .  ,  4  3  6t2  ...6  .  ,  „ 

reactions  were  included.  In  addition,  species  and  reactions  related  to 

formation  of  polycyclic  aromatics  were  also  added  Of  particular  note  is  a 

modification  in  the  dominant  bimolecular  initiation  step.  Previously  the 
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C.H.  +  C0H„  ^  n-C. H,  +  H 
i  l  l  t  4  3 

was  used  with  a  rate  constant  of  log  (k/sec  -  14.54  -  68000/4. 58T.  In 
the  present  study  the  reaction 

C-H-  +  C.H.  si  i-C.  H,  +  H 
i  l  l  i  4  3 

which  has  a  lower  endothermicity  was  employed  with  a  rate  of  log  (k/sec  - 
14 . 54 - 60500/4 . 58T .  This  larger  rate  constant  for  initiation  is  required  to 
counterbalance  termination  steps  omitted  previously. 


Formation  of  Fused  Rings 

This  work  was  guided  substantially  by  the  significant  contributions  of 
Bittner,  Howard,  and  Palmer  (7)  and  of  Frenklach,  et.  al  (3).  An 
important  conclusion  of  the  former  work  is  that  ring  growth  is  dominated  by 
addition  of  aryl  radicals  to  triple  bonds,  followed  by  addition  of  acetylene 
to  the  resultant  vinyllic  (aromatic)  radical,  cyclization,  and  loss  of  an 
H-atom. 

i.e.  , 

0  +  C,H0  =i  0CHCCCH 
4  2 

CCHCCCH  +  C2H2s:  0CHC(CCH)CHCH 
0CHC(CCH)CHCH  si  QCT  CCH  -  H+ 

If 

Other  than  reverse  processes,  the  main  competitive  process  which  can 
inhibit  ring  growth  is  the  thermal  decomposition  of  the  vinyllic  adducts 
(i.e.,  loss  of  H-atoms).  Bittner,  et.  al ,  reached  specific  conclusions 
regarding  the  importance  of  the  thermal  decomposition  of  the  vinyllic  adduct 
relative  to  its  addition  to  acetylene.  However,  these  conclusions  may  have  to 
be  re-examined,  since  calculations  were  based  on  low  pressure  flame  conditions 
rather  than  those  of  a  high  pressure  combustor. 

Frenklach,  et.  al ,  included  at  least  six  separate  reaction  sequences  for 
growth  from  a  single  to  a  fused  ring.  One  reaction  sequence,  however,  was 
found  to  dominate  although  a  second  played  a  minor,  but  contributing  role. 

The  second  (minor)  reaction  sequence  is  comparable  to  that  proposed  by 
Bittner,  et.  al  (see  above.)  although  phenyl  radicals  add  to  acetylene,  rather 
than  diacetylene.  The  dominant  reaction  sequence  was  found  to  be  initiated  by 
H-atom  abstraction  from  the  ortho  position  on  phenylacetylene ,  followed  by 
acetylene  addition  and  cyclization. 

0CCH  +  C2H2  S  0(CCH)CHCH 
0(CCH)CHCH  ^  0$ 

The  resultant  aryl  radical  can  subsequently  add  to  triple  bonds  to  continue 
growth  to  higher  order  polycyclic  aromatics.  Due  to  the  lack  of  experimental 
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rate  data  for  reactions  involving  abstraction,  addition,  cyclization,  or  ring 
fracturing  of  polycyclic  aromatics,  Frenklach,  et.  al ,  selected  generic  rate 
constants  for  classes  of  reaction. 


Rate  constants  for  reactions  involving  aromatic  species  as  listed  in 
Table  1  were  selected  using  the  same  technique  of  Frenklach,  et.  al ,  i.e., 
rate  constants  for  a  class  of  reactions  are  equated.  Previously  (3) ,  values 
for  several  of  these  classes  were  assigned  since  experimental  data  was  sparse. 
In  this  study,  we  have  where  possible  updated  the  reaction  rate  data  to  be 
consistent  with  recently  available  determinations. 


Predictions  from  the  model  (at  973K,  20%  acetylene,  and  one  atirosphere ) 
are  compared  to  the  data  from  Munson  and  Anderson  (5)  in  Fig.  1  for  acetylene 
decay  and  production  of  benzene.  In  Fig.  2,  predictions  of  production  of 
styrene,  naphthalene,  and  phenanthrene  are  shown.  Although  there  is  no 
experimental  data  from  Munson  and  Anderson  for  these  species,  the  final 
concentration  of  naphthalene  is  similar  to  other  results  (8)  near  10GC-1100K. 
The  predicted  value  for  styrene  is  about  a  factor  of  ten  high. 


Initially,  only  the  dorr.inant  and  minor  mechanism  (as  identified  by 
Frenklach,  et.  al)  were  included  in  the  reaction  sequence,  specifically 
Mechanisms  I  and  II: 


CC2H  +  r  ^  :C2H  +  RH 


(A) 


;c2h  +  c2h2 


:(c2h)chch 


C(C2H)CHCH  Si  A2 


and 


II 


CC2H  +  H 


c  4  c2h2 


CHCK 

CCHCH 

:  T CHCHCHCH 

0CHCHCHCH  si  A2H  +H 


(B) 

(C) 


£CHCH  4  C  H2 


where  A2  represents  the  1 -naphthyl  radical  and  A2H,  naphthalene. 


The  net  contribution  of  the  second  sequence  to  the  formation  of  A2H 
(or  A2)  dominates  over  that  of  the  first  by  several  orders  of  magnitude. 

For  the  specific  conditions  considered  in  this  study,  this  fact  can  be  easily 

explained. 


At  the  low  temperature  (=»1000K)  and  the  high  initial  concentrations  of 
acetylene,  forward  reaction  rates  can  be  shown  to  dominate.  Ignoring  the 
contribution  of  Reaction  C  and  taking  the  radical,  R,  in  Reaction  A  to  be  an 
H-atom,  then  the  relative  rate  is  simply  the  rate  of  H-atom  abstraction  from 
the  ring  by  H-atoms  relative  to  the  rate  of  H-atom  addition  to  the  acetylenic 
group  in  phenylacetylene 


The  value  (9)  of  k  (C^H^  +  H  -•  4 
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at  1000K  is  about  10  ^  cc/rr.ole  -  sec .  The  rate  constant  for  H-atom  addition 
to  acetylene,  according  to  Ellul,  et.  al  (10),  extrapolates  to  2  x  101 
cc/mole-sec  at  1000K.  Using  these  rate  constants  as  estimates  for  k  and 
k&,  respectively,  Mechanism  II  is  approximately  200  times  faster  than  the 
first.  The  dominance  of  the  Mechanism  II  becomes  more  apparent  when  the 
contribution  of  Reaction  C  is  considered,  since  Reaction  C  is  two  to  three 
orders  of  magnitude  faster  than  Reaction  B.  Consequently,  the  second  sequence 
is  more  than  four  orders  of  magnitude  faster  than  the  first.  This  is  opposite 
to  the  trend  observed  by  Frenklach,  et.  al .  This  strong  discrepancy  can  be 
explained  by  a  combination  of  (a)  the  higher  temperatures  of  their  study  which 
enhance  thermal  decomposition  of  the  radical  adducts;  (b)  the  lower  partial 
pressures  of  acetylene  in  their  study  (40  torr  vs.  150  torr  in  this  studv) 
which  reduce  the  rate  of  radical  addition  to  acetylene;  and  (c)  the  use  by 
Frenklach  of  a  high,  temperature  independent  rate  constant  for^geaction  A 
(with  R  as  H-atom).  Their  rate  constant  was  selected  tj^be  cc/mole-sec 

whereas,  Kiefer,  et.  al ’ s  expression  (9)  gives  1.6  x  101  cc/mole-sec  at 
1600  K. 

Consequently,  a  reanalysis  of  dominant  reactions  occurring  in  a  practical 
device  should  be  performed  for  the  specific  ambient  conditions.  A  preliminary 
analysis  of  the  kinetic  model  indicates  that  the  Mechanism  II  will  dominate  at 
temperatures  of  1500-1700K  for  high  pressure  combustors,  in  which  local 
acetylene  concentrations  may  be  at  least  an  order  of  magnitude  larger  than 
considered  in  this  and  and  previous  studies. 

Also  found  to  contribute  a  minor  but  significant  role  at  973K  are  the 
overall  reactions 

0  +  C.H.  -  A2H  +  H 

4  U 


0  C.H  +  i-C  H-  -  A2C_H  +  H 
2  4  3  2 

both  of  which  require  H-atom  shifts  prior  to  cyclization.  The  second  of  these 

could  be  part  of  a  very  attractive  sequence  (as  shown  in  Fig.  4)  depending  on 

the  concentration  of  i-C.H,  radicals. 

4  3 


Role  of  i-C.H. 

4  3 

It  is  recognized  that  thermochemistry  plays  a  significant  role  in  the 
ability  to  model  the  above  processes.  Frenklach,  et  al.  (11) 
have  demonstrated  quantitatively  that  uncertainties  in  thermochemistry 
drastically  affect  computed  results.  An  interesting  and  possibly  important 
thermochemical  aspect  of  acetylene  pyrolysis  is  due  to  the  rather  large 
difference  between  the  heats  of  formation  of  the  two  isomers  of  C.H,, 
i.e.,  HCCCHCH  and  HCCCCH2 .  The  separation  in  this  work  was  taken  to  be  10 
kcal/mole,  although  Stein  (6)  and  Bittner  (12)  give  8  and  15  kcal/mole, 
respectively.  Due  to  its  relative  stability,  the  isomer  with  the  unpaired 
electron  on  the  secondary  carbon  atom  (i-C^H^)  becomes  a  dominant  radical 
in  the  acetylene  system.  The  situation  is  exacerbated  if  an  isomerization 
step  (i-C^H^  ~  n-C^H^)  is  not  included.  Depending  on  temperature,  the 
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concentration  of  i-C.H.,  is  two  to  three  orders  of  magnitude  higher  than 
that  of  its  isomer  andean  order  higher  than  that  of  the  vinyl  radical . 
Consequently,  it  is  logical  that  i-C^H^  plays  a  significant  role  in 
termination,  and  quite  possibly  in  ring  formation  and  growth. 

The  importance  of  this  radical  to  chain  termination  or  to  ring  formation 
and  growth  is  dependent  critically  on  (a)  the  thermodynamics  of  the 
isomers;  (b)  the  isomeration  rate  (equated  in  this  work  to  a  rate 
suggested  (13)  for  the  i -propyl  —  n-propyl  isomeration);  and  (c)  rate  constants 
for  reactions  forming  and  destroying  i-C^H..  Reactions  which  dominate 
formation  of  i-C,H.  include  H-atom  abstraction  from,  vinylacetylene  by 
phenyl,  vinyl  ana  H-atom.s.  The  principal  destruction  mechanism  in  the  present 
study  is  the  isomerization  to  n-C^H  . 


Uncertainties 


Significant  uncertainties  in  both  rate  constants  and  mechanisms  still 
exist.  The  initiation  step  ir.  acetylene  pyrolysis  has  been  a  matter  of 
unresolved  dscussior.  for  nearly  thirty  years.  There  is  perhaps  just  as  much 
uncertainty  in  the  termination  step(s).  In  this  work,  recombination  of 
C  H3  +  C^H,,  H  +  C^H,,  i-C^H3  +  i-C^H  and  H  +  i-C^H 

all  contributed  to  termination,  yet  there  is  essentially  no  data  available  on 
the  absolute  value  of  their  rate  constants  (although  some  information  on 
reverse  reactions  is  available) .  The  addition  of  H-atoms  to  acetylene  is 
critically  important  to  this  mechanism,  yet  there  is  no  data  on  this  reaction 
near  or  above  1000K.  Rate  constants  for  reactions  of  similar  types  have  been 
equated  in  this  work;  however,  changes  in  rate  constants  depending  on 
molecular  size  may  be  quite  significant.  Molecular  and  ionic  processes  have 
been  ignored  in  the  present  study;  however,  there  is  as  yet  no  proof  against 
their  occurrence.  It  is  believed,  however,  that  due  to  the  reasonable 
agreement  between  the  present  model  and  experiments,  contributions  due  to  such 
reactions  are  perhaps  small. 


Conclusions 

A  chemical  kinetic  model,  revised  to  include  growth  of  aror itic  rings 
predicts  profiles  of  acetylene  decay  and  formation  of  benzene,  vinylacetylene, 
ethane,  and  hydrogen  which  are  in  agreement  with  experimental  flow  reactor 
results  near  1000K.  In  addition  the  model  predicts  the  formation  of  styrene, 
phenylacetylene ,  naphthalene  and  other  fused  rings.  An  analysis  of  the 
detailed  model  indicates  that  the  dominant  route  for  growth  from  a  single  to  a 
fused  ring  is  due  to  addition  of  phenyl  radicals  to  two  acetylenes.  Addition 
of  phenyl  to  vinylacetylene  was  proposed  and  may  play  a  significant  role 
depending  on  pressure  and  relative  concentrations.  Uncertainties  associated 
with  the  role  of  the  i-C,H,  radical  were  discussed  and  a  mechanism 

a  3 

involving  sequential  addition  of  i-C^H^  to  phenylacetylene  and  the 
resultant  products  was  proposed  as  a  conceptually  attractive  mechanism  for 
ring  growth.  Uncertainties  related  to  the  heat  of  formation  and  rate  of 
isomerization  to  n-C4H3  prevent  quantitative  predictions  as  to  the 
importance  of  such  a  mechanism. 
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TABLE  I 

REACTIONS  FOR  ACETYLENE  PYROLYSIS  NEAR  1000K 
logk  =  logA  +  nlogT  -  E/R/T/2.303  * 

Forward  Reverse 

Reactions  Rate  Constant  Rate  Constant 


logA 

n 

E 

logA 

n 

E 

1 

2C2H2=i-C4H3+H 

14.54 

0.0 

60.5 

13.49 

0.0 

0.2 

2 

H+C2H2-C2H3 

12.92 

0.0 

2.7 

12.96 

0.0 

41.6 

3 

C2H3+C2H2=n-C4H5 

12.88 

0.0 

8.0 

14.86 

0.0 

44.1 

4 

n-C4H5«C4H4+H 

13.00 

0.0 

33.0 

12.39 

0.0 

-5.4 

5 

n-C4H5+C2H2=£-C6H7 

12.18 

0.0 

5.0 

14.25 

0.0 

41.0 

6 

c-C6H7«=£-C6H7 

14.48 

0.0 

50.0 

11.36 

0.0 

0.4 

7 

C6H6+H-C-C6H7 

13.60 

0.0 

4.3 

13.12 

0.0 

24.6 

8 

C2H3+C4H4-C6H6+H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

9 

C4H4+n-C4H5-AlC2H3+H 

13.50 

0.0 

5.0 

0.00 

0.0 

0.0 

10 

A1C2H3+H=A1C2H2S+H2 

14.00 

0.0 

14.5 

13.09 

0.0 

19.7 

11 

C2H4+C6H5*A1C2H3+H 

11.57 

0.0 

2.1 

12.92 

0.0 

1.4 

12 

C6H6+H-C6H5+H2 

14.40 

0.0 

16.0 

12.39 

0.0 

9.5 

13 

C2H2+C6H5-A1C2H2 

12.00 

0.0 

4.0 

13.72 

0.0 

42.2 

14 

C2H3+C4H4=C2H4+n-C4H3 

13.00 

0.0 

10.0 

13.09 

0.0 

9.5 

15 

C4H4+C6H5*C6H6+n-C4H3 

12.48 

0.0 

0.0 

13.04 

0.0 

3.1 

16 

£-C6H5=n-C4H3+C2H2 

13.40 

0.0 

36.1 

11.11 

0.0 

0.0 

17 

C6H5=£-C6H5 

13.54 

0.0 

65.0 

10.25 

0.0 

1.7 

18 

C2H3t-C6H6=C2H4+C6H5 

13.48 

0.0 

13.0 

13.00 

0.0 

9.3 

19 

2C2H3-C4H6 

13.40 

0.0 

0.0 

17.15 

0.0 

104.9 

20 

C4H4«n-C4H3+H 

15.00 

0.0 

100.0 

12.91 

0.0 

-7.9 

21 

C6H6=C6H5+H 

15.70 

0.0 

107.9 

13.05 

0.0 

-3.2 

22 

n-C4H5+H=C4H4+H2 

13.00 

0.0 

0.0 

12.99 

0.0 

66.0 

23 

n-C4H5+H»C4H6 

13.00 

0.0 

0.0 

14.84 

0.0 

107.9 

24 

n-C4H5+n-C4H3»2C4H4 

12.70 

0.0 

0.0 

14.18 

0.0 

69.5 

25 

n-C4H5+C2H3=C4H4+C2H4 

12.00 

0.0 

0.0 

13.57 

0.0 

69.0 

26 

C2H3+H=C2H4 

12.85 

0.0 

0.0 

15.00 

0.0 

107.2 

27 

C6H5+C2H3-A1C2H3 

12.60 

0.0 

0.0 

16.11 

0.0 

106.5 

28 

C2H3+H-H2+C2H2 

13.00 

0.0 

0.0 

13.57 

0.0 

65.5 

29 

C2H4+H-C2H3+H2 

14.84 

0.0 

14.5 

13.30 

0.0 

11.7 

30 

C4H4+H-n-C4H3+H2 

13.90 

0.0 

14.5 

12.42 

0.0 

11.0 

31 

C4H6+H*n-C4H5+H2 

14.00 

0.0 

14.5 

12.76 

0.0 

10.9 

32 

C4H6+C2H3«n-C4H5+C2H4 

13.18 

0.0 

10.0 

13.52 

0.0 

9.5 

33 

C6H5+C6H6-C12H10+H 

11.80 

0.0 

11.0 

13.38 

0.0 

9.2 

34 

C2H4+M=C2H3+H+M 

16.16 

0.0 

81.8 

14.01 

0.0 

-25.4 

35 

21-C4H3-C4H4+C4H2 

11.00 

0.0 

0.0 

14.35 

0.0 

47.8 

36 

i-C4H3+H2«C2H2+C2H3 

10.70 

0.0 

20.0 

11.18 

0.0 

14.8 

37 

C4H4*i-C4H3+H 

15.20 

0.0 

95.0 

12.72 

0.0 

-2.1 

38 

C2H+C4H4-C2H2+i-C4H3 

13.60 

0.0 

0.0 

12.48 

0.0 

27.9 

39 

n-C4H3»C2H2+C2H 

14.30 

0.0 

57.0 

13.56 

0.0 

3.0 

40 

i-C4H3*C4H2+H 

12.00 

0.0 

49.0 

12.86 

0.0 

-0.2 

41 

n-C4H3-C4H2+H 

12.60 

0.0 

40.0 

13.04 

0.0 

1.4 

42 

n-C4H3»i-C4H3 

13.00 

0.0 

35.0 

12.58 

0.0 

45.7 

43 

1-C4H3+H-C4H2+H2 

13.00 

0.0 

0.0 

14.47 

0.0 

55.2 

44 

n-C4H3+H*C4H2+H2 

12.48 

0.0 

0.0 

13.53 

0.0 

65.9 

*  NOTES:  Units  for  A:  cc, moles, sec. ,  Units  for  E:  kcal/mole. 

■  represents  forward  and  reverse  directions  included  in  model, 
represents  forward  direction  only  included  in  model. 
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TABLE  I  (continued) 

REACTIONS  FOR  ACETYLENE  PYROLYSIS  NEAR  1000K 
logk  -  logA  +  nlogT  -  E/R/T/2.303  * 

Forward  Reverse 

Reactions  Rate  Constant  Rate  Constant 


logA 

n 

E 

logA 

n 

E 

45 

C4H4+H*i-C4H3+H2 

14.49 

0.0 

14.5 

12.62 

0.0 

21.9 

46 

C6H5+C4H4*C6H6+i-C4H3 

12.18 

0.0 

0.0 

12.31 

0.0 

13.8 

47 

C2H3+C4H4=C2H4+i-C4H3 

12.70 

0.0 

10.0 

12.36 

0.0 

20.2 

48 

C4H4-C2H+C2H3 

15.70 

0.0 

115.0 

12.93 

0.0 

-7.8 

49 

H2+M*2H+M 

12.35 

-.5 

92.5 

11.74 

-.5 

-11.9 

50 

C2H2+M-C2H+H+M 

16.62 

0.0 

107.0 

15.25 

0.0 

-17.9 

51 

C2H+H2-H+C2H2 

12.85 

0.0 

0.0 

13.60 

0.0 

20.5 

52 

C2H+C2H2*C4H2+H 

13.60 

0.0 

0.0 

14.78 

0.0 

15.4 

53 

C2H+C4H2-C6H2+H 

13.60 

0.0 

0.0 

14.97 

0.0 

15.1 

54 

C4H+C2H2=C6H2+H 

13.30 

0.0 

0.0 

14.91 

0.0 

8.1 

55 

C4H+H2*H+C4H2 

13.30 

0.0 

0.0 

14.30 

0.0 

13.5 

56 

C6H+H2-H+C6H2 

13.30 

0.0 

0.0 

14.53 

0.0 

10.3 

57 

C2H+C6H6«C6H5+C2H2 

13.30 

0.0 

0.0 

12.05 

0.0 

14.1 

58 

C4H+C6H6*=C6H5+C4H2 

13.30 

0.0 

0.0 

12.29 

0.0 

7.0 

59 

C2H3+C4H2-C4H4+C2H 

13.48 

0.0 

23.0 

13.70 

0.0 

5.4 

60 

2C6H5-C12H10 

12.48 

0.0 

0.0 

16.57 

0.0 

108.3 

61 

C2H+C4H4*C2H2+n-C4H3 

13.60 

0.0 

0.0 

12.91 

0.0 

17.2 

62 

A1C2H3+H-A1C2H2+H2 

14.50 

0.0 

14.5 

13.29 

0.0 

11.7 

63 

A1C2H2S-A1C2H+H 

13.00 

0.0 

45.0 

13.11 

0.0 

0.8 

64 

A1C2H2-A1C2H+H 

12.30 

0.0 

37.0 

12.71 

0.0 

0.8 

65 

A1C2H2+C2H2-A1C4H4 

12.88 

0.0 

8.0 

14.81 

0.0 

44.5 

66 

A1C4H4-A2H+H 

10.00 

0.0 

0.0 

13.32 

0.0 

18.5 

67 

C6H5+C4H4-A2H+H 

11.60 

0.0 

0.0 

17.14 

0.0 

56.4 

68 

A1C2H+C2H3-A1C2HP+C2H4 

13.70 

0.0 

13.0 

14.02 

0.0 

10.5 

69 

A1C2H+H-A1C2HP+H2 

14.40 

0.0 

16.0 

13.18 

0.0 

10.7 

70 

A1C2H+C2H-A1C2HP+C2H2 

13.30 

0.0 

0.0 

12.85 

0.0 

15.2 

71 

AlC2H+i-C4H3-A2C2H+H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

72 

A1C2H+C2H3-A2H+H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

73 

AlC2HP+C4H4=AlC2H+i-C4H3  12.18 

0.0 

0.0 

11.51 

0.0 

12.6 

74 

A1C2HP+C2H2-A1C2HV 

12.30 

0.0 

4.0 

13.70 

0.0 

41.8 

75 

A1C2HV-A2 

11.00 

0.0 

0.0 

14.14 

0.0 

53.4 

76 

A2C2H2-A2C2H+H 

12.30 

0.0 

37.0 

12.71 

0.0 

0.8 

77 

A2+C2H2-A2R5+H 

12.30 

0.0 

4.0 

0.00 

0.0 

0.0 

78 

A2C2H2+C2H2-A2C4H4 

12.88 

0.0 

8.0 

14.81 

0.0 

44.5 

79 

A2C4H4-A3H+H 

10.00 

0.0 

0.0 

13.32 

0.0 

20.7 

80 

A2H+H-A2+H2 

14.40 

0.0 

16.0 

12.88 

0.0 

10.7 

81 

A2H+H-A2P+H2 

14.40 

0.0 

16.0 

12.88 

0.0 

10.7 

82 

C2H3+A2H-C2H4+A2 

13.70 

0.0 

13.0 

13.72 

0.0 

10.5 

83 

C2H3+A2H-C2H4+A2P 

13.70 

0.0 

13.0 

13.72 

0.0 

10.5 

84 

A2P+C2H2-A2C2H2 

12.30 

0.0 

4.0 

13.70 

0.0 

41.8 

85 

A2C2H+H-A2C2HX+H2 

14.40 

0.0 

16.0 

13.49 

0.0 

10.7 

86 

A2C2H+C2H3-A2C2HX+C2H4 

13.70 

0.0 

13.0 

14.33 

0.0 

10.5 

87 

A2C2HX+C2H2-A3 

12.30 

0.0 

4.0 

0.00 

0.0 

0.0 

88 

A3+C2H2-A4H+H 

12.30 

0.0 

4.0 

0.00 

0.0 

0.0 

89 

A3H+H-A3+H2 

14.40 

0.0 

16.0 

13.18 

0.0 

10.7 

*  NOTES:  Units  for  A:  cc, moles, sec. ,  Units  for  E:  kcal/mole. 

■  represents  forward  and  reverse  directions  included  in  model. 
-  represents  forward  direction  only  included  in  model. 
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TABLE  II 

Selected  Thermodynamics  at  300K 
(from  Stein(6)) 


Species 

Identification  Heat  of  Formation 

(kcal/mole) 

Entropy 

(eu) 

C2H3 

vinyl 

68.4 

56.5 

A1C2H3 

styrene 

35.3 

82.6 

A1C2H2 

A1CHCH 

91.2 

85.2 

A1C2H2S 

A1CCH2 

83.2 

83.8 

A1C2H 

phenylacetylene 

75.2 

76.4 

A1C2HP 

A1C2H 

133.6 

79.0 

A1C2HV 

A1(C2H)CHCH 

146.6 

95.1 

A1C4H4 

A1CHCHCHCH 

104.7 

98.0 

A2 

1-naphthyl 

94.4 

83.6 

A2P 

2-naphthyl 

94.4 

83.6 

A2H 

naphthalene 

36.1 

79.7 

A2C2H 

2-naphthylacetylene  91.5 

90.9 

A2C2H2 

A2CHCH 

107.4 

99.7 

A2C2HX 

A2C2H 

149.9 

92.1 

A2C4H4 

A2CHCHCHCH 

121.0 

112.5 

A2R5 

acenaphthylene 

61.7 

87.2 

A3 

phenanthrenyl 

108.5 

96.8 

A3H 

phenanthrene 

50.1 

94.2 

A4H 

pyrene 

55.2 

96.5 
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THE  PYROLYSIS  OF  ACETYLENE  AND  VINYLACETYLENE  IN  A 
SINGLE-PULSE  SHOCK  TUBE 


MEREDITH  B.  COLKET,  III 
L  trued  Technologies  Research  Center 
Silver  Lane,  Mail  Stop  30 
East  Hartford,  CT  06108 

Acrtvlene  and  vinylacetvlene  have  been  pyrolyzed  in  a  single-pulse  shock  tube  for  the  tempera¬ 
ture  range  1 100  to  2400°K.  at  total  pressures  of  approximately  eight  atmospheres  and  for  dwell 
times  of  approximateli  700  microseconds  Initial  concentrations  of  the  hydrocarbon  in  argon 
ranged  from  about  100  ppm  to  41* .  Gas  samples  were  collected  and  analyzed  using  gas  chroma¬ 
tographs  for  hsdrogen.  and  C,  to  C,o-hvdrocarbons.  The  data  from  the  pyrolysis  of  acetylene 
exhibit  substantial  production  of  vinylacetvlene.  benzene,  and  phenylaceiylene,  but  agree  well 
with  a  detailed  chemical  kinetic  model.  Data  from  vinylacetvlene  pvrolvsis  and  thermochemical 
arguments  suggest  a  chain  mechanism  by  which  H  adds  to  vint  lacetylene  and  the  resultant 
adduct  decomposes  to  acetylene  and  a  vinyl  radical.  Rate  constants  for  the  reverse  steps  of  those 
occurring  during  vinylacetvlene  and  benzene  pvrolvsis  have  been  calculated  using  thermody¬ 
namics  and  forward  rate  constants.  These  rcserse  rate  constants  assist  in  describing  the  produc¬ 
tion  of  vinylacetvlene,  benzene,  and  phenvlacetvlene  during  acetylene  pyrolysis. 

Four  separate  reaction  mechanisms  for  the  initial  formation  of  aromatic  ringy  have  been 
identified.  The  relame  important e  of  each  step  depends  on  the  ambient  temperature  and 
relative  concentrations.  Overall  steps  can  be  written  as 

C-H,,  +  n-C<Hs  -*  C*H„  +  H 
C-H-  +  n-C,H,—  C„H, 

C,H,  *  n-C.H  —  CsHsC-H,  +  H 
C,H,  +  C,H,—  CsHs  +  H 


Introduction  cess  has  been  obtained  in  this  substantial  effort 

despite  the  lack  of  confirming  experimental 
Classman1  has  convincingly  argued  that  the  data  describing  profiles  of  intermediates  for 
chemistry  of  fuel  py  rolysis  plavs  an  important  confirmalion/supporl  of  the  proposed  model, 
role  in  sooting  diffusion  flames.  Some  pyrolysis  It  is  the  objective  of  the  present  work  to 
steps  may  be  rate-limiting  for  the  production  of  extend  these  earlier  works  in  order  to  help 
incipient  soot  panicles,  and  therefore  the  total  elucidate  pyrolysis  steps  and  chemical  mecha- 
amount  of  soot.2  Armed  with  this  information,  nisms  related  to  the  formation  and  break-up  of 
many  research  programs  have  begun  to  iden-  aromatic  rings.  A  single-pulse  shock  tube*  9  was 
tify  the  rate-limiting  processes.  For  example,  selected  for  this  study,  not  only  because  shock 
Bittner  and  Howard,  and  Bockhorn,  et  al.4  tubes  are  one  of  the  few  devices  capable  of 
have  provided  experimental  confirmation  of  a  generating  conditions  of  importance  to  soot 
variety  of  high  molecular  weight  polycyclic  formation  in  diffusion  flames,  but  also  to 
aromatics  present  in  rich  flames.  Cole,  el  al.5  extend  the  information  on  soot  production 
and  Weissman  and  Benson6  have  suggested  already  generated  through  optical  studies  in 
mechanisms  with  rates  consistent  with  thermo-  shock  tubes.10"15 

chemical  analy  sis  and  experimental  data  for  the  Preliminary  versions  of  the  modeling  work 
production  of  single-ring  aromatics.  Frenklach,  have  been  described  separately  for  acetylene14 
et  al.'  has  modeled  soot  formation  from  acety-  and  vinylacetylene.15  The  present  work  com- 
lene  pyrolysis  in  shock  tubes  using  a  detailed  bines  the  earlier  models,  including  a  recent  one 
chemical  mechanism.  The  mechanism  includes  for  benzene  pyrolysis,16  into  one  model  which 
production  of  aromatic  rings  and  continual  satisfactorily  describes  both  decomposition  (of 
growth  of  polycylic  aromatics.  Qualitative  sue-  the  parent  hydrocarbons)  and  product  forma- 
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lion  Requited  revision*  include  slight  adjust¬ 
ments  to  rate  constants  and  to  assumed  thermo- 
dtnamic  parameters 


Description  of  Facilities 


The  SPST  used  in  this  program  is  285  cm 
long  and  has  a  diameter  of  3.8  cm  (i  d  ).  The 
driver  is  88  cm  in  length  and  can  be  turned  bt 
shortening  its  length  in  3.8  cm  increments:  the 
driven  .section  is  197  cm  long.  An  11.7  liter 
"dump  tank"  is  located  in  the  driver  (lower 
pressure)  section  30  cm  downstream  of  the 
diaphragm.  Pressure  profiles  were  determined 
using  Kistler  pressure  transducers  located 
15.25  and  2.50  cm  from  the  end  wall.  Arrival 
times  were  measured  to  within  one  microsec¬ 
ond  using  digitized  pressure  traces.  Calculated 
quench  rates  are  typical!*  105  K/sec  or  higher  in 
the  rarefaction  wave.  Starting  pressures  prior 
to  filling  are  0.2g.  and  leak  rates  are  less  than 
ltc/min.  Post-shock  temperatures  were  calcu¬ 
lated  based  on  the  measured  incident  shock 
velocity  and  normal  shock  wave  equations. 

The  procedures  for  performing  an  experi¬ 
ment  are  similar  to  those  described  b\  Tsang.17 
except  for  an  automated  sampling  system.  The 
sample  is  collected  at  the  endwall  of  the  shock 
lube  using  0.0-15  inch  i.d.  tubing  heated  to  over 
85°C.  Approximately  30  milliseconds  after  the 
gas  has  been  shock  heated  and  cooled,  a 
solenoid  valve  opens  to  the  evacuated  sample 
cell  and  then  closes  after  300  milliseconds.  The 
sample  storage  vessel  is  all  stainless  steel  with  an 
internal  volume  of  25  cc. 

The  sampling  volume  is  directly  coupled  to  a 
low  volume  (<3cc).  heated  inlet  system  of  a 
Hewlett  Packard  5880  A  gas  chromatograph. 
Valves,  detectors  and  software  integration  rou¬ 
tines  as  described  previously  N  enable  this  system 
to  provide  automatic  quantitative  detection  of 
hydrogen  and  hydrocarbon  species  up  to 
Cip-hydrocarbons.  Based  on  repeated  injections 
of  calibrated  samples,  overall  accuracies  are 
estimated  to  be  three  percent.  Calibration  gases 
were  stored  in  stainless  steel  cylinders  with 
degreased  valves  and  were  heated  to  approxi¬ 
mately  60°C  prior  to  injection. 

Argon  ( 99.999 %  pure)  was  obtained  from 
Matheson  and  was  the  principal  diluent.  Com¬ 
pressed  acetylene,  also  from  Matheson,  con¬ 
tained  about  1  to  acetone  (added  for 
stability)  depending  on  bottle  conditions.  Acety¬ 
lene  was  purified  by  repeated  freezing  and 
thawing  at  liquid  N2  temperatures  and  retain¬ 
ing  only  the  middle  509t .  Final  samples  still 
contained  0.1  to  0.29?  acetone.  Vinyfacetylene 
was  obtained  from  Wiley  Organics  and  con- 
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lamed  an  unidentified  hydrocarbon  with  a 
concentration  of  approximately  8000  ppm. 


Model  Description 

Detailed  chemical  kinetic  calculations  have 
been  performed  using  CHEMK1N, 16  LSODE.19 
and  a  version  of  a  shock  tube  code  originally 
developed  by  Mitchell  and  Kee*°  but  modified 
to  include  quenching  effects  in  an  SPST. 
Quenching  rates  are  determined  using  experi¬ 
mental  pressure  traces,  assuming  adiabatic  ex¬ 
pansion,  and  using  the  equation 

dT  Ty-idP 
dt  ~  P  y  dt  ' 

Calculated  quenching  rates  vary  as  a  function 
of  shock  strength  and  time.  Initial  quenching 
rates  are  as  high  as  2xl06K/sec  for  shocks 
producing  initial  post-shock  temperatures  of 
2000°K,  but  only  259f  of  that  rate  for  shocks 
producing  reflected  shock  temperatures  near 
1200°K. 

The  chemical  kinetic  model  used  in  this  work 
is  reproduced  in  Table  1  and  is  based  on 
proposed  mechanisms  for  the  pyrolysis  of 
acetylene21  and  ethene.*2  The  reaction  set  in¬ 
cludes  the  identification  of  the  C4H3  isomers 
following  Frenklach,  et  a).'  The  radicals  are 
denoted  n-C<Hj  for  the  (normal)  isomer  with 
the  unpaired  electron  on  the  terminal  vinylic 
carbon  and  i-C*Hj  for  the  (iso)isomer  with  the 
radical  site  on  the  interior  carbon  atom.  n-C^H; 
(radical  on  end  carbon)  was  the  only  C4H5 
isomer  considered  in  this  work.  Also,  the 
mechanism  includes  the  forward  and  reverse 
processes  of  the  predominant  path  for  decom¬ 
position  of  benzene,,  i.e., 

CsH*  +  H  e  CeHj  +  Hz  (39) 
Q,H*  =t  f-CeHj  (56) 

f-C«H *  j=s  n-CjH,  +  C2Hz  (55) 

where  f-CsHj  has  been  suggested16  to  be  the 
l-hexyne-S,5-dien-6-yl  radical.  T  he  rate  con¬ 
stants  (R55.R56)  have  been  revised  slightly 
from  the  previous  work16  in  order  to  be 
consistent  with  the  reverse  processes. 

As  required,  thermodynamic  estimates  of 
some  species  were  made  using  group  additivity 
techniques. ,s  Otherwise,  data  from  readily 
available  sources  were  used.ms  Both  estimates 
of  rate  parameters  and  those  calculated  using 
equilibrium  constants  are  dependent  on  the 
selected  thermodynamics.  Of  particular  con¬ 
cern  during  the  present  research  is  the  appar¬ 
ent  uncertainly  in  the  heat  of  formation  of 
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TABLE  1 

Proposed  Sri  of  Reactions  and  Rale  Coefficients 


logk 

=  log  A 

■*  nlogT  - 

L  k'l 

2  303  * 

Reactions 

Forward 

Raie  Constant 
log  A  n  E 

Reserse 

Rate  Constant 
log  A  n  E 

Ref 

1 

C.2H4  a  M  =  C2H3  a  H  a  M 

16.16 

0.0 

81.8 

14  24 

0  0 

-23.0 

22 

2 

C2H4  a  H  -C2H3  a  H2 

14  84 

0(1 

14.5 

13  53 

0.0 

14.0 

21 

3 

H  ♦  C2H2-C2H3 

12.74 

0  0 

2  5 

13  01 

00 

43  8 

21 

4 

C2H3  f  H  =  H2  C2H2 

13.00 

0  0 

0  0 

13.34 

0.0 

63.2 

21 

5 

H2  +  M-2H  +  M 

12.35 

-.5 

92.5 

1 1.74 

-.5 

-11.9 

21 

6 

C2H2  a  M  =  C2H  +  H  +  M 

1662 

00 

107.0 

15.25 

0.0 

-17.9 

21 

7 

2C2H2«n-C4H3  +  H 

12.30 

0  0 

45.9 

1  1.67 

0.0 

-25  0 

21 

8 

i-C4H3  ♦  H2  =  C2H2  *  C2H3 

10  70 

00 

20  0 

11.41 

0.0 

17.1 

PW 

9 

C4H4-i-C4H3  +  H 

15  20 

0.0 

85.0 

12  72 

0.0 

-12.1 

PW 

10 

C4H4=n-C4H3  +  H 

15  00 

0.0 

100.0 

12.93 

0.0 

-7  7 

PW 

II 

C2H  -  C4H4  =  C2H2  -  i-C.4H3 

13.60 

0.0 

0  0 

12.48 

0.0 

27.9 

21 

12 

n-C4H3=C2H2  +  C2H 

14.30 

0.0 

57.0 

13.56 

0.0 

3  0 

PW 

13 

i-C4H3  =  C4H2  +  H 

12.00 

0.0 

49.0 

12.86 

0.0 

-0.2 

PW 

14 

n-C4H3  =  C4H2  a  H 

12.60 

0.0 

40.0 

13.04 

0.0 

1.4 

PW 

15 

n-C4H3  a  H«i-C4H3  a  H 

13.48 

0.0 

0.0 

13.06 

0.0 

10.7 

PW 

16 

i-C4H3  a  H  =  C4H2  a  H2 

13  00 

0.0 

0.0 

14.47 

0.0 

55.2 

PW 

17 

n-C4H3  a  H  =  C4H2  -  H2 

13.00 

0  0 

0.0 

14.05 

0.0 

65.9 

PW 

18 

C4H4  a  H  =  i-C4H3  a  H2 

14  49 

0.0 

14  5 

12.62 

0.0 

21.9 

26.PW 

19 

C4H4  a  H=n-C4H3  a-  H2 

13.90 

0.0 

14  5 

1245 

0.0 

11.2 

26.PW 

20 

C2H3 a  C2H2=n-C4H5 

12.04 

00 

4.0 

13.79 

0.0 

37.7 

PW 

21 

C2H  a  H2-H  a  C2H2 

12.85 

0.0 

0.0 

13.60 

0.0 

20.5 

22 

22 

C2H  *  C2H2-C4H2  a  H 

13.60 

0  0 

0.0 

14  78 

0.0 

15.4 

21 

23 

C4H2-C4H  a  H 

1489 

0(1 

120.0 

13.28 

0.0 

2.1 

22 

24 

C2H  *  C4H2=C6H2  a-  H 

13.60 

0  0 

0.0 

14  97 

0.0 

15.1 

21 

25 

C4H  a-  C2H2*C6H2  a-  H 

1 3.30 

0  0 

0.0 

14 .91 

0.0 

8.1 

21.PW 

26 

C6H2*C6H  a  H 

14.89 

0.0 

1200 

13.05 

0.0 

5.3 

22 

27 

C4H  a  H2=H  a-  C4H2 

13.30 

0.0 

0.0 

14.30 

0.0 

13.5 

22 

28 

C6H  a  H2-H  a  C6H2 

13.30 

00 

00 

14.53 

0.0 

10.3 

22 

29 

C2H  a-  H»C2  a  H2 

12  00 

0.0 

23.0 

11.55 

0.0 

3.8 

PW 

30 

C2H  a-  M  =  C2  +  H  a-  M 

16.67 

0.0 

124.0 

1561 

00 

0  4 

36 

SI 

C4H4  a  C6H5=C6H6  +  I1-C4H3 

12.00 

0.0 

0.0 

12.56 

0.0 

3.1 

PW 

32 

C4H4  a  C6H5-C6H6  a  i-C4H3 

12.00 

0.0 

0.0 

12.14 

0.0 

13.8 

PW 

33 

C2H  *  C6H6-C6H5  +  C2H2 

13.30 

0.0 

0.0 

12.05 

0.0 

14.1 

34 

34 

C4H  a  C6H6-C6H5  +  C4H2 

13.30 

0.0 

0,0 

12.29 

0.0 

7.0 

34 

35 

C2H3  a  C4H2-C4H4  a  C2H 

13.48 

0.0 

23.0 

IS  46 

0.0 

3.1 

PW* 

36 

C2H3  *  C4H4-C2H4  *  n-C4H3 

1 1.70 

0.0 

16.3 

11.55 

0.0 

13  4 

35.PW 

37 

C2H3  a  C4H4-C2H4  *  I-C4H3 

11.70 

0.0 

16.3 

11.13 

0.0 

24.1 

35, PW 

38 

C6H6-C6H5  a-  H 

16  18 

0.0 

107.9 

13.55 

0.0 

-3.0 

16 

39 

C6H6  a  H-C6H5  +  H2 

14.40 

00 

16.0 

12.39 

0.0 

9.5 

34 

40 

2C6H5-C12H10 

12.48 

0.0 

0.0 

16.57 

0.0 

108.3 

17 

41 

C6H5  +  C6H6-C12HI0  a  H 

11.80 

0.0 

11.0 

13.27 

0.0 

8.4 

37 

42 

n-C4HS  +  C6H3-CI0H8 

13.00 

0.0 

0.0 

0.00 

0.0 

0.0 

PW 

43 

C2H2 a  C6H5-C8H6  a  H 

12  00 

0.0 

4.0 

13.50 

0.0 

3.1 

S8.PW 

44 

C2H4 a  C6H5-C8H8  a  H 

11.57 

0.0 

2.1 

12.99 

0.0 

1.6 

38 

43 

C2H3  +  C4H4-C6H6  a  H 

11.60 

0.0 

0.0 

0.00 

0.0 

0.0 

PW 

46 

n-C4H5-C4H4  a  H 

14.00 

0.0 

41.4 

13.42 

0.0 

3.2 

6 

47 

n-C4H5  a  C2H2-/-C6H7 

12.81 

0.0 

9.0 

14.80 

0.0 

44.4 

6.PW 

48 

n-C4H5  a  H-C4H4  a  H2 

13.00 

0.0 

0.0 

13  03 

0.0 

66.2 

PW 

49 

C4H6  +  H«n-C4H5  a  H2 

14.00 

0.0 

14.5 

12.80 

0.0 

11.2 

26.PW 

50 

C6H6  a  C2H-C8H6  a  H 

12.00 

0.0 

0.0 

12.25 

0.0 

13.2 

PW 

51 

C4H4-C2H  a  C2H3 

16.00 

0.0 

1050 

13.46 

0.0 

-15.5 

PW 

CLSHING'Combustion  114-L 
Final  Page  No.  4 


HURON  VALLEY  GRAPHICS 

Record  No.  163  N  /F 2 


MASTER 
proof  Set 


ARTICLE  NUMBER  101 


TABLE  I  (Conti  nurd) 


Reactions 

Forw  ard 

Rate  Constant 
logA  n  E 

Reverse 

Rate  Constant 
bgA  n  E 

Ref 

52 

C4H4  +  n-C-tll  5-C8H8  +  H 

13.90 

00 

3.0 

0.00 

0.0 

0.0 

PW 

53 

•C8H8  ♦  H-C8H6  +  H  +  H2 

14.60 

0  0 

7.0 

0.00 

0.0 

0.0 

PW 

54 

C6H5 +  C2H3=C8H8 

13.00 

00 

0.0 

16.34 

0.0 

104.4 

PW 

55 

/-C6H5=n-C4H3  «■  C2H2 

14.00 

0.0 

36.0 

11.71 

0.0 

-0.1 

PW 

56 

C6H5«/-C6H5 

13  54 

0.0 

65  0 

1022 

0.0 

14 

PW 

5? 

C6H6 +  H-C-C6HT 

13.60 

0.0 

4.3 

13  12 

0.0 

24.6 

PW 

58 

C-C6H7-/-C6H7 

14.48 

0.0 

50.0 

11.36 

0.0 

0  4 

PW 

•NOTES:  Units  for  A:  cc, moles, see. 

Units  for  E  kcal/molc. 

represents  forward  and  reverse  directions  included  in  model 
"  represents  forward  direction  onls  included  in  model 
P\V  indicates  rate  evaluted  from  the  present  work. 


vinyl263'  and  related  (e.g.  n-C,Hj  and  n-C,H5) 
radicals.  Uncertainties  in  the  heats  of  formation 
and  entropies  of  these  radicals  directly  translate 
into  errors  in  proposed  rate  constants.  Ther- 
modvnamic  parameters  at  room  temperature 
for  species  included  in  this  work  are  presented 
in  Table  II. 

Results  and  Discussion 

Several  series  of  runs  of  the  pyrolysis  of 
acetylene  and  vinvlacctylene  have  been  com¬ 
pleted.  Each  series  represents  approximated  10 
to  15  separate  experiments  and  each  experi¬ 
ment  in  a  series  has  the  same  initial  concentra¬ 
tion  in  argon.  For  each  run,  a  chemical  analysis 
is  performed  of  final  product  distribution. 
Initial  post-shock  temperatures  for  runs  from  a 
given  series  typically  range  from  1100  to 
2400°K.  Detailed  chemical  kinetic  model  calcu¬ 
lations  have  been  performed  for  most  series 
and  the  numerical  results  for  selected  series  are 
compared  to  the  experimental  data  in  Figs. 
1.3-5.  Mass  balance  data  were  presented 
previouslv.1'1  A  detailed  discussion  of  the  ex¬ 
perimental  data  and  kinetic  mechanism  follows. 

Acetylene 

The  decomposition  of  3.791  acetylene  in 
argon  and  production  of  the  major  species,  i.e. 
C|Hj,  H}.  and  QH;,  agree  well  with  existing 
kinetic  models.  As  seen  in  Fig.  I,  vinvlacety- 
lene  (VA),  benrene  and  phenvlacetylene  are 
also  observed  and  may  play  a  critical  role  in 
the  growth  and  production  of  polycyclic  aro¬ 
matics.  The  bimolecular  reaction  2C?Hj  — * 
C^H,  is  thought  to  occur  at  low  temperatures. 


Extrapolation  of  low  temperature  rate  con¬ 
stants  provided  quantitative  agreement  with 
the  experimental  profile.  In  addition,  the 
activation  energy  for  this  process  (—40  kcal/ 
mole)  is  in  reasonable  agreement  with  the 
assumption  of  intermediate  formation  of  the 
C4FI,  diradical,  which  is  approximately  35  to 


TABLE  II 


Selected  Thermodynamics  at  500 K 


Species 

Heat  of  Formation 
(kcal/molc) 

Entraps 

C,jH|i. 

43.6 

93  7 

C,„H, 

360 

79.5 

C,H( 

78.3 

79.6 

97.4 

81.8 

C-CeH: 

49.9 

72.1 

C.H, 

19  9 

65.2 

C«Hj 

78.5 

69.4 

f-C*H5 

139.2 

80.0 

C«H, 

169.7 

71.1 

C«H 

233.2 

74.3 

C,H* 

26.1 

66.6 

n-C,H, 

82.5 

69.1 

C.H, 

69.4 

66.1 

n-C,Hj 

125.1 

68.7 

i-C,H, 

1 15.2 

71.6 

C,H; 

111.7 

59.9 

C,H 

179.0 

62.8 

CjH, 

12.6 

52.4 

C,H, 

65.7 

54.5 

C,H, 

54.2 

48.1 

C,H 

128.5 

49.6 

c, 

200.2 

47.7 
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43  Leal  more  energciic  than  two  molecules  of 
aceivlene>  Lnforiunatelv ,  this  mechanism 
does  not  provide  a  plausible,  parallel  path  lor 
the  formation  of  benzene,  vet  the  similaritv  in 
profiles  of  benzene  and  VA  strongly  suggest  a 
parallel  mechanism.  Stein"  suggested  an  alter¬ 
native  chain  process  initialed  by  H-a:r> m  addi¬ 
tion  to  acetylene  to  fprm  VA.  i.e  , 

C>H->(  +  Hi  —  C->Hs(  +  C«H  )  — 
n-C.|Hj(-H)  -* 

This  sequence  is  identical  to  the  first  steps  in 
the  ring  formation  process  proposed  by  Frenk- 
lach,  ct  aid  The  final  steps  of  the  proposed 
process  include 

C4H<(-H)  —  n-C4Hj(  +  C2H2)  — 
f-CfiHs  — *  phenyl 

followed  by  subsequent  formation  of  benzene 
or  phenylacetvlene. 

Detailed  chemical  kinetic  modeling  of  the 
present  experimental  data,  however,  strongly 
suggests  that  earls  (i.e.  low  temperature)  ben¬ 
zene  formation  arises  principalis  from  acety¬ 
lene  addition  to  the  normal-butadieny  I  radical, 
not  the  n-C*Hj  radical,  i.e., 

n-C4H5(  +  C2H2)  —  f-C«H;  -* 
c-C  6H,(-H)  -*  QH6 

where  f-Cs Hr  is  defined  to  be  the  1,3.5-hexa- 
trien-l-yl  radical. 

Preliminary  calculations  using  this  mecha¬ 
nism  have  also  been  performed  for  comparison 
to  the  flow  reactor  data  on  acetylene  pyrolysis 
obtained  by  Munson  and  Anderson . 11  Semi- 
quantitative  agreement  for  the  production  of 
benzene  and  VA  was  achieved  and  thus  the 
bimolecular  reaction  involving  intermediate 
formation  of  a  C4H4  diradical  is  not  necessary  to 
describe  lower  temperature  acetylene  pyrolysis. 

Above  1500°K,  cyclic  compounds  are  formed 
principalis  by  acetylene  addition  to  n-C4H3  in 
agreement  svith  the  proposal  by  Frcnklach,  et 
al.‘  Previously  it  was  assumed,  hosvever,  that 
n-C4Hj  was  formed  through  path  A,  svhereas 
(R-12)  was  the  dominant  route  under  the 
present  conditions 

n-C4Hj(-H)  —  QHd-H)  —  n-C4H,  A 
C2H2  +  C*H  -*  n-C*Hj  (-12) 

This  conclusion  arises  from  results  dis¬ 
cussed  later  in  this  paper  which  show  that  VA 
decomposition  principally  involves  H-atom  ad¬ 
dition  to  C4H4,  not  H-atom  abstraction.  As 
temperature  increases,  the  abstraction  route 
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will  become  more  significant,  however,  the 
concentration  of  C2H  will  also  increase  so 
(R-12)  still  remains  competitive.  Furthermore 
(R-12)  is  a  straightforward  process  for  pro¬ 
ducing  the  normal  radical  and,  in  Fig.  1, 
assists  in  adequately  describing  the  production 
of  benzene  and  phenylacetylene.  Demonstra¬ 
tion  of  the  relative  sources  of  n-C«Hj  can  be 
seen  in  Fig.  2,  where  reaction  sources  and 
sinks  are  plotted  as  a  function  of  time  for  two 
diflerent  temperatures. 

Figure  1  shows  that  the  model  overpredicts 
observed  benzene  and  phenylacetylene  concen¬ 
trations  above  1600‘K.  Coincidentally,  this  is 
the  same  temperature  at  which  a  significant 
deficit  in  recovered  mass  first  appears?4  Some 
phenylacetylene  and  benzene  above  1600*K  is 
probably  converted  to  higher  molecular  weight 
species  which  are  not  observed  and  not  ac¬ 
counted  for  bv  the  model. 


Vinylacetylene 

Mixtures  of  0.01,0.1 15  and  1.09?  VA  in  argon 
were  pyrolyzed  over  the  temperature  range 
1100  to  2500°K.  The  data  for  the  lowest  and 
highest  concentrations  are  presented  in  Figs.  3, 
4  and  5.  The  other  species  not  shown  in  these 
figures  were  observed  at  concentrations  less 
than  2^  of  the  parent.  These  species  include 
methane,  ethylene,  allene/methylacetylene,  sev¬ 
eral  unidentified  Cy  species,  C^H*.  toluene,  and 
a  Cg-hydrocarbon  (possibly  indene).  Data  on  the 
series  with  intermediate  concentration  have 
been  presented  previously.15  The  most  signifi¬ 
cant  information  is  that  at  low  temperatures  and 
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Fig.  I .  Experimental  and  Model  Results  for  Pyroly¬ 
sis  of  3.79E  Acetylene  in  Argon.  Dwell  Time  — 
700  x  10'*  sec.  Total  pressure  «■  8  atm. 
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Fic.  2.  Net  Formation  Rate  for  Production  'De¬ 
struction  of  the  n-C,H,  Radical  2a)  at  1600k,  2b)  at 
1800k. 

Contributing  reactions: 

R7  C:H;  +  C7H;=nC<H,  +  H 
R12  n-C,H,»C.H3  +  C,H 
RH  n-C<H5**C.|H;  +  H 
R19  C«H,  +  Hs±n-C,H,  +  Hj 
R55  l-C*H,jan-C,H,  +  C;H: 

Quenching  wave  arrives  at  700xl0‘e  sec.  Total 
pressure  “  8  atm.  Pvrolysis  of  3.7(?  acetylene. 

low  initial  concentrations,  acetvlene  is  the  pre- 
dominant  product,  while  diacetylene  (as  well  as 
Ho)  is  produced  at  levels  approximately  ten 
times  less  than  that  of  acetylene.  Aromatics  are 
produced  at  higher  temperatures  and  concen¬ 
trations,  but  not  in  quantities  sufficient  to 
encourage  rapid  growth  of  polvcydics  and 
severe  mass  imbalance.  RelativeK  little  data  have 
been  presented  previously  for  comparison  to 
this  SPST  data.  Lundgard's31  review  has  identi¬ 
fied  one  previous  work  which  overlaps  the 
present  temperature  range.  Yampol’skii,  et  a!.,! 
examined  VA  pyrolysis  between  1023  and 


MT1AL  POST-SHOCK  TEMPERATURE  (K) 


Fic.  3.  Experimental  and  model  results  for  100 
ppm  Vinylaceiylene  in  Argon.  Dwell  time  “ 
700  x  10'6  total  pressure  m  8  atm. 


MT1AI  POST-SHOCK  TEMPERATURE  (K) 


Fic.  4.  Experimental  and  Model  Results  for  Ali¬ 
phatic  Species  during  Pyrolysis  of  1%  Vinylaceiylene 
in  Argon.  Dwell  lime  -  700x  10'Sec,  total  pressure 
•>  8  atm. 
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Fig.  5.  Experimental  and  Model  Results  for  Aro¬ 
matic  Products  during  Pyrolysis  of  1%  Vinylacetylene 
in  Argon.  Dwell  time  “  700 x  10'*sec,  total  pressure 
-  8  atm. 
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1273°K  and  produced  qualitaii' el \  similar  r<- 
sulis  10  tin  SI’S  1  wink.  Tltf  main  Ion  tempera¬ 
ture  aliphatic  found  b\  Yampolskii  was  acelv- 
lent-  will)  traces  of  methane  and  eihene,  al¬ 
though  diacetvlene  was  not  deiected.  Benzene 
was  also  observed  previous!) .  as  well  as  a 
poll  mil  whose  precursor  mas  lie  sisrene  The 
overall  decomposition  of  VA  observed  bs  Yarn- 
polsku.et  al.,  k=  1 .6*  10"  exp(-52.800caLRT) 
sec* 1  is  approximately  five  times  higher  than  the 
SPST  results. 

Kinetic  modeling  requires  the  identification 
of  the  initiation  step  and  a  chain  process  which 
describes  the  predominant  foini.ition  of  acety¬ 
lene  with  minor  production  ol  diacetvlene, 
stirene,  benzene  and  phemlacetclene.  The 
three  possible  initiation  processes  are 

R9  C4H4  — •  i-C-iHs  ■+  H  -  98  heal  mole 

RIO  Ci Hi  -*  n-CiH,  +  H  -  1 1 1  heal  mole 

R51  C4Hi  — *  C2H  +  CoH,  —  1 22  kcabmole 

where  the  estimated  endoihermicities  mav  each 
be  in  error  b\  as  much  as  10  heal  mole  due  to 
uncertainties  in  heats  of  formation  of  the 
hvdrocarbon  radicals.  In  this  work,  the  princi¬ 
pal  initiation  process  was  assumed  to  be  Reac¬ 
tion  9,  which  required  a  rate  expression  of 
1 0 1 3  -  exp  (-42800T)  sec'1,  although  the  A- 
factor  seems  high  and  Ea(n.  low,  for  the  C-H 
bond  scission.  It  is  important  to  note  that  the 
model  results  were  relatively  insensitive  to  the 
absolute  magnitude  of  the  initiation  rate  or  its 
temperature  dependence.  A  factor  of  three 
change  in  the  initiation  rale  resulted  in  ap¬ 
proximately  a  20-309i  change  in  the  overall 
rate  of  decomposition  and  the  formation  of 
products.  The  explanation  for  this  phenome¬ 
non  is  that  the  decomposition  of  VA  in  the 
range  1200-1400°k  is  controlled  by  a  chain 
mechanism: 


REPLACE 


CiHissi-CiH,  ->•  H 

1 5.2 

83000 

H  +  CiHi»n-C4Hs 

ISO 

IS70. 

n-CiH|**CjHj  ♦  C,H, 

14  0 

43360 

*  C,H; 

ISO 

cc. mole. sec 

43710. 

calmole 

where  the  last  three  reactions  account  for  the 
majority  of  the  decomposition  of  the  parent 
and  the  predominance  of  acetylene  in  the 
low-temperature  products.  This  chain  is  not 
only  important  specifically  to  understanding 
the  pyrolysis  of  vinylacetylene;  but,  as  discussed 
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earlier  in  this  paper,  also  to  the  mechanism  of 
ring  formation  during  acetylene  pyrolysis.  Ex¬ 
planations  for  elimination  of  other  mechanisms 
are  presented  below  . 

Alternative  chain  mechanisms  for  decompo¬ 
sition  of  VA  would  include  H-atom  abstraction 
from  CiHi  resulting  in  CiCCH  :CH2. 
HCiCC:CH2,  or  HCiCCH:CH.  Formation  of 
the  first  of  these  can  be  neglected  from  thermo- 
dvnamic  considerations.  In  addition,  direct 
formation  of  two  acetylene  molecules  from  this 
structure  seems  unlikely  The  second  isomer, 
i-CiHj,  should  decompose  principally  into  dia- 
reiy  lent,  since  breakage  of  the  C.-C  bond  would 
have  a  high  activation  barrier  due  to  the 
formation  of  vinvlidene  as  an  intermediate. 
The  last  isomer.  n-C4H3,  may  decompose  into 
C4 Ho  +  H  or  C2H  +  C2H2.  Estimated  rates 

Log  k 

11)00°  K)  Log,oA  £ 

R12 

n-CiHjj^CjH  +  CjH-  4.72  14.3  57000. 

R14 

n-CiH,  ssCiH..  +  H  5.88  12.6  40000. 

suggest  that  Reaction  14,  i.e.  production  of 
C4Ho,  is  faster  by  a  factor  of  14  near  1300°K 
where  VA  decomposition  is  observed.  At 
higher  temperatures  the  two  rates  approach 
one  another.  These  relative  rate  estimates  are 
supported  by  experimental  decomposition  data 
of  benzene  pvrolysis  from  Kern,  et  al.55  Those 
researchers  found  that  the  initial  ratio  of 
products  C2H,/CiH,  at  1700°k,  is  approxi¬ 
mately  two  to  one.  Assuming  the  overall  decom¬ 
position  path 

UHsl-Hl^QHs^n-CiHs  +  QtH2. 

then  the  C2H2/C«H2  product  ratio  suggests  that 
the  branching  ratio  k|4/kl2  is  also  approximately 
two  to  one.  The  calculated  ratio  using  the  above 
rates  is  three  (3)  at  I700'K.  Thus,  if  n-C«Hj  is 
the  principal  intermediate  from  VA  decompo¬ 
sition,  then  the  initial  production  rate  of  C4H2 
would  be  similar  to  or  higher  than  that  of  C2H2 
during  VA  pyrolysis.  Such  a  prediction  is  not 
substantiated  by  the  SPST  experiments. 

H-atoms  may  add  to  V A  at  locations  other 
than  to  the  secondary  acetylenic  carbon.  Addi¬ 
tion  to  the  other  carbons  would  form  H2C: 
CCHCH2,  HC^CCHCHj,  or  HC-CCHjCHj. 
However,  each  of  these  radicals  would  be  ex¬ 
pected  to  reform  VA,  decompose  into  products 
other  than  two  acetylenes  or  involve  energetic 
intermediates.  Thus  it  appears  that  the  decom¬ 
position  chain  involving  n-CaHj  will  dominate. 
Detailed  chemical  kinetic  calculations  using  the 
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specific  rale  constants  in  Table  I  are  consistent 
with  this  analysis. 

Production  of  diacetvlene  is  described  by 

R  -  CgH.,  —  i-QH,  +  RH 
R  +  C4H4  —  n-C4H,  +  RH 

followed  by  decomposition  of  CyH,  into  dia¬ 
cetylene  plus  H-atoms.  R  mas  be  an  H-atom  or 
a  hydrocarbon  radical. 

As  shown  in  Fig  5.  reasonable  agreement  with 
experimental  profiles  of  aromatic  species  were 
obtained  using  the  reaction  sequence  of  Table  1. 
The  apparent  underprediction  of  VA  at  ele¬ 
vated  temperatures  is  presumably  an  experi¬ 
mental  problem  due  to  partial  gas  sampling  of 
boundary  layers  which  contain  the  unheated 
parent  hydrocarbon.  Above  1600°K,  phenyl  is 
formed  from  the  n-C.,Hj  +  G.-H?  recombination 
and  is  followed  b\  conversion  to  benzene  or  phe- 
nylacetylene.  At  lower  temperatures,  the  pre¬ 
dominate  formation  routes  of  aromatics  are 

R45  C4H4  +  QjHy  —  C6H: 

(1 ,3,5-hexatnen-3-vl) 

Q,H;  —  c-C6H: 

C-QH;  —  QH,,  +  H 

for  benzene  formation,  and 

R32C4H4  +  n-C4H5  —  CeHo 
CflHg  — ►  C-C$Hg 
c-CgHq  — *  CsHf  +  H 

to  produce  styrene.  PA  is  produced  principally  b\ 

R53C„H8  +  H  —  CgH;  +  Ho 
CeHr  —  CfH6  +  H 

near  I500'K  and  below  ,  but  at  higher  tempera¬ 
tures,  acetylene  addition  to  phenyl  also  occurs. 
Reactions  45,  52,  and  53  were  assumed  to  be 
overall  processes  and  nonreversible.  For  the 
conditions  (radical  concentrations  and  tempera¬ 
tures)  at  which  these  reactions  contribute,  pre¬ 
liminary’  calculations  confirmed  the  validity  of 
these  assumptions.  Despite  the  good  agreement 
with  the  experimental  data,  there  are  several  un¬ 
settling  features  of  these  proposals.  Reaction  45 
requires  a  1,3  or  1,4  H-atom  shift  prior  tocydi- 
zation  and  Reaction  47,  which  has  been  pro¬ 
posed  previously5  should  be  a  preferred  route 
for  aromatization.  However,  k4:  =»  10Mcc/mole- 
sec  would  be  required  to  explain  the  experimen¬ 
tal  data.  This  rate  is  significantly  higher  than  the 
rate  in  Table  I,  2*  10'  'cc/mole-sec,  required  by 
the  acetylene  modeling,  and  previous  determi¬ 
nations  [7.5x  JO10  (Ref.  5)and  2 .7*  IO"(Ref.  6) 
at  1300-K]. 


Another  concern  of  the  modeling  results  is 
that  the  value  of  kj?  required  to  match  the 
experimental  data  is  approximately  100  times 
higher  than  a  previous  determination  using 
data  from  a  low-  pressure,  premixed  butadiene 
flame.  Rate  determinations  from  theSPST  data 
are  subject  to  errors  due  to  boundary-layers, 
quenching  effects,  and  modeling  complexity. 
However,  the  flame  determination  is  expected 
to  be  a  lower  limit  since,  when  evaluating  the 
rate  constant,  Cole,  et  at.5  assumed  that  the 
n-CyH5  radical  was  the  dominant  C<Hs  species. 
Other  isomers,  particularly  H2C:CHCCH2. 
which  is  less  reactive  and  more  stable  by  3 
kcal/mole,  can  contribute  significantly  to  mass 
53. 

The  assumption  of  the  predominance  of 
n-CyHj  in  the  flame  work  would  only  partially 
explain  the  difference  between  the  two  evalua¬ 
tions  Uncertainties  in  heats  of  formation  and 
entropies  for  n-C«Hs  and/or  vinyl,  which  may 
affect  the  n-C4Hs  concentration  via 

C?Hj  +  CjHjrin-CyHs 

mas  also  help  to  alleviate  the  difference,  but 
probably  by  no  more  than  a  factor  of  five. 
Alternatively,  reactions  not  considered  in  the 
present  work  may  describe  production  of 
aromatics  or  significantly  enhance  the  concen¬ 
tration  of  n-CyHj. 


Conclusions 

Single-pulse  shock  tube  data  have  been 
obtained  for  the  pyrolyscs  of  acetylene  and 
vinylaceiylene.  The  data  have  been  used  to 
support  previous  proposals  for  acetylene  and 
benzene  pyrolysis  with  some  revisions.  A  de¬ 
composition  model  for  vinylaceiylene  involving 
H-atom  addition  to  vinylaceiylene  has  been 
proposed  in  order  to  explain  product  forma¬ 
tion.  Severe  constraints  imposed  by  simultane¬ 
ously  modeling  decomposition  of  each  of  the 
parent  hydrocarbons  as  well  as  five  or  more 
product  species  have  been  satisfied  with  the 
proposed  model.  Vinylaceiylene,  benzene,  and 
phenylacetylenc  are  observed  during  acetylene 
pyrolysis.  The  sequence 

CjHj  +  H*t  CjHj 
C,Hj  +  C?H,  j*  n-C«Hj 
n-QHjisCsH,  +  H 
C*H2  +  n-C«Hj  -*  f-C«H: 

( 1 ,3,5-hexatrien- 1  -yl) 
f-CeH,  -►  c-C,H, 
c-CeH,  -  C«H«  +  H 
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describes  the  low  temperature  formation  of 
vinv facets lene  and  oenzene  ( 1  1  ()<•- I400'K  i 
Re  verse  rates  calculated  from  thcrmodv  namics 
for  the  first  three  steps  describe  the  decomposi¬ 
tion  of  vinvlacetv lene  At  temperatures  above 
1500°k,  phenvl  is  formed  bv 

Gfl -  +  HS  GH  +  Hj 
GH.  +  C.«H  n-C(Hi 
GH 2  +  n-GH?  f-GHv 
( 1-hexv  ne-3,5-dien-6-yl) 
f-GHjS^phenvl 

and  is  followed  bv  formation  of  benzene  or 
phcnvlacetylene.  Reverse  rates  calculated  from 
thermodynamics  are  consistent  with  a  previous 
proposal  for  decomposition  of  phenvl  during 
benzene  pvrolvsis. 

Aromatic  formation  during  v  mv  facet  v  lene 
pvrolvsis  at  low  temperatures  is  described  bv 

C5H5  +  C4H<  —  GH; 

( 1 ,3.5-hexatrien-3-vl) 

GH:  —  c-GH- 
c-GH;  -  GH,  +  H 
and  n-GH;,  +  C,H,  —  <-GH, 
f-GEly  — *  C-CgHg 
c-C,H,  —  GH„  «■  H 

The  proposed  rate  constant  for  the  first  step  is 
low,  but  the  reaction  requires  an  H-atom  shift. 
The  required  rate  for  the  addition  of  the  normal 
butandienv  I  radical  to  vinv  lately  lene  is  approxi¬ 
mately  100  times  higher  than  a  previous  mea¬ 
surement.  Consequentlv ,  it  is  possible  that  calcu¬ 
lated  concentrations  of  n-GHs  are  inaccurate  or 
that  other  mechanisms  are  significant. 
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COMMENTS 


F.  Temps,  MPI  F.  Stromungsforschung. 

W’est- German).  Would  the  author 
please  comment  on  the  possible  role  of 
vinvlidcne  (singlet  or  triplet)  or  triplet 
acetylene  on  the  acetylene  pyrolysis? 

Author's  Reply  V'inylidene  appears  to  be  a  strong 
candidate  for  initiation  of  acetvlene,  since  it  lies  only 
about  44  kcaPnole  above  the  ground  state  of  acety¬ 
lene  and  since  the  intermediate  transition  state  is  only 
two  (2)  to  four  (4)  kcal/mole  above  vinylidene'.  A 
subsequent  step: 

H.C  -  C:  +  HCCH  —  H,CCCHCH  —  i-C,H,  +  H 


the  sequence  involving  vinylidene  is  still  about  13  kcal 
higher  than  the  activation  energy  of  the  bimolecular  in¬ 
itiator  (Reaction  7)  used  in  the  present  work. 

Vinylidene  could  also  play  a  role  if  the  above  inter¬ 
mediate,  HjCCCHCH,  undergoes  a  1,4  H-atom  shift 
to  produce  vinyiacetylene.  Although  this  route  was 
considered  in  the  present  work,  there  is  no  apparent 
parallel  pathway  to  benzene,  which  adequately  de¬ 
scribes  experimental  profiles.  It  is,  of  course,  possible 
that  this  mechanism  contributes  to  the  production  of 
vinyiacetylene.  A  possible  role  of  triplet  acetylene  was 
not  considered  in  the  above  analysis. 


is  endothermic  (overall)  by  about  15  kcal/mole  (or  69 
kcal/mole  above  two  acetylene  molecules).  The  overall 
endothermicityof  this  sequence  islowerthan  other  pro¬ 
posed  bimolecular  initiators  which  formCjH  +  C?H)Or 
n-C4Hj  +  H.  which  are  about  86  and  79  kcal  endother¬ 
mic,  respectively.  Unfortunately,  the  endothermicity  of 
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1.  Carrington,  T.,  Hubbard.  L.M..  Schaefer, 
H.F.,  and  Miller,  W.H.:  J.  Chem.  Phys.  50(9),  1 
(1984). 


CUSHING/Combustion  114-L 
Final  Page  No.  17 


Proof  Se^i 

- - - MS 


c-.-i 


.»] 


HURON  VALLEY  GRAPHICS 
Record  No.  549  N  ( F2 


yvvv 

’j-a*  , 
.  N.  V 

«\>v- 

s.%"  v 

ssa 


■ 

i 

i 


-  <WW W!  Vl*-ATiVUU«.1  IQ1VHL1  V  JtT  Af.  .V  Vt’VTXT^"rX.wXT  *'  '  X."  X.’ AT  V." 


5  £  $/e:n,  Chemira '  Kinehc '  />;r,  S’ottonnl  bureau  of 
Standards.  L’.S.A  How  do  \ou  interpret  the  unusualK 
high  rate  constant  for  stvrene  formation- 

Author  s  Hep! i  As  summarized  in  the  conclusions  of 
ihe  manuscript,  the  rale  expression  lc,s  -  7.9*  10” 
expl-30(HI  RTl  comole-sec.  is  about  a  factor  of  100 
loo  high,  w  hen  compared  to  rate  constants  of  compa¬ 
rable  addition  processes  The  interpretation  is  that 
either  other  mechanisms  dominate  stvrene  formation 
or  that  rate  coefficients  and  or  thermodsnamics  for 
several  of  the  reactions/spccics  are  in  error.  The  most 
likelv  alternative  mechanism  is  a  dimerization  pro¬ 
cess;  the  most  likels  errors  associated  with  rate 
coefficients  and  thermodynamics  probably  involve 
vinsl  and  butadiene!  radicals 

There  is  presently  insufficient  information  avail¬ 
able  to  ascertain  which  of  the  above  (or  if  all) 
contribute  to  the  discrepance  Consequently ,  the  data 
are  presented  as  is 


K  H-  Homann.  Technuche  Hoehirhult  Darmstadt  Meet 
Germany  Did  sou  consider  a  dimerization  step  as  the 
first  reaction  in  the  pyrolysis  of  C;H:  and  C<H4*  In  a 
recent  studs  of  C(H?  pyrolysis  at  100eC  we  found  that 
ihe  homogeneous  initiation  step  was  2  C,H;,  — »  C»H(. 
Addition  of  H  atoms  to  the  reaction  svstem  had  no 
effect  on  the  products'. 

1.  Homann  K.H.  and  Piooll  L’.V.t  Ber  Bun- 
senges  Phys.  Chem  1986.  in  press. 

Author’s  Reph  In  an  earls  analysis  of  this  work,  the 
dimerization  of  acetslene  was  considered  as  the 
predominant  mechanism  for  vinylacetvlene  forma¬ 
tion.  The  rate  constant  (k  “  3.6«  10"  exp(-44 100'RT) 
cc/mole-sec)  reported  be  Bradles  and  Kistiakowskv1 
was  found  to  describe  nicely  the  loss  temperature 
formation  of  C«H..  Coincidental!) ,  the  overall  activa¬ 
tion  energy  of  this  process  is  similar  to  the  endother- 
micits  for  the  formation  of  HCCHCHCH,  thediradi- 
cal  intermediate.  The  endothermicits  of  this  process 
was  estimated  to  be  35  to  43  kcal/mole  by  interpreting 
results  from  Kollmar.  et  al.!  Despite  the  good 
agreement  with  experiment  and  reasonable  mecha¬ 
nistic  explanations  (see  also  the  response  to  Prof.  F. 
Temps  regarding  the  role  of  vinylidene)  several 
problems  arose.  First,  the  Bradles  and  Kistiakowski 
data  were  shown’  to  be  misinterpreted,  leaving  doubt 
regarding  their  reported  rate  constants.  Secondly, 
there  appears  to  be  no  parallel  pathway  to  the 
formation  of  benzene  at  low  temperature  which 
describes  the  experimental  results.  In  support  of  the 
chain  mechanism  during  acetylene  pyrolysis  is  that 
production  of  C,H,  and  C«Ht  can  be  explained  using 
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the  proposed  chain  which  includes  rate  constants  very 
much  consistent  with  previously  proposed  rate  con¬ 
stants  and  thermochemistry.  Additionally,  many 
other  products  were  not  reported  due  to  their 
relatively  low  concentrations;  but  their  existence  is 
strongly  indicative  of  radical,  chain  mechanism.  Fur¬ 
ther  support  for  the  chain  mechanism  is  drawn  from 
the  fact  that  the  proposed  mechanism  qualitatively 
predicts  the  product  distributions  observed  during 
low  temperature  pyrolysis  of  acetylene*  (allowing  for 
a  small  change  in  the  initiation  rate  and  for  pressure 
dependence  of  addition  reactions). 

Finally,  Callear  and  Smith5  have  shown  convinc¬ 
ingly  that  processes  very  similar  to  those  described  in 
this  work  occur  at  room  temperature  when  H-atoms 
are  added  to  acetylene.  The  conclusion  of  the  present 
work  is  that  a  radical,  chain  mechanism  adequately 
describes  the  observed  dominant  products  using 
mechanisms  and  rates  which  are  reasonable.  Dimeri¬ 
zation  processes  cannot  as  yet  be  ruled  out  absolutely 
(and  they  may  play  a  contributing  role),  but  it  does 
not  seem  necessary  to  invoke  their  existence  to 
explain  the  observed  results. 

As  strong  a  case  against  dimerization  processes 
cannot  be  made  in  the  case  of  vinylacetylene,  since  it 
is  necessary  to  add  a  reaction  (#45)  which  had  not 
been  proposed  previously  and  (most  disconcerting)  to 
use  an  unusually  high  rate  for  the  addition  of  n-C«Hj 
to  C,H«  (R52)(see  response  to  question  by  Dr.  S. 
Stein).5  However,  with  the  lack  of  evidence  of 
dimerization  processes  occurring  during  high  tem¬ 
perature  acetylene  processes,  the  author  favors  a 
chain  mechanism  for  explaining  the  product  distribu¬ 
tion.  The  viny  lacetylene  data  can  be  fitted  reasonably 
if  one  assumes  dimerization  processes  dominate 
aromatic  formation.  For  example,  by  reducing  It*.  k», 
kj,.  and  k5J  by  factors  of  2,  2.  100,  and  2  respectively 
and  including: 

2C,H,  -•  C,H,(a)  and  2C«H,  —  C.H,  •+  2H  (b) 

(as  irreversible  steps)  with  respective  rate  constants 
of: 

k,  -  9x10”  exp(-40000/RT) 


It*  -  Jx  101*  exp(-35000/RT)  cc/mole-sec 

then  experimental  profiles  can  be  predicted  reasona¬ 
bly  (although  benzene  a  overpredicted  at  low  tem¬ 
peratures).  k.,  in  comparison,  is  approximately  two  to 
three  times  higher  than  an  extrapolation  of  Lund- 
gard  and  Heicklen's(6)  expression: 

k  -  3  x  |0"  exp(- 29200  cal/mole/RT)  cc/mole/sec 

for  the  net  formation  of  “dimerized"  products, 
formed  via  a  combination  of  molecular  and  radical 
processes.  The  radical  process  at  low  temperatures 
was  assumed  to  involve  formation  of  a  diradical 
intermediate,  although  possible  contributions  from 
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lilt’  chain  process  as  suggested  in  this  pa|>ot  should 
1101  tx  necessanh  ignored 
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• 

S  J  Harm.  General  Motors  Research  Labs,  L’.S.A.  Art 
C.;H;  and  C(H;  equilibrated  in  your  system?  Ai  what 
temperatures?  If  they  art  equilibrated,  then  you 
cannot  obtain  mechanistic  information  from  the 
contractions.  If  they  are  not  equilibrated,  are  thes  far 
enough  away  from  equilibrium  so  that  you  hate 
confidence  that  the  mechanism  is  controlling  the 
products? 

Author's  Reply.  For  the  dwell  times  (500-800  micro¬ 
seconds)  of  these  pyrolysis  experiments,  the  acety¬ 
lenes  are  equilibrated  only  above  1500-1700  K, 
depending  on  the  concentration  and  identits  of  the 
reactant.  Major  conclusions  of  the  present  work  are 
dependent  strongly  on  the  experimental  result  that 
products  of  Ion  temperature  (1 100-1300  K)  pyrolysis 
of  vinylacetylene  are  principally  acetylene,  diacety¬ 
lene.  and  hydrogen  with  the  product  ratios  C?H;' 
C,H:  and  ,  approximate!'  ten  to  one.  Pertur¬ 

bations  of  these  product  ratios  by  equilibration 
reactions  would  significantly  alter  the  conclusions  of 
this  paper.  The  "non-equilibration"  of  acetylene  at 
low  temperatures  and  short  limes  can  be  demon¬ 
strated  by  the  following  modeling  study  and  by 
comparison  to  experiments  on  benzene  pyrolysis 

In  an  artificial  modeling  exercise  of  C(H,  pyrolysis, 
rate  coefficients  for  Reactions  46  and  20  were 
decreased  while  coefficients  of  Reactions  13  and  18 
were  increased  by  factors  of  ten  ( 10)  or  more.  The  net 
result  was  to  convert  the  overall  decomposition 
pathway  from: 

C,H,  -  2CjH3 
to 

C,H,  —  C«Hj  +  Hj 

The  modeling  results  below  1400  K  indicated  final 
C«Hj/CsHj  ratios  of  approximately  one  hundred 
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(100)  io  one  which  is  29:  orders  of  magnitude  above 
equilibrium  Consequently,  this  modeling  exercise 
indicates  the  relative  unimportance  of  equilibrating 
reactions  under  the  low  temperature  conditions  at 
which  the  decomposition  mechanism  was  deduced. 

Experimental  evidence  for  non-equilibration  at 
low  temperatures  can  be  inferred  from  recent  data  on 
benzene  pyrolysis.  At  1500K,  several  hundred  de¬ 
grees  higher  than  the  temperatures  at  which  mecha¬ 
nistic  information  was  obtained  in  the  present  work, 

Colket1  found  an  initial  C?Hy,C,H.  ratio  of  approxi- 
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inateh  three  (3)  to  one  (I).  a  value  lower  than  the 
equilibrated  ratio  of  about  five  (5)  to  one  ( 1 )  (depend¬ 
ing  on  assumed  thermochemistry ).  The  experimental 
ratio  is  similar  to  the  value  (two  to  one)  obtained  by 
Kern,  et  al.2.  at  yet  higher  temperatures  (1704  K). 
Consequently,  both  modeling  and  experimental  data 
support  arguments  that  early  (i.e.  initial)  formation  of 
acetylenic  products  are  not  equilibrated. 

REFERENCES 

1.  Colket,  M  B  .  “Pyrolysis  of  C*Hk."  Preprints  of 
the  Division  of  Fuel  Chemistry.  American 
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2.  Ref.  33  of  manuscript. 

• 

J  H  Kiefer,  Unit  of  Illinois  at  Chicago.  L'.S.A.  In 
answer  to  Dr.  Harris'  question,  current  models  of 
high-temperature  acetylene  pvrolvsis  require  C»H 
radicals  for  C<H,  equilibration,  and  the  formation  of 
this  radical  seems  unlikely  at  these  temperatures.  I 
also  have  questions  for  the  speaker.  First,  for  your 
chain  production  of  CjH.  to  proceed  sufficiently  I 
suspect  you  may  well  need  the  high  rate  of  vinyl 
dissociation  you  assumed.  Otherwise  the  fast  process 
CjH>+H  — »  CjH;+H;  would  very  effectively  termi¬ 
nate  the  chain.  Do  you  think  you  could  stand  the  low 
rate  of  vinyl  dissociation  of  your  butadiene  work'? 
Finally,  if  the  heat  of  formation  of  C<Hj  were  about 
130  kcai/mol.  as  has  been  suggested  to  me  by  Karl 
Melius,  could  your  model  still  provide  sufficient 
conversion  of  C,H,? 

REFERENCE 

1.  Kiefer,  J.K.,  Wei,  H.C.,  Kern  R.D.,  and  Wu, 
C.H.:  Int.  J.  Chem. Kin.,  17,  225  <1985). 

Author's  Reply.  The  philosophy  for  selection  of 
reactions  and  corresponding  rate  constants  presented 
in  Table  1  was  to  use  well-established  mechanisms  and 
rate  constants  while  selecting  other  mechanisms  and 
rate  constants  as  required  by  the  experimental  re¬ 
sults.  but  consistent  (in  most  cases)  with  thermoche¬ 
mistry.  Vinyl  decomposition  was  assumed  to  be 
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pressure  independent  and  iis  rale  consiam  was  not 
varied  in  the  analysis  described  bs  the  manuscript  A 
preliminary  examination  ol  the  effect  of  the  use  of  a 
lower  rate  for  sinsl  decomposition  was  examined  bv 
using  Kiefer  s  rate  expression  w  hich  is  more  than  100 
times  slower  than  that  used  in  the  present  analysis  In 
the  modeling  of  vinvlacetvlene.  the  overall  decompo¬ 
sition  rate  of  the  parent  decreased  by  only  20(7. 
although  the  “quasi-steady"  concentration  of  vinyl 
radicals  increased  bs  nearls  a  factor  of  100  The 
overall  route  of  s  ms  lacetvlene  decomposition  changed 
onh  slightlv  from  that  described  in  the  manuscript. 
Main  problems  that  arise  are:  (!)  diacetylene  and 
ethene  are  overpredicted  by  about  a  factor  of  two  and 
five,  respectively;  and  (2)  earls  formation  of  both 
vinvlacetvlene  and  benzene  from  acetylene  py  rolysis 
are  underpredicied  each  bs  about  an  order  of 
magnitude.  Some  features  suggest  that  an  intermedi¬ 
ate  value  for  sins!  decomposition  might  be  reason¬ 
able.  The  data  on  vinylacetylene  and  acetylene  from 
this  work,  however,  do  not  suppori  the  low  rate 
constant  as  determined  in  your  analysis  of  ms 
butadiene  data. 

If  CiHt  radicals  have  a  heat  of  formation  of  130 
kcal/mole.  it  would  tend  to  support  present  arguments 
that  H-atoms  add  to  vinylacetylene  rather  than  ab¬ 
stract  H-atoms  from  s  ins  lacetylene  Bs  using  this  high 
heat  of  formation  in  modeling,  however.  two  major 
problems  will  arise.  First,  the  effective  activation 
energy  of  the  initiation  process  (C<H(  — >  i-CiH,  +  H) 
was  85  kcal'mole  in  this  work  This  low  value  is 
questionable  as  it  stands,  but  would  be  nearly  impossi¬ 
ble  to  reconcile  if  the  overall  endothermicin  were  1 13 
kcal'mole.  A  second  problem  results  from  the  fact  that 
the  high  heat  of  formation  (vs.  125  kcal/mole  assumed 
in  this  work)  would  substantially  reduce  the  “quasi- 
steady  state”  concentration  of  n-C|Hj.  Thus,  the 
formation  of  phenvl  via  n-C<H,  addition  to  acetylene 
would  be  substantially  reduced  and  benzene  and 
phenylacetylcne  significantly  underpredicted. 


W’.C.  Gardiner ,  Unhvrsity  of  Texas  L.S.A.  Inclusion 
of  irreversible  steps  in  polvmer  formation  mecha¬ 
nisms  has  to  be  regarded  with  special  care,  particu¬ 
larly  in  the  case  of  soot  formation,  where  the  diversity 
of  pathways  and  unsuspected  reversibility  of  quite 
exoergic  steps  make  the  points  of  no  return  first 
appear  at  molecular  sizes  well  beyond  those  con¬ 
sidered  in  this  paper*. 

In  order  to  make  valid  comparisons  between  com¬ 
puted  and  measured  yields  of  one-ring  and  two-ring 
aromatic  compounds,  the  reaction  mechanism  must  be 
extended — under  prov  ision  for  reverse  reaction — to 
still  larger  species,  as  a  minimum  to  accnapthalene 
For  this  reason  it  would  appear  that  the  rate  constants 
used  in  this  paper  to  model  formation  of  aromatic 
rings  have  to  be  regarded  as  artificial 
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and  Stein.  S.E.:  20th  Symp.  (Int.)  on  Combus¬ 
tion,  p  887.  The  Combustion  Institute,  Pitts¬ 
burgh,  1985 

Author t  Reply.  In  general,  your  comment  is  very 
reasonable  and  should  be  considered  for  long  reac¬ 
tion  times  or  conditions  which  allow  conversion  of 
significant  matter  to  larger  species  However,  for  the 
conditions  at  which  the  modeling  was  applied  in  this 
work,  this  effect  is  relatively  unimportant.  In  the 
present  study,  the  kinetics  analysis  was  performed 
for  short  reaction  times  (700  microseconds)  and  low 
temperatures  (<1600  K)  for  which  early  formation 
kinetics  (of  aromatics)  could  be  observed.  Under 
these  conditions,  relatively  small  fractions  of  the  ob¬ 
served  aromatics  undergo  growth  processes.  For  ex¬ 
ample,  in  the  case  of  3.5(7  acetylene  pyrolysis.  10 
ppm  of  benzene  was  observed  at  1 100  K.  the  lowest 
temperature  at  which  the  mixture  was  pyrolyzed. 
Concentrations  of  phenylacets  lene  at  this  tempera¬ 
ture  can  be  estimated  to  be  I ‘A  to  2  orders  of  magni¬ 
tude  lower  and  do  not  achieve  a  comparable  concen¬ 
tration  until  1500  K  (for  the  same  dwell  time).  Based 
on  data  using  higher  initial  concentrations  of  acety¬ 
lene  (4 .9*7? ).  naphthalene  concentrations  are  more 
than  an  order  of  magnitude  below  those  of  phenyla- 
cetvlene  at  these  temperatures.  An  analysis  of  mass 
balance  indicates  good  recovery  of  initial  carbon  be¬ 
low  1650  K.  These  results  may  not  be  proof,  but  thes 
indicate  strongly  that  the  processes  observed  in  this 
studs  are  descriptive  of  the  initial  formation  steps  (of 
aromatics)  and  perturbation  due  to  lack  of  inclusion 
of  subsequent  growth  processes  are  negligible. 

• 

Af.  Frenklach,  Pennsylvania  Slate  Untv.,  L.S.A .  Our 
mechanism  (1,2)  to  which  the  author  refers,  was 
composed  of  a  very  large  variety  of  reaction  pathways 
including  all  those  discussed  by  the  author.  No 
assumptions  about  the  importance  of  individual 
elementary  reactions  in  the  mechanism  were  made; 
our  conclusions  about  dominant  reaction  pathways 
leading  to  polycyclic  aromatic  hydrocarbons  and  soot 
were  drawn  instead  from  the  results  of  computer 
experiments  described  in  the  cited  references.  Thus, 
contrary  to  the  authors's  assertion,  the  pathway 

n  G,H,  C4H4  a  H 
C,H,  +  H  -*  n  -  C,Hj  +  H,  (1) 

was  not  assumed  but  deduced  to  be  the  dominant 
pathway  for  the  production  of  n  -  C<H)  for  the 
conditions  of  “soot”  formation  in  shock-tube  pyrolysis 
of  acetylene.  Although  we  found  the  reaction 

CtHj  +  C,H-*n  -  C«Hj  (2) 
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10  be  fast  at  vers  short  reaction  times,  it  quicklx  gave 
was  to  the  sequence  (1)  that  became  dominant  (or 
sikii  appearance  The  competition  between  (I)  and 
(2)  is  sensitive  to  experimental  conditions;  eg., 
reaction  (2t  plavs  a  more  prominent  role  in  the 
presence  of  oxx  gen  (3t  Hence,  differences  in  experi¬ 
mental  conditions  and  different  assumptions  about 
thermochemical  values  are  the  likely  sources  of  the 
different  relatise  fluxes  of  ill  and  (2)  noted  in  this 
paper. 
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Author's  Reph .  As  suggested  by  Dr  Frenklach.  the 
manuscript  should  be  changed  to  read:  "Previousls.it 
sx as  drdurtd  .  .  .  that  n-C4H3  was  formed  through  path 
A.  svhereas  (R-12)  was  the  dominant  route  under  the 
present  conditions.”  In  both  the  present  and  previous 
work,  it  should  be  remembered  that  the  authors  haxe 
been  forced  to  estimate — or  assume— values  of  abso¬ 
lute  rate  constants  lor  which  either  no  data  or  little 
data  are  available.  Oftentimes,  substantial  uncertaintv 
exists  in  estimation  techniques.  Consequently,  this 
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procedure,  bs  its  xers  nature,  leads  to  assumptions 
about  the  importance  of  indix-idual  elementary  reac¬ 
tions.  Conclusions  reached  or  deduced  in  any  work 
usually  are  only  as  good  as  the  initial  assumptions. 


A.Af.  Doan,  Exxon  Res.  tf  Eng.  Co.,  V.S.A.  Given  the 
fast  rate  constant  for  CjH5  a-  C,H,  (with  no  activation 
energy )  that  was  needed  to  explain  C*Ht  formation,  is 
it  not  possible  to  conclude  that  the  rapid  rate  of 
aromatics  formation  is  still  a  puzzle? 


Author's  Rrply.  The  rate  constant  for  L«s  was 
determined  to  be  4  x  10"  cc/mole-sec,  which  is  about 
two  orders  of  magnitude  below  collision  frequency 
and  is  similar  in  magnitude  to  rate  constants  deter¬ 
mined  for  other  addition  reactions.  The  only  unusual 
feature  regarding  Reaction  45  is  that  the  linear 
intermediate.  1 ,3.5-hexatrien-3-yl,  must  undergo  an 
H-atom  shift  and  convert  to  the  1,3.5-hexatrien-l-yl 
radical  prior  to  cydization.  It  is  believed,  however, 
that  this  H-atom  shift  is  not  rate  controlling.  The  size 
of  the  rate  constant,  in  the  opinion  of  the  authors 
therefore,  does  not  appear  to  be  unreasonably  large. 
The  activation  energs  could  not  be  determined  in  this 
experiment;  and  so  none  was  assigned,  although  it 
would  not  be  unreasonable  to  assume  an  effective 
activation  energy  of  at  least  several  kcal/mole. 

Certain  features  of  aromatics  formation  certainly 
can  be  considered  to  be  unresolved,  or  at  least  under 
discussion.  See.  for  example,  the  other  questions  and 
responses  to  this  paper.  Certainly  more  experimental 
and  modeling  work  needs  to  be  performed.  One 
critical  area  of  research  is  the  accurate  determination 
of  enthalpies  and  entropies  (and  their  temperature 
dependencies)  for  important  radical  intermediates. 
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SINGLE-PULSE  SHOCK  TUBE  EXAMINATION  OF  HYDROCARBON 
PYROLYSIS  AND  SOOT  FORMATION 


M.  B.  Colket ,  III 

United  Technologies  Research  Center 
East  Hartford,  CT  06108 


A  single -pulse  shock  tube  coupled  to  a  capillary  gas 
chromatograph  has  been  used  to  observe  the  pyrolysis 
of  a  variety  of  hydrocarbons.  Product  species,  mass 
balances,  and  H/C  ratios  of  uncollected  material  have 
been  determined.  Relative  "sooting"  tendencies  have 
been  found  to  be  toluene  >  benzene  >  allene  >  vinyl- 
acetylene  acetylene  >  1 , 3 -butadiene cyclohexene 
for  similar  carbon-atom  densities.  A  detailed 
chemical  kinetic  model  has  been  constructed  for  the 
pyrolysis  of  acetylene  and  describes  profiles  of 
major  species  as  well  as  the  production  of 
vinylacetylene ,  benzene  and  phenylacetylene . 


Introduction 


Shock  tubes  are  frequently  used  for  studying  kinetics  and  mechanisms 
of  both  hydrocarbon  pyrolysis  and  soot  formation  at  conditions 
characteristic  of  diffusion  flames.  Single-pulse  shock  tubes  (SPST) 
have  been  used  perhaps  most  successfully  for  measuring  rates  of  primary 
decomposition  processes  when  fractional  decomposition  of  the  parent 
hydrocarbon  is  small  .  Its  applications  to  the  study  of  soot  may  be 
limited  by  wall  quenching  and  finite  quenching  in  the  rarefaction  wave; 
but  the  potential  for  providing  detailed  species  information  with  the 
aid  of  capillary  gas  chromatography  has  motivated  the  present  research. 


The  ultimate  goal  of  the  experimental  work  is  to  provide  data  from 
which  mechanisms  of  hydrocarbon  pyrolysis  and  soot  formation  can  be 
obtained.  As  part  of  that  goal  the  relative  tendency  of  different 
hydrocarbons  to  soot  in  a  SPST  is  measured  in  this  work.  ^ In  addition, 
an  objective  of  this  paper  is  to  extend  a  recent  proposal  and  present 
some  results  related  to  soot  formation  during  the  pyrolysis  of 
acetylene.  Specifically,  experimental  and  modeling  results  on  the 
production  of  vinylacetylene  (VA) ,  benzene,  and  phenylacetylene  (PA)  are 
shown. 


Description  of  Facilities 


The  SPST  used  in  this  program  is  258  cm  long  and  has  a  diameter  of 
3.8  cm  (i.d.).  The  driver  is  88  cm  in  length  and  can  be  tuned  by 
shortening  its  length  in  3.8  cm  increments;  the  driven  section  is  197  cm 
long.  Pre-stressed  aluminum  diaphragms  were  used  to  control  the  rupture 
pressure.  An  11.7  liter  "dump  tank"  is  located  in  the  driver  (lower 
pressure)  section  30  cm  downstream  of  the  diaphragm.  The  pumping 
station  is  located  at  this  tee.  Pressure  profiles  were  determined  using 
Kistler  pressure  transducers  located  15.25  and  2.50  cm  from  the  end 
vail.  Arrival  times  were  determined  to  within  one  microsecond  using 


ss 
I'Si 


digitized  pressure  traces.  Calculated  quench  rates  are  typically  10 
K/sec  in  the  rarefaction  wave.  Starting  pressures  prior  to  filling  are 
0.2// and  leak  rates  are  less  than  l/Z^min.  Post-shock  temperatures  were 
calculated  based  on  the  measured  incident  shock  velocity  and  normal 
shock  wave  equations . 

The  procedure j  for  performing  an  experiment  are  similar  to  those 
described  by  Tsang  ,  except  for  an  automated  sampling  system.  The 
sample  is  collected  at  the  endwall  of  the  shock  tube  using  0.045  inch 
i.d.  tubing  heated  to  over  85  C.  Approximately  30  milliseconds  after 
the  gas  has  been  shock  heated  and  cooled,  a  solenoid  valve  opens  to  the 
evacuated  sample  cell  and  then  closes  after  300  milliseconds.  All 
valves  and  lines  in  this  system  have  been  degreased  and  rinsed  in 
toluene  to  minimize  condensation  on  the  walls.  The  sample  storage 
vessel  is  all  stainless  steel  with  an  internal  volume  of  25  cc.  Total 
internal  volume  of  the  remainder  of  the  sampling  system  is  less  than  3 
cc . 

The  sampling  volume  is  directly  coupled  to  a  low  volume,  heated 
inlet  system  of  a  Hewlett  Packard  5880  A  gas  chromatograph.  Typical 
injection  pressures  are  0.5  atmospheres  and  are  measured  to  within  two 
percent  using  a  calibrated  pressure  transducer.  One  gas  sampling  valve 
injects  samples  onto  a  25  m  x  0.33  mm  i.d.  CP  Sil  5  CB  (1  micron  film 
thickness)  fused  silica  capillary  column  (from  Chrompack  Inc.)  followed 
by  quantitative  hydrocarbon  analysis  using  a  flame  ionization  detector. 

A  second  valve  leads  to  a  silica  gel  packed,  stainless-steel  column  and 
a  single-filament,  modulated  thermal  conductivity  detector  for  analysis 
of  hydrogen.  Valves,  detectors  and  software  integration  routines  enable 
this  sytem  to  provide  automatic  quantitative  detection  of  hydrogen  and 
hydrocarbon  species  up  to  C^-hydrocarbons .  Based  on  repeated 
injections  of  calibrated  samples,  overall  measurement  accuracies  are 
estimated  to  be  approximately  three  percent.  Calibration  gases  were 
stored  in  stainless  steel  cylinders  with  degreased  valves  and  were 
heated  to  approximately  60C  prior  to  injection. 

Argon  (99.999%  pure)  was  obtained  from  Matheson  and  was  the 
principal  diluent.  Other  chemicals  were  at  least  99.5%  pure  according 
to  gas  chromatographic  analysis,  except  for  vinylacetylene ,  allene,  and 
acetylene.  VA  contained  approximately  8000  ppm  of  an  unidentified 
C,. -hydrocarbon.  Allene  was  obtained  from  several  different  vendors 
yet  was  typically  87  to  93%  pure  with  a  variety  of  other  C^-hydro- 
carbons  as  Impurities.  Bulb-to-bulb  distillation  enhanced'7 the  purity 
level  to  approximately  97%.  Similarly,  the  initial  concentration  of 
acetone  in  acetylene  was  reduced  from  approximately  2%  of  the  initial 
acetylene  to  about  1000  ppm. 

Results 

Toluene,  benzene,  cyclohexene,  1,3-butadiene,  VA,  allene,  and 
acetylene  have  been  pyrolyzed  in  a  bath  of  argon.  For  each  hydrocarbon, 
a  series  of  shock  tube  runs  were  performed.  Denslti^g  behind  the 
reflected  shock  wave  range  were  approximately  1.8x10  carbon 
atoms/cc.  Reflected  shock  temperatures  and  pressures  ranged  from 
approximately  1100  to  2500  K  and  5  to  8  atmospheres,  respectively. 

Dwell  times  were  typically  550  to  800  microseconds  although  for  the  bulk 
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of  runs,  dwell  times  were  close  to  700  microseconds.  Typical  results 
are  shown  in  Figs.  la  and  lb  for  the  pyrolysis  of  1.17%  benzene  in 
argon.  Presented  in  these  figures  are  the  final  product  distribution 
(for  18  separate  shocks)  after  the  test  gas  has  been  shock  heated  and 
then  quenched  by  the  rarefaction  wave.  The  general  features,  i.e., 
production  of  aromatics  over  the  1500  to  1900  K  temperature  range,  and 
the  dominance  of  hydrogen  and  'equilibrium'  concentrations  of  the 
acetylenes  above  2200  K  is  apparent  in  virtually  all  fuels.  For  the 
non-aromatic  fuels,  equilibrium  concentrations  are  generally  achieved 
above  1900  to  2000  K.  A;  temperatures  below  1500  to  1600  K,  the  product 
distribution  is  controlled  by  the  kinetics  of  the  decomposition  of  the 
hydrocarbon  parent . 

Figure  1  exhibits  a  relatively  large  amount  of  hydrogen  produced  at 
low  temperatures.  Modeling  demonstrates  that  this  'low  temperature' 
hydrogen  is  principally  formed  via  the  sequence: 


C  H  --->  C,H5  +  H 
CbH°  +  C  H  --->  C  H 

h6+^c.h  ^-- ->  c,h/I  ft? 
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A  comparable  amount  of  biphenyl  should  be  produced,  yet  is 
unobserved  due  to  its  apparent  loss  during  the  sampling  process  or  in 
the  inlet  system  to  the  gas  chromatograph.  The  amount  of  carbon 
converted  to  biphenyl  at  1400  K  would  be  approximately  a  few  percent  of 
the  total  and  is  consistent  with  the  measured  mass  balance.  Total 
hydrogen  and  carbon  mass  balances  have  been  determined  for  the  data  in 
Fig.  1  and  are  plotted  in  Fig.  2.  Also  shown  are  calculated  H/C  ratios 
of  the  missing  material.  (Uncertainties  in  this  calculation  are 
obviously  large  when  either  hydrogen  or  carbon  are  nearly  fully 
recovered.)  At  low  temperatures,  the  ratio  is  in  reasonable  agreement 
with  the  assumption  that  biphenyl  (H/C  -  0.825)  or  other  small 
multi-ringed  compounds  are  produced.  At  higher  temperatures,  the  H/C 
ratio  is  very  low,  suggesting  the  formation  of  very  large  polynuclear 
aromatics  (PAH)  or  even  soot.  The  low  ratio  of  0.1  observed  in  these 
benzene  experiments  is  atypical  since  for  most  other  hydrocarbons,  the 
H/C  ratio  of  the  missing  material  ranges  from  0.2  to  0.3  above  1700  K. 
This  difference  between  benzene  and  other  hydrocarbons  may  not  be 
significant  in  light  of  the  potential  errors  in  calculating  the  H/C 
ratios  of  the  missing  material. 


The  fractional  mass  recovered  is  a  strong  function  of  both  the 
concentration  and  the  density  of  the  parent  hydrocarbon.  The  effect  of 
hydrocarbon  type  on  mass  balance  is  shown  in  Fig.  3.  Assuming  that  the 
mass  loss  is  indicative  of  soot  productions,  the  relative  tendency  to 
soot  is  approximately  toluene  ,>  benzene  >  allene  >  VA  ^  acetylene  > 

1, 3 -butadiene  cyclohexene.  The  single-pulse  shock  tube  data  is  in 

excellent  agreement  with  the  relative  tendencies  to  soot  and  ^ith  the 
temperature  dependence  as  determined  using  optical  absorption  .  The 
detailed  species  profiles  offer  some  insight  to  the  relative  sooting 
tendencies.  Production  of  aromatics,  particularly  PA,  appears  to  be 
critical.  Allene,  for  example,  produces  large  quantities  of  benzene  and 
PA.  The  conversion  of  allene  to  benzene  is  as  high  as  25%  of  the 
initial  carbon  (and  4%  to  PA).  If  some  benzene  is  converted  to  PAH  and 
possibly  soot,  then  the  actual  production  of  benzene  may  be  even  larger. 
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Thus  the  "sooting"  tendency  of  allene  is  quite  similar  to  those  of  the 
aromatics.  In  comparison,  only  9%  of  the  initial  butadiene  carbon  is 
converted  to  benzene  and  the  production  of  PA  is  delayed  (relative  to 
its  production  during  benzene  or  allene  pyrolysis)  by  100  to  150  K.  The 
similarity  in  mass  balance  between  acetylene  and  VA  at  elevated 
temperatures  can  be  explained  by  the  dominant  overall  chain  process 
+  H  -->2C  +  H  apparent  from  the  product  distribution 

during  VA  pyrolysis.  Also  of  interest  are  significant  concentrations  of 
PA,  which  are  presumably  formed  via  C.H^  addition  to  VA,  followed  by 
cyclization  and  H-atom  loss.  The  production  of  PA  near  1400-1700  K  may 
.lead  to  the  significant  mass  deficit  observed  in  this  regime.  The 
similarity  between  cyclohexene  and  butadiene  is  explained  by  the  rapid 
unimolecular  decay  c -C^H^q- ->1 , 3-C^H,  4  C.H^  and  the 
relatively  low  sooting  tendency  for  fref.  3). 

Detailed  chemical  kinetic  modeling  has  been  performed  using  Chemkin 
and  a  modified  shock  tube  code  originally  developed  at  Sandia  National 
Laboratories.  The  modified  code  includes  the  capability  of  modeling  the 
quenching  effect  of  the  rarefaction  w^vg  in  an  SPST.  The  basic  ^inetic 
mechanism  is  similar  to  modifications  ’  of  a  proposed  mechanism 
and  as  shown  in  Fig.  4  describes  well  the  distribution  of  mono-,  di-, 
and  triacetylene.  The  production  of  VA,  benzene  and  PA  is  indicated 
below : 


Log  A 

E(  cal/mole) 

H  4  C  H  — >C  H 

12.74 

2500. 

C-H-  %  C-H. -  - ->NC  h 

13.00 

8000. 

♦  h4  5 

13.70 

41400. 

nc;h*  +  cJL--->c.h,  4  h 

12.60 

9000. 

nc‘hj  ♦  eh2.. ->c‘h‘ 

12.48 

5000. 

W*  c2V-->csV4  H 

12.32 

5600. 

A  in  units  of  cc,  moles,  sec. 


The  selection  of  this  reaction  sequence  and  rate  constants  has  been 
based  on  a  combination  of  £hg  gPST  experimental  results  and  pioneering 
work  at  other  laboratories*'  .  The  experimental  and  modeling 
results  compare  well  at  low  temperature,  but  suggest  some  revision  may 
be  required  at  temperatures  above  1800  K. 
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Fig.  la  SPST  Pyrolysis  of  1.17%  Benzene  at  6  atmospheres  and  for 
700  microseconds:  distribution  of  aliphatic  products. 
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lb  SPST  Pyrolysis  of  1.17*  Benzene  at  6  atmospheres  and  for 
700  microseconds:  distribution  of  aromatic  products. 
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Fig.  2  Fractional  recovery  of  initial  hydrogen,  carbon,  and  benzene 

during  pyrolysis  of  1.17*  Benzene.  Also  shown  is  the  calculated 
hydrogen-to-carbon  ratio  of  the  uncollected  mass. 
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Fig.  4  Experimental  and  model  results  for  the  pyrolysis  of  3.7%  acetylene. 
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The  addition  of  alkyl  hydrocarbons  to  aromatic  rings  is  believed  to  be 
important  in  the  growth  to  polyaromatic  hydrocarbons  in  diffusion  flames. 

To  examine  this  problem,  ve  have  pyrolyzed  1,3-butadiene  in  the  presence  of 
benzene  using  a  single-pulse  shock  tube.  Butadiene  was  selected,  since  the 
vinyl  radical  is  an  important  intermediate  during  its  pyrolysis  and  since 
1,3-butadiene  is  known  to  pyrolyze  at  temperatures  significantly  below 
those  at  which  benzene  pyrolyzes. 

Four  series  of  experiments  vere  performed  over  the  temperature  range  of 
approximately  1100  to  2000  K,  for  dwell  times  of  approximately  700  micro¬ 
seconds,  and  total  pressures  of  approximately  five  to  seven  atmospheres. 

The  concentrations  in  this  series  of  experiments  vere  (1)  0.1752  butadiene; 
(2)  0.175Z  butadiene,  1.05Z  benzene;  (3)  0.175Z  butadiene,  4. 61  benzene; 
and  (4)  1.1Z  benzene  vith  the  balance  argon.  The  results  of  the  first  of 
these  experiments  have  been  described  previously  (Colket,  1983).  In 
comparison  to  the  pure  butadiene  pyrolysis,  the  presence  of  benzene  does 
little  to  perturb  the  overall  pyrolysis  reactions  below  about  1450  K  except 
that  aignificant  concentrations  of  toluene,  styrene,  and  phenylacetylene 
are  formed.  Above  1450  K,  benzene  begins  to  decompose,  a  significant 
deficit  from  the  initial  mass  is  observed,  and  direct  comparisons  to  the 
pure  pyrolysis  is  questionable.  The  production  of  each  of  the  alkylated 
aromatics  in  mixtures  (2)  and  (3)  can  be  compared  vith  their  production 
during  the  pure  pyrolyses  of  butadiene  or  benzene  in  Figures  1-3.  These 
figures  clearly  show  significantly  enhanced  formation  of  alkylated 
aromatics  due  to  the  presence  of  benzene.  Types  of  reactions  vhich  may  be 
responsible  include: 


When  R  represents  CH.,  C.H.,  or  CjH,  the  respective  product  species  are 
toluene,  styrene,  and  phenylacetylene.  It  is  probable  that  reactions  (1) 
and  (3)  can  be  neglected  in  the  case  when  R  -  C,H  since  the  ratio  CgR/CgHg 
can  be  shown  to  be  very  lov  near  1400  X.  Alternatively,  some  of  the 
phenylacetylene  may  be  produced  from  the  decomposition  of  styrene. 


CONCENTRATION 
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A  surprising  experimental  result  is  that  an  increase  in  the  initial 
benzene  concentration  by  a  factor  of  about  four  does  not  result  in  a 
proportional  increase  in  the  concentrations  of  the  alkylated  aromatics. 

Very  little  increase  is  observed  in  the  case  of  the  toluene  or  pbenyl- 
acetylene  profiles  (below  1450  K)  and  atyrene  production  is  increased  only 
by  a  factor  of  1.5  to  2.  The  most  logical  explanation  for  this  result  is 
that  the  presence  of  benzene  must  bring  about  a  significant  decrease  in 
concentrations  of  the  alkyl  radicals,  furthermore,  type  2  reactions  are 
probably  not  dominant  since  RH  (i.e.,  CH^,  C^H^,  and  C^R^)  concentrations 
are  not  significantly  affected  by  the  presence  (or  concentration)  of 
benzene,  yet  an  increase  in  initial  benzene  concentrations  probably  results 
in  a  corresponding  (proportional)  increase  in  phenyl  concentrations. 

A  significant  decrease  in  alkyl  radical  concentrations  can  occur  if 
Reactions  1,  3  or  R  ♦  C.H.  — —  RB  ♦  C,l,(4)  are  fast.  Since  concen¬ 

trations  of  alkylated  aromatics  are  about  an  order  of  magnitude  less  than 
those  of  C.-  and  C--  hydrocarbona ,  Reactions  1  and  3  can  be  neglected 
relative  to  Reaction  4.  Ignoring  thermal  decomposition  and  recombination 
of  radicals.  Reaction  4  thus  should  be  compared  to  R  ♦  C.B;  —^■RH  ♦ 

(5). 

Reaction  4  should  effectively  compete  with  Reaction  5  since  the  C^H^/ 
C.H. ratios  for  mixture  numbers  (2)  and  (3)  are  about  6  and  25,  and  since 
Reaction  4  ia  only  about  2  keal/mdle  greater  than  Reaction  5  (for  the 
prefered  route  (Cole,  1982)  of  extraction  of  a  primary  hydrogen). 

Assuming  that  Reaction  4  ia  the  dominant  loss  mechanism  for  methyl  and 
tinyl  radicals.  Reaction  1  is  the  dominant  formation  route  for  toluene  and 
styrene,  and  that  phenylacetylene  is  formed  principally  from  styrene 
-decomposition,  then  ratios  for  k./k.  can  be  estimated  to  be  the  ratio  of 
final  concentrations  of  products f  i#ts9  For  both  sethyl 

and  vinyl,  kj/k^  vas  found  to  be  approximately  .3±.15. 

In  addition  some  detailed  chemical  kinetic  modeling  has  been  performed 
using  the  shock  tube  version  of  the  chemkin  code  (Mitchell  and  Kee,  1980) 
but  modified  to  account  for  the  quenching  processes  in  a  single-pulse  shock 
tube.  The  chemical  kinetic  mechanism  proposed  by  Kiefer,  et  al,  1984  for 
butadiene  kinetics  vas  used  except  that  tbe  thermal  decomposition  kinetics 
of  benzene  as  determined  by  fujii  and  Asaba  (1973)  were  employed.  It  vas. 
found  that  the  profiles  for  styrene  in  Figs.  1-3  could  be  approximately  fit 
regardless  of  which  reaction  mechanism  was  nsed.  The  determined  rate 
constants  for  the  vinyl  radical  when  Reactions  1,2,  and  3  were, each  need 
individually  are  k.<  3.5  x  10  ,  k.<  1.2  x  10* ,  and  k-<  1.8  x  10  cc/mode- 
aec  and  should  be  Considered  to  be  upper  limits.  ThJse  values,  however, 
are  critically  dependent  on  tbe  reaction  kinetics  of  the  phenyl  and  vinyl 
radicala.  Additional  modeling  efforts  are  in  progress  to  identify  which, 
if  any,  of  the  above  steps  dominate. 
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1  PRODUCTION  OF  TOLUENE 

00.17 %  ChH6,  4.6?  CgHg;  A  0.175?  C^Hg,  1.0?  CgHg 
00.17?  Cj,Hg;  V  1.1?  CgHg.  Dwell  times  are  approximately 
700  microseconds. 
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PYROLYSIS  OF  VINYLACETYLENE 
by 

M.  B.  Colket,  III 

United  Technologies  Research  Center 
East  Hartford,  CT  06108 

Eastern  Section  of  the  Combustion  Institute 
Philadephia,  Pa,  November  4-6,  1985 

Mixtures  of  0.01,  0.115,  and  1.0%  vinylacetylene  in  argon  have  been 
pyrolyzed  in  a  single-pulse  shock  tube  (SPST).  Reaction  conditions  varied 
from  approximately  1100  to  2500  K  and  6  to  8  atmospheres  and  dwell  times  were 
about  700  microseconds.  Gas  samples  were  collected  and  analyzed  using 
capillary  gas  chromatography.  Logarithims  of  the  final  (quenched) 
concentrations  of  stable  products  are  plotted  versus  initial  post-shock 
temperature  in  Figs.  1-4.  Principal  decomposition  products  at  low 
temperatures  and  low  initial  concentrations  are  acetylene,  diacetylene  and 
hydrogen,  although  acetylene  is  by  far  the  dominant  product.  (Hydrogen 
profiles,  although  not  shown  in  the  figures,  are  similar  to  or  slightly  higher 
than  the  profiles  for  C^^.)  With  increasing  temperature  tri-  and 
tetra-acetylenes  are  produced.  At  higher  concentrations  other  products  are 
observed  including  benzene,  styrene,  phenylacetylene  and  species  not  shown  in 
the  figures,  e.g.  methane,  ethylene,  methylacetylene ,  a  variety  of  C,. 
compounds,  C^H^,  toluene,  an  unidentified  Ch  compound,  and  naphtalene. 

The  carbon  balance  was  ±5%  for  the  series  or  100  and  1150  ppm  experiments , 
although  the  1150  ppm  series  did  not  extend  to  elevated  temperatures  where  the 
largest  mass  deficits  were  observed.  The  mass  deficit  for  the  1%  series  of 
experiments  was  20-25%  at  temperatures  above  2000  K.  The  "bell -shape"  which 
has  been  observed  for  aromatics  was  not  apparent.  Instead,  the  mass  deficit 
is  relatively  linear  over  the  temperature  range  1400-2000  K  and  is  similar  in 
magnitude  and  shape  to  that  observed  during  acetylene  pyrolysis  at  similar 
concentrations  (Ref.  1). 

A  detailed  chemical  kinetic  model  was  constructed  and  model  calculations 
were  obtained  using  a  shock  tube  code  (Ref.  2)  which  has  been  modified  at  UTRC 
to  include  the  effects  of  the  quenching  wave  in  the  SPST.  The  chemical 
kinetic  model  principally  consists  of  a  model  for  acetylene  pyrolysis  which 
has  been  based  on  previous  work  and  has  been  used  to  describe  recent  SPST 
results  (briefly  described  in  Ref.  1).  Several  reactions  specific  to  the 
vinylacetylene  system  were  added  to  this  acetylenic  mechanism.  There  are 
three  possible  Initiation  processes: 


(1)  C.H - M-C.H.  +  H  -98  kcal/mole 

(la)  C.H, - >n-C,Hf  +  H  -111  kcal/mole 

(lb)  C*H* - >C2H  +  C2H3  -122  kcal/mole 


where  the  estimated  endothermicities  may  each  be  in  error  by  as  much  as  10  to 
15  kcal/mole  due  to  uncertainties  in  heats  of  formation  of  the  hydrocarbon 


radicals.  In  this  work,  it  was  assumed  that, the  principal  initiation  process 
was  Reaction  1  with  a  rate  expression  of  101  '  exp(-42800/T)  sec* 


although  the  A-factor  seems  high  for  the  C-H  bond  scisson.  It  is  important  to 
note  that  the  model  results  were  relatively  insensitive  to  the  absolute 
magnitude  of  the  initiation  rate  or  its  temperature  dependence.  A  change  of  a 
factor  of  three  in  the  initiation  rate  resulted  in  approximately  a  20-30% 

G-} 


change  in  the  overall  rate  of  decomposition  and  the  formation  of  products. 

The  explanation  for  this  phenomena  is  that  the  decomposition  of  vinylacetylene 
in  the  range  1200-1400  K  is  controlled  by  a  Rice -Herzfeld  mechanisms: 


loS10A  E 


(1) 

initiation 

c,  H- 

- >i-C. H,  +  H 

15.2 

85000 

(2) 

chain 

H  +  C4H4' 

— >n-<B 

13.13 

1379 

(3) 

chain 

n-<H‘ 

--->C,H,  *  C,H0 

14.82 

45060 

(4) 

#  r  \ 

chain 

<X- 

13.0 

46000 

(5)  chain  termination 


where  Reactions  2-4  account  for  the  majority  of  the  decomposition  of  the 
parent  and  the  predominance  of  acetylene  in  the  low- temperature  products. 

(The  Arrhenius  factors  are  given  in  units  of  cc,  moles,  sec  and  cal/mole . ) 
Termination  is  rapidly  equilibrated  with  initiation  and  the  overall  rate  is 
dependent  at  most  on  the  square  root  of  the  initiation  rate.  Formation  of 
diacetylene  is  explained  by  the  chain: 

(6)  H  +  C  H - >i - C ,  H -  +  H.  14.5  14500. 

(7)  i - C,Ytl - >C.*C  iH  13.0  52500. 

4  3  4  2 

although  some  is  produced  via  initiation  (Reaction  1)  followed  by 

Reaction  7.  ^ 

Production  of  benzene  and  styrene  near  1400  K  is  described  by  Reactions  2 
and  3  followed  by: 


log10A  E 


(8) 

C.H,  +  C.H,---->C,H_ 

c-C6H7"">C6H6  *  H 

11.30 

0 

and 

(9) 

C4H4  +  n'C4H 5  — ->  CgH9 
C8H9--">c-C8H9 

13.6 

3000 

cC8H9  >C8H8  +  H 

Phenylacetylene 

is  formed  principally  by 

(10) 

CrH_  +  H---->C.H_  +  H, 

❖S-'XVe  S  * 

14.48 

7000 

at  low  temperatures,  although  with  Increasing  temperature  acetylene  addition 
to  phenyl  radicals  competes.  Reactions  8-10  were  assumed  to  be  overall 
processes  and  non- reversible  to  simplify  the  calculations.  The  value  of  kg 
appears  to  be  1  to  2  orders  of  magnitude  higher  than  expected  for  a  complex 
process  Involving  addition  and  cyclization.  This  rate,  however,  as  well  as 
many  others  in  this  system  of  reactions  is  highly  dependent  on  the  selected 
thermodynamics  and  competitive  rate  proceses  both  of  which  critically  affect 
the  quasi-steady  state  concentrations  of  the  radicals.  In  regards  to 
thermodynamics,  heats  of  formation  at  300  K  used  in  the  present  modeling 
effort  for  certain  species  are  reproduced  in  Table  1  and  are  literature 


values.  Of  note  is  that  the  heats  of  formation  of  diacetylene  and  vinyl -like 
radicals  are  in  general  disagreement  with  values  estimated  using  group 
additivity  techniques.  Of  specific  interest  is  the  necessity  of  using  a  high 
value  for  the  of  diacetylene  in  order  to  predict  properly  equilibrium 
values  of  diacetylene.  A  vaule  near  108  kcal/mole,  which  is  consistent  with 
group-additivity  techniques,  results  in  a  significant  over  prediction  of  the 
equilibrium  diacetylene  concentrations. 

The  comparison  of  the  experimental  data  and  the  model  (including 
chemistry,  thermodynamics  and  quenching  processes)  is  shown  in  Figs.  1-4  and 
exhibits  an  excellent  comparison  for  a  wide  range  of  species  inspite  of  a 
variation  of  two  orders  of  magnitude  in  the  initial  concentration  of 
vinylacetylene . 

TABLE  I 

Selected  Thermodynamics  at  300  K 

Heat  of  Formation  Entropy 

Species  (kcal/mole)  (cal/mole/k) 


C  H 

127.0 

49.6 

c9h 

63.4 

54.5 

C4H2 

111.7 

59.9 

i^H 

115.2 

71.6 

n'C4H3 

128.1 

67.7 

CA 

69.4 

66.1 

n-C^H 

82.5 

69.1 

C6H5 

78.5 

69.4 

C6H6 

19.9 

65.2 

C8H6 

C8H8 

78.3 

79.6 

35.2 

82.4 

Acknowledgments 

This  work  has  been  supported  by  the  Air  Force  Office  of  Scientific 
Research  (AFSC)  under  Contract  No.  F49620-85-C-0012 .  The  United  States 
Government  is  authorized  to  reproduce  and  distribute  reprints  for  governmental 
purposes  notwithstanding  any  copyright  notation  hereon. 

References 

1.  Colket,  M.B.,  "Single-Pulse  Shock  Tube  Examination  of  Hydrocarbon 
Pyrolysis  and  Soot  Formation"  in  Proceedings  of  the  Fifteenth 
International  Symposium  on  Shock  Waves  and  Shock  Tubes,  to  be 
published  1985. 

2.  Mitchell,  R.E.  and  Kee,  R.J.,  "A  General-Purpose  Computer  Code  for 
Predicting  Chemical  Kinetic  Behavior  Behind  Incident  and  Reflected 
Shocks"  SAND82-8205,  Sandia  National  Laboratories,  March,  1982. 


LOG  (F«*) 


n::u.  pcr-s-::> 


1  E**r  t  R-*UK£ !  * ; 


INiTIA..  POS'-SHOC*  TtMPtSfJti.P 


FIGURE  3  11  VI*UA:ETUC*£  PTRO:rs;$ 

till  vi.  ■oOf’,  70C  ric-ojfccnes,  iiipmucs 


FIGURE  4  It  FIRTlMCTUEKI  FTROtTSIS 
tall  *j.  most'  .  70C  »1Crot*COW>».  »ro«*t1ci 


100C  1100  1400  1400  1000  WOC  MOO  HOC 

INITIAL  FOST-SNOCK  TfWERATUREI*) 


G-4 


Pyrolysis  of 


by 

M.  B.  Colket,  III 

United  Technologies  Research  Center 
Silver  Lane,  East  Hartford,  CT  06108 


Introduction 

Several  mechanisms  for  the  pyrolysis  of  benzene  at  elevated  temperatures 
have  been  previously  proposed.  Unfortunately,  none  of  these  mechanisms  is 
entirely  satisfactory  when  both  experiments  and  thermodynamics  are  considered. 
Yet,  knowledge  of  benzene  decomposition  will  increase  the  understanding  of  the 
breakup  and  formation  of  other  aromatic  compounds.  Consequently,  a 
single-pulse  shock  tube  (SPST)  investigation  of  the  pyrolysis  of  benzene  has 
been  performed  over  the  temperature  range  of  1200  to  2400K.  In  addition, 
thermocheroical  estimates  and  detailed  chemical  kinetic  modeling  have  been 
performed  to  evaluate  the  previously  proposed  mechanisms  for  benzene 
pyrolysis . 


Description  of  Facilities 
and  Model 

The  3.8  cm  (i.d.)  single-pulse  shock  tube  (SPST)  used  in  this  experiment 
utilizes^the  "magic  hole"  technique  for  quenching  pyrolyzed  samples  at  rates 
above  10  K/sec.  SPSTs  were  developed  by  Glick,  Squire,  and  Hertzberg  (1) 
and  the  UTRC  facility  has  been  described  by  Colket  (2).  Gas  samples,  after 
dwell  times  of  approximately  700  microseconds,  were  automatically  collected 
and  analyzed  for  reactant  and  products  using  heated  gas  sampling  valves  and  a 
Hewlett  Packard  5880A  gas  chromatograph.  With  a  CP  Sil  5  CB  (from  Chrompack, 
Inc.)  capillary  column  and  a  silica  gel  packed  column,  H^  and  hydrocarbons 
up  to  C^q  were  identified  and  quantitatively  analyzed. 

Argon  (99.999%  pure)  was  obtained  from  Matheson  and  LC-grade  benzene  was 
obtained  from  the  Burdick  and  Jackson  Laboratory.  The  initial  mixture 
concentration  was  130  ppm  benzene  in  argon  and  was  prepared  gravimetrically . 
Gas  chromatographic  analysis  indicated  that  impurities  included  unidentified 
C^,  C^  and  Cy  hydrocarbons  as  well  as  toluene,  although  the  total 
concentration  of  impurity  was  less  than  0.2%  of  the  Initial  benzene. 

Detailed  chemical  kinetics  calculations  are  performed  using  CHEMKIN(3), 
LS0DE  (4),  and  a  version  of  a  shock  tube  code  (5)  which  has  been  modified  to 
Include  the  quenching  effects  in  a  SPST.  Quenching  rates  varied  with  shock 
strength  and  wera  calculated  using  measured  pressure  traces  and  assuming 
Isentroplc  expansion.  The  modified  code  also  allows  monitoring  of 
time -dependent  contributions  from  each  reaction  to  the  formation  and/or 
destruction  of  each  species. 


Proposed  Mechanisms 


Mechanisms  for  Che  decomposition  of  benzene  are  listed  in  Table  I.  It  is 
imporCant  Co  note  ChaC  most  are  not  single-step  mechanisms.  Radical 
Intermediates  equilibrate  rapidly  and  their  concentrations  are  sufficiently 
low  to  rei.der  them  (nearly)  undetectable  in  many  systems.  Experimental 
results  from  the  present  and  previous  high  temperature  works  (6,7)  suggest 
that  primary  stable  products  are  acetylene,  diacetylene,  and  hydrogen  and  the 
initial  production  rate  of  acetylene  is  two  to  three  times  higher  than  that  of 
diacetylene.  Unfortunately,  this  information  is  of  minimal  use  in  sorting  out 
the  mechanisms  since  benzene  pyrolyzes  at  high  temperatures  (relative  to 
pyrolysis  temperatures  of  other  hydrocarbons)  and  all  intermediates  shown  in 
Table  11  rapidly  pyrolyze  to  the  approximate  mixture  of  acetylene  and 
diacetylene  that  has  been  observed. 

Mechanism  A  is  the  generally  accepted  reaction  sequence;  however,  only 
Reaction  1  is  understood.  The  breaking  of  the  C-H  bond  is  believed  to  be  the 
initiation  step  in  both  pyrolytic  and  some  oxidation  studies.  Its  rate  has 
been  determined  from  D-atom  production  (in  pyrolysis  of  C,D,)  (8)  and  from 
detailed  modeling  studies  of  both  pyrolysis  (9)  and  oxidation  (10).  An  RRKM 
fit  has  been  performed  by  Kiefer,  et  al  (9)  using  available  experimental  data 
and  is  consistent  with^thermodynamics .  Their  kjf,  is  given  by  101  ‘ 
exp (-118  kcal/RT)  sec  .  Knowledge  of  the  mechanism  for  phenyl 
decomposition  is  substantially  less  than  that  for  Reaction  1.  Some 
information  is  available  on  the  overall  rate,  k^;  yet  the  details  of  the 
ring  fracturing  process  have  not  been  defined.  The  thermochemical  estimate  by 
Fujii  and  Asaba  (11)  has  been  until  recently  the  most  often  quoted  rate  for 
this  process.  More  recent  modeling  (8,9)  of  benzene  pyrolysis  has  produced 
similar  rates  as  Ref. 11;  however,  since  Reaction  1  is  rate  limiting  over  most 
of  the  temperatures  regimes  examined,  the  modeling  results  for  k^  are 
expected  to  be  lower  limit  estimates  with  large  uncertainties  in  "the 
teijgerature  dependence.  Thj  Fujii  and  Asaba  (11)  estimate,  k,  -  3.16  x 
10  exp  (-86  kcal/RT)  sec  ,  was  based  on  an  "old"  value  for'5 the  heat  of 
formation  for  i-C,H-  of  102  kcal/mole  (12);  however,  more  recent  estimates 
using  group  additivity  (13)  and  BAC-MP4  (14)  techniques  fix  this  value  closer 
to  115  kcal/mole.  In  addition,  n-C^Hj,  which  has  sAH^  of  126  kcal/mole, 
is  the  preferred  isomer  for  phenyl* decomposition.  Consequently,  Reaction  3  Is 
over  100  kcal/mole  endothermic.  If  this  process  involves  direct  CjH^ 
elimination,  then  the  previous  modeling  efforts  that  produced  lower  l^ig 
rates  strongly  suggest  that  the  high  pressure  A-factor  is  at  least  10  '  . 

This  value  is  orders  of  magnitude  higher  than  what  would  be  expected  for  this 
multibond  process.  Consequently,  it  must  be  assumed  that  Reaction  3 
represents  an  overall  process  or  that  processes  involving  radicals  other  than, 
or  in  addition  to,  phenyl  are  Important  to  benzene  decomposition. 

One  alternative  route  for  benzene  decomposition  is  direct  C„H- 
elimination  (Mechanism  B)  and  has  received  recent  support  (6,8)T  however, 
Kiefer,  et  al  (9)  have  shown  it  is  not  necessary  to  Invoke  this  step  to 
describe  experimental  profiles.  In  addition,  it  should  be  noted  that  the 
direct  elimination  is  a  multibond  process  which  should  be  unlikely,  and 
presumably  involves  the  intermediate  formation  of  the  C.H,  diradical.  For 
this  intermediate,  the  effective  activation  barrier  to  Reaction  4  may  be  as 
high  as  180  kcal/mole. 
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Another  possible  pyrolysis  route  involves  the  production  of  benzyne,  which 
subsequently  decomposes.  Unimolecular  decomposition  of  phenyl  (Reaction  5)  to 
benzyne  is  endothermic  by  93  kcal/mole,  which  is  nearly  as  much  as  Reaction  3. 
Similar  thermochemical  arguments  can  eliminate  this  decomposition  route.  The 
other  route,  Reaction  6,  is  a  radical  termination  step  and  would  produce  an 
overall  slowing  of  benzene  pyrolysis.  Knudsen  cell  pyrolysis  experiments  (7) 
have  shown  the  formation  of  a  C,H^  compound.  Assuming  this  compound  is 
benzyne  or  another  product  of  phenyl  decomposition,  Smith  and  Johnson  (7) 
argued  that  is  an  important  intermediate  during  benzene 

decomposition,  especially  at  elevated  temperatures.  At  least  some  of  the 
product,  however,  may  be  composed  of  the  chained  isomers  of  C,H, .  These 
isomers  may  be  produced  via  hydrogenation  of  triacetylene,  which  was  also 
observed  in  significant  concentrations.  Consequently,  no  clear  evidence  of 
the  Importance  of  Mechanism  C  is  apparent. 

Mechanism  D  appears  to  be  very  attractive,  since  an  extrapolation  of  the 
rate  coefficient  (15)  for  H-atom  addition  to  benzene  (Reaction  8)  is 
approximately  an  order  of  magnitude  higher  (Ref.  9)  than  H*atom  abstraction 
(Reaction  2)  at  1600  K.  The  addition  reaction,  however,  competes  with  its 
reverse  reaction.  Using  thermodynamics  for  c-C^H^  derived  from 
measurement  of  the  forward  and  reverse  rates  of°Reaction  8  (15),  Mechanism  D 
can  be  shown  to  have  a  negligible  impact  on  benzene  pyrolysis.  However,  a 
different  set  of  thermodynamics  for  c-C^H,  reported  in  Ref.  15  suggests 
that  Mechanism  0  may  play  an  Important  role  at  low  temperatures.  Detailed 
modeling  calculations  using  Mechanism  D  were  limited,  due  to  what  appears  to 
be  rather  large  uncertainties  in  both  the  heat  of  formation  and  entropy. 
Nevertheless,  there  are  attractive  features  of  this  decomposition  mechanism, 
especially  at  low  temperatures,  and  it  should  be  explored  further.  It  is 
worth  noting  that  Reaction  8  is  not  sufficiently  energetic  (only  16  to  26 
kcal) ,  that  it  can  be  followed  immediately  by  Reaction  9  (71-81  kcal 
required).  Instead,  c-C^H^  will  collisionally  thermallze  prior  to  its 
decomposition  to  products.  Only  a  minimal  acceleration  in  rate  due  to  the 
formation  of  an  excited  complex  can  be  expected.  Further  exploration  of  this 
route  (Mechanism  C)  should  be  performed. 

With  no  fully  satisfactory  alternatives,  and  the  expectation  that  phenyl 
must  decompose  to  aliphatics  at  sufficiently  elevated  temperatures  (when 
production  of  phenyl  via  Reaction  1  is  fast),  it  is  prudent  to  re-examine 
Mechanism  A.  One  can  separate  Reaction  3  into  the  following  sequence: 

C,H.~  WA  3a 

1-C^H5  -  n-C^H3+C2H2  3b 

This  sequence  is  the  reverse  of  the  processes  suggested  for  the  formation 
of  phenyl  during  acetylene  pyrolysis  (17)  and  has  been  shown  to  compare 
favorably  with  recent  experimental  data  (2).  If  one  assumes  that  the  ring 
breaks  at  the  single  bond  as  shown 
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then  this  is  a  relatively  straightforward  process  analogous  to  the 
decomposition  of  nC,H5;  i.e.  CH : CH . CH : CH.  •+  C.H.+C.H..  Breakage 
of  either  of  the  other  two  single 'bonds  would  require  a  more  complex  process 
involving  H-atom  shifting  or  the  formation  of  energetic  intermediates.  The 
resultant  linear  compound  would  be  expected  to  decompose  to  acetylene  and 
n-C^Hj  or  re-cyclize  to  phenyl.  Detailed  chemical  modeling  was  performed 
using  the  reaction  sequence  in  Table  II  coupled  with  an  acetylenic  mechanism 
similar  to  that  reported  (17,  16).  A  comparison  of  the  model  and  single -pulse 
shock  tube  data  Is  presented  in  Figure  1  for  a  series  of  shocks  at  an  initial 
'concentration  of  130  ppm  benzene  in  argon.  In  this  figure,  final 
concentrations  produced  after  a  dwell  time  of  700  microseconds  followed  by 
quenching  are  plotted  as  a  function  of  initial  post-shock  temperature.  Total 
pressure  for  these  experiments  is  approximately  seven  atmospheres. 

The  comparison  is  quite  reasonable.  The  higher  fractional  decomposition 
at  low  temperatures  observed  in  the  experiment  may  be  due  to  impurities  from 
the  wall  initiating  the  reaction  at  low  temperature.  The  residual  benzene 
observed  experimentally  at  elevated  temperatures  may  be  caused  by  sampling  a 
portion  of  the  boundary-layer  near  the  walls  of  the  shock  tube.  The  mechanism 
proposed  in  Table  II  Is  essentially  consistent  with  Kiefer  et  al’s  except  that 
Reaction  3  has  been  separated  to  3a  and  b.  An  effective  rate  constant  for 
Reaction  3  can  be  estimated  by  assuming  a  steady-state  concentration  for 
l-C,Hc,  i.e. 

6  5  ,[c6h5] 
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This  curved  evalution  is  depicted  as  a  solid  line  In  Fig.  2  over  the 
temperature  regime  where  sensitivity  to  this  rate  was  observed.  The  curve  is 
extrapolated  to  both  low  and  high  temperatures  to  facilitate  comparisons  to 
0te?fec¥?veat*°ns  ’  re^<,t^ve^y  high  values  obtained  In  this  work  for  kjand 

V ,  are  not  surprising  since  the  rate  constant  determined  in  these 

experiments  are  expected  to  be  close  to  the  high  pressure  limiting  value. 

Conclusions 


Previously  proposed  mechanisms  for  the  pyrolysis  of  benzene  at  high 
temperatures  have  been  reviewed  using  detailed  chemical  modeling  and  a  brief 
thermochemical  examination.  The  most  widely  accepted  sequence, 

C,H,  *■*  C,H,(+H)  -*  C.H.+C.H.,  is  slightly  modified  to 

C6H6"*  C6H5(+H>  ^  l*C6H5  n'C4H3+C2H2  *nd  the  *tructur* 

of  the  linear  C,H.  intermediate  Is  proposed.  Forward  and  reverse  rates 

are  consistent  with  thermochemistry  and  experimental  data.  A  benzene 


decomposition  route  involving  c-C^H^  and  1-C,H^  was  examined,  yet 
sufficiently  accurate  thermochemistry  was  not  available  to  reach 


aufflclently 

conclusions. 
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TABLE  I 

Previous  Mechanisms  for 
Benzene  Pyrolysis 

t 

i 

A 

C6H6  -  C6H5+H 

(1) 

• 

> 

H+ 

C6H6  -  C6H5+H2 

(2) 

c6h5  -  c4h3+c2h2 

O) 

! 

' 

X 

B 

C6H6  -  WC2H2 

(4) 

C 

C,H,  -  C,H,+H 

6  5  6  4 

(5) 

1 

H+ 

C6H5  C6  VH2 

(6) 

» 

l 

C6H4  -  C4H2+C2H2 

(7) 

1 

8 

\ 

5 

t 

D  H+ 

C6H6“*  c‘C6H7 
c-C6h7  -  1-C6H7 

(8) 

(9) 

!> 

l.C6H7^  n-C4H5+C2H2 

(10) 

i 

1 

i 

TABLE  11 

Proposed  Reaction  Sequence 
for  Benzene  Pyrolysis 

Forward  rate 
iogiQA  E 

cal/mole 

Reverse  Rate 

lo®10A  n  E 

cal/mole 

l. 

C6H6*  C6H5+  H 

16.18  107,900. 

10.05 

0.98 

-5690. 

i 

2. 

h+c6h6*c6h5+  r2 

14.40  16,000. 

8.35 

1.12 

6420. 

. 

3a. 

c6h5  ^i-c6h5 

14.00  65,000. 

13.11 

-0.68 

3300. 

3b. 

l-C6H5t  n-C4H3+C2H2 

15.34  38,000. 

5.97 

1.97 

-3610. 

’ 

11. 

n-C4H3  *H+C4H2 

12.43  37,000. 

11.29 

0.44 

-2790. 

12. 

n-C4H3  -  C2H+C2H2 

14.60  54,000. 

7.89 

1.66 

-3120. 

13. 

c2h+c6h6-*-c2h2+c6h5 

13.30  0 

12.00 

0 

12520. 
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FICURE  1.  1J0PPH  BENZENE  PYROLYSIS 
COMPARISON  OP  MODEL  AND  EXPERIMENT 


FIGURE  2.  ARRHENIUS  PLOT  OF  k3 
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